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1.0  INTRODUCTION 

1.1  Background 

A  need  exists  to  improve  the  procedures  employed  by  Coast  Guard 
personnel  to  position  aids  to  navigation.  Presently,  sophisticated  position¬ 
ing  systems  such  as  the  Global  Positioning  System,  and  Inertial  Positioning 
Systems  are  not  justified  for  widespread  Coast  Guard  use.  Until  such  systems 
are  justified,  alternative  procedures  must  be  explored.  Presented  here  are 
various  mathematical,  statistical  and  operational  procedures  which  will  be 
useful  in  efforts  to  position  aids  to  navigation  both  accurately  and  precisely. 

1.2  Discussion  of  the  Problem 


A  major  step  taken  by  the  Coast  Guard  to  improve  the  accuracy  in 
positions  of  aids  to  navigation  was  the  commitment  made  to  replace  graphical 
plotting  procedures  by  analytical  procedures  using  horizontal  control  coordi¬ 
nates  and  principles  of  geodesy  (reference  20).  In  this  commitment,  the  Coast 
Guard  displayed  an  intention  to  progress  from  procedures  based  on  common 
navigational  standards  to  procedures  based  on  the  more  stringent  geodetic 
survey  standards.  The  expression  used  herein  for  the  culmination  of  this 
progression  is  "hydrodetic  procedures1'. 

Graduation  to  hydrodetic  procedures  requires  effort  on  many  fronts. 
Among  these  efforts  are: 

a.  training  of  personnel 

b.  identification  of  error  sources 

c.  setting  standards  for  aid  positions 

d.  developing  procedures  to  meet  the  standards 

e.  developing  verification  and  recording  procedures 

Employment  of  hydrodetic  procedures  can  be  tedious  and  is  subject  to 
many  potential  error  sources.  When  computations  are  performed  in  the  field, 
this  drawback  boldly  presents  itself.  Proper  implementation  of  hydrodetic 
procedures  requires  planning,  operational  and  administrative  procedures.  Each 
is  discussed  in  more  detail  in  the  following  subsections  with  reference  to 
similar  procedures  that  are  employed  in  geodetic  survey. 

1.2.1  Planning  Procedures 

Implementation  of  hydrodetic  procedures  must  include  analyti¬ 
cal  planning  to  define  the  positioning  requirements  for  each  fixed  or  floating 
aid.  A  probability  expectation  approach  is  desired  to  allow  designation  of 
the  best  measurement  combinations  available  in  a  given  situation  and  the 
expected  result  when  using  any  one  of  the  combinations.  In  addition,  the 
planning  process  should  include  identification  of  the  situations  that  are 
prone  to  difficulty  and  find  remedies  if  possible.  The  planning  phase  is  of 
prime  importance  in  positioning  aids.  In  geodetic  survey,  such  phrases  as 
strength  of  figure,  standards  of  accuracy,  general  specifications,  and  recom¬ 
mended  spacing  are  associated  with  the  planning  task.  In  essence,  the  plan¬ 
ning  phase  is  based  on  the  determination  of  expected  accuracy  and  precision, 
including  propagation  of  errors.  If  the  proper  procedures  are  followed  with 
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the  expected  variability,  then  the  position  determined  should  meet  the  stan¬ 
dards  specified, 

1.2.2  Positioning  Procedures 


Following  planning,  the  tools  required  to  analytically 
position  an  aid  are  necessary.  The  adaptation  of  hydrodetic  procedures 
requires  a  firm  foundation  of  mathematical  models  and  statistical  analysis. 
Analytical  procedures  presently  being  tested  need  documentation  and  verifica¬ 
tion.  Guidelines  and  checks  for  handling  both  usual  and  unusual  situations 
are  needed,  as  are  out! ier  identification  procedures.  The  goal  is  to  provide 
a  consistent  set  of  equations  that  may  be  applied  at  various  levels  of  the  aid 
positioning  macro-structure  to  best  achieve  the  standards  established  during 
planning.  In  geodetic  survey,  such  terms  and  phrases  as  base  measurement, 
side  check,  number  of  readings,  rejection  limit  and  closure  are  associated 
with  the  positioning  task.  In  essence,  the  positioning  phase  is  based  on  a 
statistical  evaluation  of  the  accuracy  and  precision  in  a  set  of  measurements 
to  see  if  specified  standards  have  been  achieved;  and  if  not,  to  provide  a  set 
of  guidelines  to  achieve  them. 

1.2.3  Administrative  Procedures 

For  legal  and  future  planning  purposes,  the  data  accumulated 
during  any  positioning  evolution  must  be  verified  and  recorded.  Although  this 
report  does  not  focus  on  this  aspect  of  the  aid  positioning  problem,  it  is 
worthy  of  mention  because  planning  and  operational  procedures  are  dependent 
upon  the  constraints  imposed  by  administrative  procedures.  In  geodetic 
survey,  the  terms  and  phrases  adjustment,  geodetic  data  publications,  horizon¬ 
tal  control  data  and  published  data  are  associated  with  the  administrative 
task .  In  essence,  the  administrative  phase  is  one  of  data  accumulation  for 
future  use. 

1 .3  Summary  of  Previous  Work 

The  capability  for  using  hydrodetic  procedures  has  to  some  extent 
been  extended  to  all  Coast  Guard  aid  servicing  units.  The  Aids  to  Navigation 
Manual  -  Positioning  (CG-222-5)  has  been  revised  to  include  many  procedures 
that  can  be  called  hydrodetic  procedures.  Specifically,  creation  of  gradient 
diagrams,  procedures  for  selection  of  landmarks,  and  collection  of  horizontal 
control  data  are  all  in  progress.  These  procedures,  however,  are  still  far 
inferior  to  those  employed  in  geodetic  survey.  There  is  no  comprehensive  set 
of  standards  and  no  analytical  capability  to  measure  the  adequacy  of  a  posi¬ 
tioning  effort. 

Select  ships  are  involved  in  providing  a  data  base  for  use  in 
determining  the  degree  of  improvement  in  accuracy  and  precision  achieved  using 
selected  hydrodetic  procedures.  This  effort  is  documented  within  work  unit 
2702.2.2.5  Calculator-Assisted  Procedures  (CAP)  of  the  ANPAR  Project  Area 
(reference  T}"!  A  secondary  objective  of  CAP  is  to  create  a  data  base  for 
further  study  of  aid  positioning.  Feedback  collected  can  be  an  assist  in 
making  an  analytical  positioning  tool  which  is  easy  to  use. 
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Past  work  reported  within  work  unit  2702.1.2  Systems  of  the  ANPAR 
Project  Area  has  been  directed  toward  literature  research,  error  modeling,  and 

the  mathematical  concepts,  all  of  which  form  the  basis  for  this  report 

(references  1  and  2). 

1.4  Symbology,  Terminology,  and  Conventions 

A  list  of  symbols  has  been  provided  and  each  symbol  is  defined  where 
first  encountered.  When  possible,  the  symbols  have  been  selected  to  agree 
with  previously  defined  symbology.  In  some  cases,  this  leads  to  one  symbol 
representing  more  than  one  quantity  at  different  locations  in  the  report.  If 

confusion  occurs  as  to  the  meaning  of  a  specific  symbol,  refer  to  the  list  of 

symbols  which  designates  the  appropriate  meaning  for  a  specific  section. 

Terminology  is  in  accordance  with  the  Aids  to  Navigation  Manual 
(CG-222-5)  where  possible.  Terms  that  might  cause  confusion  are  underlined 
and  defined  in  a  glossary.  Terms  that  are  defined  in  Definition  of  Terms  used 
in  Geodetic  and  Other  Surveys,  (S.P.  #242),  U.S.  Dept  of  Commerce,  are  defined 
by  paraphrasing.  An  important  convention  used  herein  is  that  all  bearings , 
azimuths  and  directions  are  taken  clockwise  from  true  North. 
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2.0  APPROACH  -  SCOPE 


One  objective  of  this  report  is  to  develop  a  comprehensive  description 
for  hydrodetic  surveys.  By  its  nature,  a  complete  description  of  a  hydrodetic 
survey  requires  frequent  use  of  mathematics  beyond  the  level  of  the  average 
Coast  Guardsman.  This  imposes  serious  constraints  on  a  logical  presentation 
of  a  hydrodetic  survey  from  its  rudimentary  stages  to  rigorous  mathematical 
detail.  The  procedures  must  be  immediately  accessible  to  field  personnel  but, 
simultaneously,  the  procedures  must  be  supported  by  experimental  evidence  or 
mathematical  derivation.  In  order  to  meet  the  opposing  requirements,  all 
mathematics  beyond  basic  algebra  have  been  located  within  the  appendices  in 
the  same  order  as  in  the  body  of  the  report.  Frequent  reference  to  the  glos¬ 
sary  is  helpful  in  understanding  mathematical  terminology.  The  following 
subsections  are  general  descriptions  of  the  major  topics  discussed  in  the  body 
and  appendices  of  this  report. 

2.1  Geodetic  Survey 

A  discussion  of  geodetic  survey,  as  related  to  the  overall  goals  of 
hydrodetic  survey,  is  introduced  so  that  the  reader  may  draw  a  parallel  be¬ 
tween  hydrodetic  and  geodetic  survey  throughout  the  remainder  of  the  report. 

It  is  understood  that  some  of  the  geodetic  survey  terminology  is  introduced 
prematurely,  but  the  meanings  will  become  apparent  as  the  reader  progresses 
through  the  parallel  hydrodetic  survey  sections  that  follow  their  introduc¬ 
tion.  Following  geodetic  survey  is  a  description  of  the  reference  system 
employed  in  geodetic  calculation  of  distances  and  directions.  Some  approxima¬ 
tions  are  made  and  where  appropriate  the  accuracy  of  each  approximation  is 
discussed.  The  distances  and  directions  are  necessary  to  both  the  planning 
and  the  positioning  sections  of  this  report. 

2.2  PI anninq 

The  procedures  necessary  for  planning  a  positioning  effort  are 
introduced.  These  procedures  include  determination  of  expected  results  and 
suggested  procedures  to  achieve  the  desired  results.  The  procedures  are 
necessary  to  define  the  measurements  -  and  the  corresponding  lines  of  position 
(LOPs)  -  that  might  be  available  to  the  positioning  team  when  they  undertake 
the  positioning  effort.  Where  possible,  the  procedures  are  developed  so  that 
all  user  parameters  are  in  dimensions  recognizable  to  the  average  positioning 
team  member.  That  is,  the  standards  to  be  achieved  and  the  measurements 
needed  to  achieve  them  are  presented  in  a  specified  measurement  space  for  the 
user. 


2.3  Positioning 

Equations  to  compensate  for  systematic  error  sources  that  are  known 
to  the  user  are  presented.  Equations  are  derived  for  the  transformation  of 
measurements  made  by  the  user  into  understandable  measures  of  the  precision 
and  accuracy  of  the  position  determined.  An  approach  which  makes  use  of 
bivariate  statistics  is  presented  for  positioning  calculations. 


4 


2.4  Standards 


Section  4.0  conveys  a  variety  of  position  error  measures  that  could 
be  used  for  quantitative  standards  for  the  accuracy  and  precision  with  which 
the  positions  of  aids  to  navigation  will  be  established.  Standards  must  give 
clear  guidance  to  the  positioning  team,  and  must  also  reflect  the  needs  of  the 
users  of  the  waterways.  The  standards  for  positions  of  fixed  and  floating 
aids  will  be  different.  The  results  of  the  positioning  effort  come  in  three 
forms:  (1)  measures  of  accuracy,  (2)  measures  of  precision,  and  (3)  measures 
which  indicate  accuracy  and/or  precision.  The  purpose  of  this  report  is  not 
to  evaluate  the  various  measures  but  to  define  them,  explain  the  computations 
needed  to  calculate  them  and  discuss  any  obvious  drawbacks  in  procedures  based 
on  their  use. 
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3.0  GEODETIC  SURVEY  (references  21  and  22) 


3.1  General 


Geodesy  may  be  defined  as  that  science  concerned  with  the  exact 
positioning  of  points  on  the  surface  of  the  earth  and  determination  of  the 
exact  size  and  shape  of  the  earth.  The  spherical  concept  of  the  shape  of  the 
earth  offers  a  simple  surface  which  is  mathematically  easy  to  deal  with.  Many 
navigational  computations  use  it  as  a  surface  representing  the  earth.  While 
the  sphere  is  a  close  approximation  to  the  true  figure  of  the  earth  and 
satisfactory  for  many  purposes,  to  the  geodesists  interested  in  the  measure¬ 
ment  of  long  distances,  a  more  exact  figure  is  necessary.  Since  the  earth  is 
in  fact  flattened  slightly  at  the  poles  and  bulges  somewhat  at  the  equator, 
the  geometrical  figure  used  in  geodesy  to  most  nearly  approximate  the  shape  of 
the  earth  is  an  ellipsoid  of  revolution.  The  ellipsoid  of  revolution  is  a 
figure  which  would  be  obtained  by  rotating  an  el  1 i pse  about  its  shorter  axis 
(figure  2-1) . 

Any  point  on  the  ellipsoid  is  designated  in  terms  of  latitude, 
longitude  and  height,  with  zero  starting  references  for  each.  The  three 
values,  in  addition  to  the  specifications  of  the  ellipsoid  itself,  are  called 
geodetic  datum.  Instead  of  describing  the  position  of  a  point,  P,  in  terms  of 
latitude,  longitude  and  height  as  is  normally  done  in  geodetic  surveying,  it 
can  be  described  by  use  of  a  Cartesian  coordinate  system  in  x,  y  and  z.  The 
latter  is  used  for  equations  in  this  report.  The  coordinate  systems  are 
converted  from  one  to  another  by  mathematical  conversion  formulae  (appendix 
A.  1) . 


The  customary  ellipsoidal  earth  model  has  its  semi-minor  axis,  b, 
parallel  to  the  rotational  axis  of  the  earth.  The  size  of  such  an  ellipsoid 
is  usually  given  by  the  length  of  the  two  semi -axes  or  by  the  semi -major  axis, 

a,  and  the  eccentricity,  e.  Any  ellipsoid  which  is  very  close  in  shape  to  a 
sphere  is  called  a  spheroid.  The  figure  of  the  earth  considered  as  a  mean 
sea-level  surface  extended  continuously  through  the  continents  is  called  the 
qeoid.  The  ellipsoid  (or  spheroid)  is  chosen  such  that  it  closely  approxi¬ 
mates  the  geoid  in  the  region  of  interest.  Many  different  ellipsoids  have 
been  defined  for  the  various  regions  of  the  earth.  The  ellipsoid  used  by  the 
Coast  Guard  and  herein  is  the  Clarke  Spheroid  of  1866  (Ref  6).  The  Clarke 
Spheroid  has  a  major  semi-axis,  a,  of  6,378,206.4  meters,  a  minor  semi -axis, 

b,  of  6,356,583.8  meters,  and  an  eccentricity,  e,  of  8.2271854  x  10~2.  One 
datum  which  uses  the  Clarke  Spheroid  as  the  mathematical  model  of  the  earth, 
is  North  American  Datum  of  1927  (NAD27)  (reference  6).  The  NAD27  also  serves 
as  the  datum  for  nautical  charts  on  which  aid  positions  are  depicted. 

The  primary  procedures  for  positioning  buoys  for  many  years  have 
been  graphical  using  the  three-arm  protractor  and  charts  (details  are  unimpor¬ 
tant  here).  It  will  suffice  to  state  that  the  procedures  are  inaccurate  due 
to  many  causes.  If  positioning  an  order  of  magnitude  better  than  that  of  the 
navigator  is  required,  then  the  Coast  Guard  must  complete  the  conversion  to 
hydrodetic  procedures  using  the  NAD27  or  an  equivalent  datum. 


a  =  ONE-HALF  OF  THE  MAJOR  AXIS  =SEMI-MAJOR  AXIS 


b  =ONE-HALF  OF  THE  MINOR  AXIS  =  SEMI-MINOR  AXIS 

2  *2_  b2 

e  =  2 

a 


FIGURE  2-1.  Elements  Of  An  Ellipse 
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3.2  Survey  Procedures 


3.2.1  General 

In  attempting  to  develop  adequate  hydrodetic  survey  proce¬ 
dures,  a  thorough  research  of  proven  geodetic  survey  procedures  is  the  logical 
first  step.  Where  possible,  direct  adoption  of  applicable  geodetic  procedures 
is  the  most  efficient  way  to  reach  high  standards  of  accuracy  and  precision. 

Unfortunately,  such  an  efficient  route  is  not  amenable,  except  for  fixed  aids 
to  navigation.  Within  the  U.S.  Coast  Guard,  resection  is  the  primary  method 
of  survey  upon  which  hydrodetic  procedures  are  based.  Resection  is  the  deter¬ 
mination  of  the  horizontal  position  of  a  survey  station  by  observed  directions 
from  the  station  to  points  of  known  position.  The  primary  source  for  geodetic 
survey  standards  and  specifications,  reference  25,  is  void  of  marine  geodetic 
standards  for  resection.  Standards  and  general  specifications  for  hydrograph¬ 
ic  survey  are  published  by  the  Department  of  Commerce  in  the  Hydrographic 
Manual  which  was  prepared  by  M.J.  Umbach  in  1976.  Unfortunately,  these  stan¬ 
dards  and  specifications  are  not  directly  usable  by  current  Coast  Guard  posi¬ 
tioning  teams  in  their  present  form.  This  is  primarily  due  to  inadequate 
surveying  expertise  and  surveying  equipment  throughout  the  Coast  Guard. 

r 


In  order  for  the  reader  to  fully  appreciate  the  complexity  of 
an  ordinary  acceptable  geodetic  survey,  the  following  section  on  selected 
geodetic  survey  procedures  is  provided.  As  hydrodetic  procedures  for  planning 
and  positioning  are  introduced,  it  is  worthy  to  note  the  frequent  parallels 
with  the  stages  of  geodetic  survey. 

3.2.2  Procedures 

Although  accuracy  and  precision  are  often  used  interchange¬ 
ably  in  everyday  conversation,  they  have  distinctly  different  meanings. 
Accuracy  relates  to  the  quality  of  a  result,  and  is  distinguished  from  pre¬ 
cision  which  relates  to  quality  of  the  operation  by  which  the  result  is 
obtained.  In  other  words,  precision  is  the  degree  of  refinement  in  the 
performance  of  an  operation  or  in  the  statement  of  the  result.  Precision  is 
of  no  significance  unless  accuracy  is  also  obtained. 

A  primary  reference  concerned  with  the  accuracy  and  precision 
of  geodetic  survey  is  the  national  "Classification,  Standards  of  Accuracy,  and 
General  Specifications  of  Geodetic  Control  Surveys,"  published  by  the  Depart- 
ment  of  Commerce.  T'ne  purpose  of  this  section  is  to  discuss  the  various 
requirements  in  that  reference. 

3. 2. 2. I  Classification 

Surveys  are  classified: 

First  Order 

Second  Order  -  Class  I 
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Second  Order  -  Class  II 

Third  Order  -  Class  I 

Third  Order  -  Class  II 

with  First  Order  providing  the  best  accuracy  and  most  precision  with  subse¬ 
quent  orders  becoming  less  accurate  and  less  precise.  The  standards  are 
provided  for  three  basic  methods  of  survey: 

a.  tri i  angulation 

b.  trilateration 

c.  traverse 

Any  one  could  be  used  in  this  section  but  analogies  exist  such  that  only 
tri angulation  will  be  discussed. 

The  surveyor  recognizes  the  existence  of  error.  The 
procedures  specified  are  directed  at  eliminating  errors  or  at  least  reducing 
the  effects  of  those  that  are  not  compensatable.  For  a  survey  to  qualify  as  a 
given  order,  procedures  and  results  must  be  both  precise  and  accurate  within 
the  specifications.  The  accuracy  of  a  survey  is  expressed  as  the  dimension¬ 
less  quantity,  closure.  Closure  is  a  ratio  that  is  calculated  by  comparison 
of  a  calculated  value  against  some  "true"  value.  It  is  not  in  dimensions  of 
distance.  It  is  based  on  neither  an  assumed  nor  an  expected  discrepancy,  but 
on  an  actual  comparison  with  "truth."  It  is  important  to  clarify  the  concept 
of  "truth"  as  it  applies  to  surveying.  If  any  order  survey  starts  from  and 
finishes  at  higher  order  stations,  the  starting  and  finishing  points  are 
considered  "truth."  "Truth"  is  often  obtained  astronomically  and  in  more  and 
more  cases,  by  satellite.  If  a  survey  starts  and  ends  from  stations  of  the 
same  order  as  the  survey,  those  stations  are  best  described  as  "tentative 
truth,"  subject  to  subsequent  adjustment.  A  survey  of  a  desired  order  cannot 
start  or  end  on  lower-order  stations. 

3. 2. 2. 2  Specifications  (Tri angulation) 

The  first  specification  is  entitled  recommended 
spacing  of  principal  stations.  Higher  order  stations  are  general ly  spaced 
further  apart  than  lower  order  stations.  The  data  in  any  local  survey  is 
carried  through  the  principal  stations  to  the  global  survey  network  of  which 
the  principal  station  must  be  a  part.  In  order  to  conduct  an  accurate  survey, 
it  is  important  to  tie  into  the  global  survey  network  at  least  as  often  as  the 
recommended  spacing. 


The  second  specification  is  entitled  strength  of 
figure.  To  avoid  a  complete  discourse  on  this  specification,  it  is  basically 
an  exponential  figure-of-meri t  representing  the  precision  of  computed  lengths 
in  a  tri angulation  net  as  determined  by  the  size  of  the  angles,  the  number  of 
conditions  to  be  satisifed  and  the  distribution  of  base! ines  and  points  of 
fixed  position.  The  specifications  place  desirable  limits  of  the  strength-of- 
figure  factors  for  each  triangle  in  the  survey  chain.  Since  the  measure  is 
exponential  in  nature,  the  propagation  of  error  through  a  chain  of  triangles 
is  additive.  Desirable  and  maximum  limits  are  placed  on  strength-of-f igure 
before  it  is  necessary  to  tie  into  an  existing  base! i ne  or  to  measure  one. 
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The  two  specifications  above  are  applicable  to  the 
planning  of  the  survey.  They  are  based  on  a  measure  of  expected  accuracy  and 
precision,  including  propagation  of  errors.  The  upper  limits  imposed  on 
individual  triangles  and  on  the  chain  are  determined  by  the  ultimate  accuracy 
which  must  be  achieved. 


The  third  specification  is  the  first  applicable  to 
actual  conduct  of  a  survey.  This  is  entitled  base  measurement  and  requires 
that  the  starting  basel ine  (and  any  other  baseline  which  is  calculated  from 
the  points  being  surveyed)  be  measured  in  length  by  acceptable  procedures  so 
that  the  standard  error  of  the  multiple  measurements  meet  specified  stan¬ 
dards.  8ase  measurement  also  provides  a  validity  check  for  the  baseline. 

The  fourth  specification  is  horizontal  directions 
(or  angles)  which  specify  the  type  of  instrument,  the  number  of  observations, 
and  the  rejection  limits  for  those  observations.  For  example,  First  Order 
requires  a  0.2"  instrument  and  16  positions  with  a  rejection  limit  of  4  seconds 
from  the  mean.  A  position  constitutes  a  direct  and  a  reverse  measurement 
between  the  two  sighted  points.  This  entails  four  pointings  and  readings  per 
position  (total  of  64  for  First  Order).  Each  position  entails  the  determina¬ 
tion  of  the  mean  angle  taken  from  the  direct  and  reverse  readings  (for  First 
Order,  a  total  of  32  angles). 

The  fifth  specification,  triangle  closure,  consists 
of  two  parts:  They  are  "average,  not  to  exceed"  and  "maximum,  seldom  to 
exceed."  The  three  angles  measured  in  each  triangle  are  added  together.  The 
difference  between  their  total  and  180  degrees  plus  spherical  excess,  is  the 
triangle  error  of  closure.  Since  the  triangle  cannot  be  uniquely  defined  if 
misclosure  exists,  it  is  removed  by  applying  it  in  equal  portions  to  all 
angles  (residual  error  averaging).  If  misclosure  is  excessive,  the  measure¬ 
ments  are  repeated. 


The  sixth  specification  is  entitled  side  checks  and 
concerns  sides  that  are  common  to  two  different  tri angulation  chains.  Side 
equations  are  used  to  ensure  that  the  length  and  direction  are  the  same  using 
either  chain. 


The  seventh  specification  is  entitled  astro  azimuths 
and  specifies  the  frequency  (figure  spacing)  at  which  astronomical  sightings 
must  be  made.  Requirements  are  imposed  on  astronomical  sightings  that  are 
similar  to  those  imposed  on  horizontal  sightings. 

The  eighth  specification  is  vertical  angle  observa¬ 
tions  and  sets  forth  the  number  of  vertical  observations  the  allowable  spread 
and  the  number  of  figures  allowed  between  known  elevations.  This  specifica¬ 
tion  is  needed  to  reduce  horizontal  control  dafa  to  the  reference  ellipsoid. 
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3. 2. 2. 3  Accuracy 

The  final  requirement  is  that  of  closure  in  length 
and  position.  Having  satisfactorily  met  the  first  eight  specifications,  this 
specification  concerns  the  comparison  of  computed  lengths  of  bases  with  their 
measured  lengths  or  with  known  lengths  from  finishing  on  existing  control  or 
the  comparison  of  computed  positions  with  their  known  positions.  This  mis- 
closure,  if  within  tolerances  specified,  indicates  the  achievement  of  desired 
accuracy  at  the  specified  endpoints;  having  also  met  the  stringent  procedural 
and  variability  contraints  from  start  to  finish. 

3. 2. 2. 4  Summary 


a.  Instruments  used  must  meet  rigid  specifications 
on  inherent  accuracy,  precision  and  resolution. 

b.  Procedures  specified  for  use  of  the  instruments 
involve  the  averaging  of  multiple  readings  to  obtain  a  final  value,  provide  a 
means  of  averaging  out  systematic  errors  and  force  the  standard  errors  in  the 
measurement  to  be  small  and  within  calculable  limits  through  the  inclusion  of 
rejection  limits. 

c.  Residual  error  averaging  is  used,  where  resi - 
duals  are  acceptably  small,  to  ensure  consistency  of  figures. 

d.  Propagation  of  errors  is  controlled  by  frequent 
validation  against  "truth"  (azimuth  and  baseline). 

e.  Accuracy  is  not  assumed  until  misclosure  to 
truth  is  adequately  small,  precise  procedures  have  been  followed,  and  vari¬ 
ability  limitations  at  all  levels  have  been  met. 

f.  The  geometry  of  a  survey  is  planned  such  that 
positive  results  can  be  expected. 

g.  Lease  squares  adjustment  is  performed  after 
completion  of  steps  a.  through  f.  above. 

3.3  Geodesics 


A  geodesic  is  the  line  of  shortest  distance  between  any  two  points 
on  any  mathematically  defined  surface.  Both  planning  and  positioning  require 
accurate  calculation  of  the  length  and  azimuth  of  geodesics  between  points  on 
the  reference  spheroid.  The  reason  for  this  requirement  will  become  evident 
in  later  sections  where  lines  of  position  are  determined.  Procedures  for 
length  and  azimuth  computation  are  readily  available  in  the  literature  (refer¬ 
ences  8,  23,  and  24  to  name  a  few).  The  computations  can  be  simplified  by 
using  various  approximations  in  the  calculations. 

Three  procedures  for  approximating  geodesic  length,  one  procedure 
for  approximating  geodetic  azimuth,  and  the  accuracy  of  each  approximation  are 
presented  herein.  The  first  approximation  to  the  length  of  the  geodesic  is 
made  by  approximating  the  shape  of  the  reference  spheroid  by  a  sphere;  the 
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radius  of  which  is  calculated  by  the  arithmetic  average  of  the  radii  of 
curvature  of  the  geodesic  at  the  endpoints.  The  length  of  the  chord  between 
the  two  endpoints  of  the  geodesic  is  calculated  and  used  to  find  the  length  of 
the  great  circle  arc  that  connects  the  endpoints  and  lies  on  the  approximating 
sphere. 

The  second  approximation  to  the  length  of  the  geodesic  is  made  using 
only  the  endpoint  from  which  the  survey  originates  to  define  the  radius  of 
curvature  of  the  geodesic.  The  remaining  calculations  are  identical  to  those 
in  the  first  approximation  described  above. 

The  third  and  final  approximation  to  the  length  of  the  geodesic  is 
made  by  calculation  of  the  chord  length  between  the  two  endpoints  of  the 
geodesic. 

The  accuracy  of  the  approximations  decrease  from  the  first  to  the 
third  approximation  and  is  dependent  on  geodesic  length,  azimuth  and  latitude 
in  each  case. 

The  azimuth  of  the  geodesic  at  the  point  of  interest  is  approximated 
by  calculating  the  angle  between  the  chord  and  the  meri di an  (line  of  constant 
longitude)  as  projected  onto  a  plane  tangent  (touching  only  at  the  point  of 
interest)  to  the  reference  spheroid.  The  forward  azimuth  is  of  the  geodesic 
at  the  survey's  origin  and  the  back  azimuth  is  of  the  geodesic  at  the  opposite 
end  of  the  geodesic  (see  figure  3-2).  For- short  geodesics,  the  forward  and 
back  azimuths  do  not  differ  significantly  from  180°.  However,  the  differ¬ 
ence  between  forward  and  back  azimuth  for  long  geodesics  can  differ  signifi¬ 
cantly  from  180°  due  to  convergence  of  the  meridians.  The  amount  it  differs 
from  180°  is  called  convergence  (also  conversion  angle  by  navigators). 

The  following  sections  provide  additional  detail  on  each  approxima¬ 
tion  and  provide  a  description  of  the  accuracy  of  each.  The  mathematical 
details  are  contained  in  appendix  A-l. 

3.3.1  Approximating  Length  of  Geodesics 
3. 3. 1.1  Approximation  GL2 

The  procedure  for  approximating  geodesic  length 
using  both  endpoints  of  the  geodesic  is  as  follows: 

a.  define  reference  spheroid 

b.  define  geodetic  coordinates  of  endpoints  of 
geodesic 

c.  transform  geodetic  coordinates  of  endpoints  to 
Cartesian  coordinates 

d.  calculate  length  of  chord  connecting  endpoints 
found  in  (c.)  above  by  using  Pythagorean  relation 


e.  calculate  radius  of  curvature  of  geodesic  at 
both  endpoints  and  find  arithmetic  average 

f.  define  approximating  sphere  with  radius  equal  to 
radius  of  curvature  of  geodesic  found  in  (e.) 

g.  calculate  length  of  arc  on  sphere  defined  in 
(f.)  and  subtended  by  the  chord  defined  in  (d.) 

The  mathematical  details  of  this  set  of  calculations 
are  found  in  appendix  A-l.  Example  results  of  geodesic  length  approximations 
are  presented  in  table  3-1  at  two  latitudes  (0°N,  45°N)  and  fifteen 
geodesic  length  and  azimuth  combinations. 

From  the  table,  it  is  evident  that  approximation  GL2 
exceeds  all  accuracy  requirements  on  geodesic  length  for  aid  positioning 
operations.  For  distances  greater  than  1000  kilometers,  an  accuracy  of  within 
ten  meters  can  be  expected.  Within  normal  positioning  scenarios  (geodesics 
less  than  20  km)  the  accuracy  is  within  millimeters. 

3. 3. 1.2  Approximation  GL 1 

The  steps  in  approximating  geodesic  length  using  one 

endpoint  are  as  follows: 

a.  define  reference  spheroid 

b.  define  geodetic  coordinates  of  endpoints  of 
geodesic 

c.  transform  geodetic  coordinates  of  endpoints  to 
Cartesian  coordinates 

d.  calculate  length  of  chord  between  points 

e.  calculate  radius  of  curvature  of  geodesic  at 
endpoint  of  interest 

f.  define  approximating  sphere  with  a  radius  equal 
to  the  radius  of  curvature  of  geodesic  found  in 
(e.) 

g.  calculate  length  of  arc  on  sphere  defined  in 
(f.)  and  subtended  by  the  chord  defined  in  (d.) 

The  mathematical  details  of  this  set  of  calculations  are  found  in  appendix 
A-l.  Example  results  of  geodesic  length  approximations  are  presented  in  table 
3-1  at  two  latitudes  (0°N,  45°N)  and  fifteen  geodesic  length  and  azimuth 
combinations. 


From  the  table,  it  is  evident  that  approximation  GL 1 
exceeds  all  accuracy  requirements  on  approximating  geodesic  lengths  in  aid 
positioning.  For  distances  greater  than  1000  kilometers  an  accuracy  of  within 
40  meters  can  be  expected.  Within  normal  positoning  scenarios  the  accuracy  is 
within  mil  1 imeters. 

3. 3. 1.3  Approximation  GLC 


The  procedure  for  approximating  geodesic  length  by 
chord  length  is  as  follows: 


14 


ERRORS  IN  APPROXIMATIONS  TO  GEODESIC 
LENGTH  AND  AZIMUTH 


o  o  o  o  o 


o  o  o  o  o 


o  o  o  o  o 


S*  ~Q  ‘o  ”o  ~0 

o  o  o  o 

So  "o  o  "o 
o  o  o  o 
o  o  o  a  o 
o  o  o  o  o 

O'  O'  O'  O'  O' 


8  8  8S8 


ooooo 


o  o  o  o  co 


ooooo 


o  o  o  o  o 


88888 

o  o  o  o  o 

o  o  o  o  o 

cc  CO  ao  CO  00 


3  0  0  0  CM 


c’l  f'l  O'  3 


CO  CO  CO  CD  M 

«r  <T  XT  C  N 

o  o  o  o  o 

C  C  XT  XT  <T 

m  to  m  rn  co 


V  V  CM  O' 
•  3  -  O 
cm  -  cm  O 

v  v  p** 


r  s  •  ift  tr 

cm  *r  v  cm  o 

.  .  O  -  o 

xt  N*  cm 
O  XT  »>*  XT  O 


o  o  o  o 


CM  *■—  CO  •— * 


v  o  m 

CSJ  (\J  CSJ  XT  xff 


0  0  0  3  3 


CM  CM  C\j  CM 
CM 


O  O  O  O 


m  m  m 


CM  CM  CM  CM  CM 


3  3  O  3  O' 


O  O  O  O  O 


O  O  O  O  O 


8  8  888 

3  8  8  8  *8 

o  o  o  o  o 

O  3  O  O  O 

CO  CO  cc  cc  co 


CC  CO 
1-0  ro 
O  O 


to  O 
CM  in 


wrt 

Irt 

£ 

S  3  8  8  3 

8  8  8  8  8 

CD  CO  D  CO  xT 

XT  xy  tt  xy  — < 

O  O  ro  —  CM 

O  O  O  CM  ro 

8  8  8  8  8 

co  o» 
ro  ro 

-3 

3 

3  0  3  0  0 

rv>  fs.  r*«» 

8  3  8  8  8 

OOOOO 

OOOOO 
xr  v  xe-  *y  in 

OOOO 

O  O  O  O  ro 
r^,  r»*,  rv*. 

8  8  8  8  8 
c  0  0  0  C 

O  O 

«T  xy 

CM  CM 

O  o  o  o  o 


vO  cm 

rs.  xr 

3  O'  .  vO 
xr  -T  •  CM  CM 

in  •  m  cm  ro 

•  in  in  r-> 

iTI  lO  Li  — «  •—> 

co  co  cc  —  — 


•  rv  cnj  in  xr 

-t  «r  xr  o  o 

CO  CO  co  —  — 


=  s  s 

O'/)  ' 

CM  CM  •  O' 
—  —  •  O'  O' 

vn  .  ^  o'  — « 

•  in  vr»  co  — 

-x  -  -  >o  lO 

vo  vo  o  in  vn 

CM  CM  CM  •— x  «— x 


XT  *T  «T  3 


in  vn 

CM  n*. 

vO  CO  •  XT 
CM  CM  *  ro  — 

CM  •  CM  *T  O' 

•  CM  CM  —  o 

3  in  in  —  — 

co  co  co  —  —x 


•  3 

•v.  — ^ 

CM  CM 
CM  CM 


w  O  O  O 

Q) 

E 


0  0  0  —  0 


O  O' 

^  XT  C  .  rn 

u  vn  •  m  cm 

a>  •  vn  vn  rn 

333  — 

CD  co  co  co  — 

E  -t  —  —  — 


r»*  r-«. 
m  vn 
cm  cm  •  co 


co 

•  cm  cm  vn  in 

^  ^  O  O 

CO  CO  CO  -  - 


vn  xt  —i  ■ — > 


•  CM 

ro  CM 
O  — 
•  in  m  o'  m 

—  —  —  VO  vo 

vo  vo  vo  vn  in 

CM  CM  CM  —  — • 


O 

vn  o  O' 
r**.  r-*  O  • 
cm  —  •  •  m 

rr  •  —  O  O 

•  ^  ^  in  o 

—  —  —  CO  CO 

ro  ro  ro  CO  CC 

—  —  «-h  r**  p*». 


in 
vn  r^ 
r->  — 

vO  00  *00 

cm  cm  •  3  — 

CM  •  CM  ^  ro 

•  CM  cm  —  CM 

in  in  in  —  — 

CO  ao  35  -  - 


— «  O' 

— •  o 


O  O  O  — 


■I  + 1 


O  O  O  —  — 


CM 

r-*. 

CO  CO 

O  O'  .00 
<S!  ^  V  .  ro  ro 

i-  vn  »  ro  cm  cm 

a»  »  in  in  ro  ro 

w  3  3  3  —  — 

^  oo  ao  ao  —  — 


O'  — 
vn  m 

CvJ  CM  •  rs. 


«  cm  cm  vn  vn 

^  ^  V  o  o 

CO  CO  CO  —  — 


CM  CM  •  CO 

—  —  •  m  V 

in  .^*o  — 

»  in  m  O'  vn 

-  -  -x  <D 

vo  vo  vo  in  m 

cm  cm  cm  —  — 


CM  —  •  •  no 

V  •  —  o  o 

•  ^  ^  vn  o 

—  —  —  CO  CO 

ro  -o  ro  cq  rp 


VO  CO 
ro 

vo  co  •  in 
CM  CM  •  in  CM 

CM  •  CSJ  XT  n 

•  CM  CM  — •  CM 

m  o'  in  «  ^ 

CO  CO  CC  —  — X 


—  O' 

•  o 


< 


* 

s  s  r  :  O 

8  8  8  8- 

8  8  8  -S  -S 

c 

3  :  :  x  0 

8  8  8  8- 

8 

8  3  2  8  1 

Q-088 

OOOOO  00  —  00 

8  *3  2  8  § 

OOOOO  S 

8  8  8  2.  2 

8  8  8  2  2 

OOOOO 
—  O  O  —  — 

0 

0 

o  o  o  o  o 


S3: 


g  3  8  8  -8 

•o  3  2  8  ? 

8882,2 


88  88. 1 

8 

& 


:  (  x  i  © 

<3  0  00  0 


—  008 

0  —  00 

OOo 

1,00 

1,01 

t.io 

1,00 

0 

0 

0 

0 

0 

s'= 

0000 

0  0  —  — 

8  8  8  2.  2 

8  S 

% 


o  o  o  o 


o  o  o  o  o 


o  o  o  o  o 


m  u->  in  in  m 


in  m  vn  vn  m 


in  vn  vn  in  vn 


15 


6‘Z$  IboZ  „Z‘lS,lbo/0  ,,0-6b,Motbl  „/  'OP  ,9So0fC  C0I  *E  Ob’ZEIBin  6  EC’IBbOZCl  0I<  Sb'68bOZEI  „00,00o0I  ..OO.OOoOl 
„8'/P,Zfr  „8‘/*\Zb  .Z'EC.lEoPH  .S' IZ .bloSZC  P  ZS'l/RSfl  Z-0l*Z  EI  P/8SEI  1+  Wb/BSEl  ..OO.OOol  .OO.OOol 

„6‘bO,Z  .6bO,Z  „Z'6Z,ZEobH  „9'K ,  bbot-ZC  2-01*e  OZ'OO III  D  IZ’OOZZZ  0  IZ'OO/ZZ  ..OO.OIoO  „00,0Io0 


r 


a.  define  reference  spheroid 

b.  define  geodetic  coordinates  of  points  at  end  of 
geodesic  of  interest 

c.  transform  geodesic  coordinates  of  points  to 
Cartesian  coordinates 

d.  calculate  length  of  chord  between  points 

The  mathematical  details  of  this  set  of  calculations  are  found  in  appendix 
A-l.  Example  results  of  geodesic  length  approximations  are  presented  in  table 
3-1  at  two  latitudes  (0°N,  45°N)  and  fifteen  geodesic  length  and  azimuth 
combinations . 


From  the  table,  it  is  evident  that  approximation  GLC 
exceeds  all  accuracy  requirements  on  geodesic  lengths  for  normal  positioning 
problems. 


The  accuracy  of  approximation  GLC  is  within  centi¬ 
meters  for  geodesics  as  long  as  20  kilometers;  those  that  are  normally  encoun¬ 
tered  in  positioning  aids  to  navigation.  However,  its  use  for  geodesics  of 
greater  than  100  kilometers  can  cause  errors  as  large  as  five  meters.  GLC 
accuracy  is  unacceptable  for  long  geodesics  (  1000  km)  where  the  approximation 
is  only  within  thousands  of  meters  of  the  correct  length. 


3.3.2  Approximating  Azimuth  of  Geodesics  (GAC) 


follows: 


The  procedure  for  approximating  geodesic  azimuth  is  as 


a.  define  reference  spheroid 

b.  define  geodetic  coordinates  of  points  at  ends  of 
geodesic  of  interest 

c.  transform  geodetic  cordinates  of  points  to  Cartesian 
coordinates 

d.  calculate  length  and  three-dimensional  direction  of 
chord 

e.  calculate  length  and  direction  of  projection  of 
chord  onto  plane  tangent  to  spheroid  at  point  where 
azimuth  of  geodesic  is  to  be  calculated 

f.  define  direction  of  meridian  on  tangent  plane  to  the 
north 

g.  calculate  difference  between  direction  of  meridian 
and  direction  of  the  projected  chord  found  in  (e.) 

h.  convert  difference  angle  to  azimuth 

The  azimuth  at  the  opposite  end  of  the  geodesic  can  be  approximated  in  an 
identical  manner  but  with  the  endpoint  coordinates  interchanged. 

The  mathematical  details  of  azimuth  approximation  are 
contained  in  appendix  A-l.  Examples  of  the  accuracy  of  the  azimuth  approxima¬ 
tions  are  displayed  in  table  3-1  under  various  conditions.  To  check  the 
accuracy  of  azimuth  calculations  under  various  conditions,  the  forward  and 
back  azimuths  were  approximated  by  approximation  GAC  and  their  difference  was 
compared  to  the  exact  conversion  angles  using  a  precise  formula.  The  conver¬ 
sion  angles  agree  to  within  one  half  of  a  minute  for  geodesics  less  than  1500 
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kilometers  in  length  and  are  exact  to  a  tenth  of  a  second  for  lengths  up  to 
150  kilometers.  This  accuracy  exceeds  any  requirements  on  normal  aid  posi¬ 
tioning. 


3.3.3  Accuracy  of  Approximations 

Any  of  the  three  geodesic  length  approximations  are  more 
than  adequate  for  computations  normally  used  in  positioning  aids  to  naviga¬ 
tion,  where  chord  lengths  infrequently  exceed  20  kilometers. 
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4.0  PLANNING 


Planning  must  be  performed  prior  to  the  execution  of  positioning  effort 
to  establish  authoritative  standards  and  to  ensure  the  best  possible  combina¬ 
tion  of  measurements  is  selected.  The  planner  should  study  all  available 
measurement  combinations  mathematically  and  specify  reasonably  achievable 
standards  applicable  to  each  aid  location.  For  example,  it  would  be  of  no  use 
to  plan  a  positioning  effort  by  resection  methods  if  the  nearest  landmarks  are 
100  kilometers  away.  Likewise,  it  would  not  be  wise  to  set  low  standards  for 
aids  that  mark  dangerous  areas. 

The  remainder  of  this  section  is  divided  into  two  vastly  different 
parts.  The  first  part  relates  to  general  requirements  prior  to  positioning 
and  the  second  concerns  calculation  of  lines  of  position. 

4.1  Requirements  Prior  to  Aid  Positioning 

The  objectives  prior  to  aid  positioning  are  to:  (1)  identify  the 
measurement  combinations  which  are  predicted  to  provide  the  best  aid  position¬ 
ing  data  (e.g.,  accuracy  and  precision),  (2)  to  project  expected  results  of  a 
positioning  effort  using  the  selected  combinations,  and  (3)  to  provide  alter¬ 
native  procedures  in  the  event  problems  are  encountered. 

In  order  to  accomplish  the  first  objective,  the  positioning  process 
is  mathematically  modeled  using  assumed  measurement  variances  (appendix  B) . 

The  result  is  a  pi anninq  model .  The  planner  then  decides  which  position  error 
measures  are  most  critical  to  the  aid  location.  These  measures  may  be  of 
accuracy  and/or  precision.  If  the  aid  is  to  mark  a  very  distinct  obstacle  in 
a  channel,  then  it  would  be  wise  to  make  every  effort  to  place  the  aid  with  a 
high  probability  of  marking  that  obstacle,  i.e.,  accurately  and  precisely.  If 
the  aid  only  marks  the  side  of  a  channel,  then  reasonable  error  in  the  direc¬ 
tion  of  traffic  flow  may  be  acceptable. 

Let  the  position  error  measure  be  called  S.  Assume  that  there  are 
many  available  measurement  combinations  at  a  particular  station.  By  assuming 
measurement  variances,  expected  values  of  S  can  be  calculated  using  the 
planning  model.  The  measurement  combination  that  provides  the  best  position 
error  measures,  S,  is  the  first  priority  for  use  in  positioning.  A  priority 
listing  of  the  combinations  can  be  created  for  each  station. 

Objective  (2)  above  can  be  achieved  by  a  thorough  evaluation  of  the 
circumstances  due  to  the  location  of  the  aid.  The  planner  must  consider 
signal  availability  and  quality,  the  importance  or  criticality  of  the  aid, 
water  depth,  watch  circle  and  past  history  of  the  aid.  With  all  of  this  in 
mind,  the  planner  can  establish  practical  standards  for  the  aid.  Objective 
(3)  can  be  met  partially  by  use  of  the  planning  model.  Assuming  that  the 
expected  results  are  not  achieved  on  the  first  try,  alternative  procedures 
should  suggest  which  may  include  (but  certainly  are  not  limited  to):  checking 
instruments  for  error,  making  redundant  measurements,  checking  for  blunders, 
redefining  geometry  to  the  next  on  the  priority  list,  checking  landmarks, 
adding  measurements  or  simply  repeating  the  entire  positioning  evolution.  The 
planner  must  be  made  aware  of  the  problems  incurred  so  that  future  planning 
can  solve  the  problems  or  avoid  similar  occurrences  in  later  positioning 
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efforts.  The  planning  model  should  be  employed  to  determine  which  of  the 
above  procedures  makes  the  most  significant  improvement  in  position  error. 

For  example:  If  a  measurement  is  to  be  repeated,  which  one  should  be  repeated? 

4.2  Lines  of  Position 


Appendices  B  and  C  provide  a  mathematical  discussion  of  planning  and 
positioning  models.  The  starting  point  of  each  appendix  is  a  set  of  lines  of 
position  (LOP)  and  the  gradient  vector  of  each.  The  gradient  vector  is  de¬ 
fined  by  the  transverse  displacement  of  the  respective  line  of  position  per 
unit  positive  change  in  the  measurement  used  to  define  the  line  of  position 
and  is  in  the  direction  of  the  transverse  displacement.  The  differences  be¬ 
tween  observed  measurements  and  the  measurements  expected  at  the  designated 
position  are  also  needed  for  positioning. 

The  procedures  discussed  in  section  3.3  are  used  to  compute  in¬ 
verses  ,  from  which  values  of  measurements  that  determine  lines  of  position 
which  pass  through  the  designated  position  are  obtained.  To  accomplish  this 
aim,  the  following  data  are  required: 

a.  geographic  coordinates  of  designated  position 

b.  geographic  coordinates  of  observable  landmarks  or  other 
horizontal  control 


Performing  the  computations  to  find  expected  values,  prior  to  making 
actual  observations,  is  called  precomputation. 


The  discussion  of  measurement  precomputation  is  grouped  within  four 
categories.  They  are: 


Category 

time  difference  (range) 
bearing 

time  difference 
horizontal  (difference) 


Type  LOP 

small  circle 
great  ellipse 
hyperbolic 

angle  small  circle 


Example  Intrument 
radar 

gyrocompass 

LORAN 

sextant 


The  mathematical  model  is  linearized  so  each  conic  section  (circle,  ellipse, 
hyperbola)  is  assumed  straight  in  the  area  of  interest. 


4.2.1  Range  by  Time  Difference  Measurement 

The  expected  range  is  the  length  of  the  geodesic  between  the 
desired  geographic  location,  Pj,  and  the  geographic  location  of  the  land¬ 
mark,  The  angle, 9,  with  which  the  line  of  position  crosses  the 
meridian,  passing  through  the  desired  position,  is  found  from  the  azimuth  /S, 
of  the  geodesic  at  Pj:  (figure  3-3) 

(  j8  +  90  0  <  0  <  90 

9  ■  {  p  -  90  90  <  /?  5270  (4-1) 

1/3+  90  270  <  p  <  360 
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Landmark  pgd  Landmark 


LORAN  Horizontal  Difference  Angie 

Figure  3-3  POSITIVE  GRADIENT  DIRECTIONS  (PGD) 


The  gradient  magnitude  of  any  range  line  of  position  is  c/2,  where  c  is  the 
velocity  of  light  in  meters  per  second  (EM  waves).  That  is,  the  LOP  is  dis¬ 
placed  c/2  meters  transversely  for  a  one  second  difference  in  total  time  of 
travel  to  and  from  the  landmark.  Of  course,  range  measurement  is  not  in 
seconds.  Seconds  are  converted  to  meters  for  instrument  readout.  The  gradi¬ 
ent  magnitude,  in  terms  of  the  readout,  is  always  one  (1)  meter  per  meter 
difference  provided  there  are  no  systematic  errors.  The  direction, Y,  of  the 
gradient  vector  is  determined  from  the  azimuth  of  the  geodesic  at  ,  by, 

'Y  =  /S  +  180°  0°  £  <  360°  (4-2) 

4.2.2  Lines  of  Position  by  Bearing  Measurement 

The  expected  bearing  measurement  is  the  azimuth  of  the 
geodesic  at  the  desired  geographic  location.  The  angle,  &,  with  which  the 
line  of  position  crosses  the  meridian  passing  through  the  desired  location  is 
identical  to  the  forward  azimuth,/?.  Because  bearing  is  normally  taken  in 
units  of  degrees,  the  gradient  magnitude  is 


m 


ir  S 


deg 
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where  S  is  the  length  of  the  geodesic  in  meters.  The  positive  gradient 
direction  is  determined  using  the  azimuth  of  the  geodesic  at  Pj  by:  (fi 
3-3)  / 


gure 


nr  - 


/3  +  270 
/S  -  90 


qo<  /5  <-9qo 
90°£/S  <  360° 


(4-3) 
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4.2.3  Hyperbolic  Lines  of  Position  by  Time  Difference  Measurement 

Time  difference  (TD)  measurements  can  provide  the  hyperbolic 
loci  of  constant  measurement  points,  as  typified  by  LORAN  (references  10, 

11).  The  use  of  LORAN  in  positioning  aids  is  generally  confined  to  a  backup 
for  periods  of  adverse  weather  conditions.  After  calibration  and  control 
procedures  have  been  performed,  buoy  positioning  in  the  repeatable,  augmenteo 
or  hybrid  modes  can  be  acceptable.  (Presently,  lORAN  use  for  aid  positioning 
is  never  acceptable  without  calibration  and  rarely  acceptable  even  with 
calibration).  Each  calibration  procedure  (mode)  is  described  adequately  in 
references  4  and  9  and  will  only  be  described  briefly  here.  It  is  sufficient 
to  discuss  LORAN  as  it  relates  to  precomputations  compatible  with  this  report. 

The  assumption  is  made  for  precomputation  of  TD's  that 
signals  propagate  along  exactly  defined  paths  (geodesics  from  transmitter  to 
receiver)  at  a  constant  velocity.  This  constant  velocity  is  best  estimated  by 
the  signal  velocity  over  a  salt  water  path.  Of  course,  errors  in  precomputa¬ 
tion  of  TO's  result  as  this  assumption  does  not  truly  reflect  actual  condi¬ 
tions  over  land  paths  or  when  local  anomoiie*  in  propagation  paths  exist.  If 
systematic  errors  are  reasonably  constant  over  an  area,  the  precomputation 
procedure  provides,  quite  accurately,  offsets  in  time  differences  between 
points  of  known  geographic  location.  The  precomputation  is  based  on  ground 
wave  and  in  no  way  considers  the  random  signal  fluctuations  in  the  propagation. 
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Precomputation  of  time  differences  at  a  desired  buoy  location 
vary  with  the  different  modes  of  LORAN  operation  (reference  9  ch.  7).  In  the 
repeatable  mode,  the  precomputed  TDs  are  not  used  as  they  are  inferior  to 
historical  TDs  obtained  when  the  aid  was  positioned  by  more  accurate  methods. 
The  augmented  mode  is  an  extension  of  the  repeatable  mode  in  which  differen¬ 
tial  corrections  to  the  TDs  are  observed  by  another  accurately  located  source 
in  the  region  and  communicated  to  the  positioning  team  at  the  time  of  the 
positioning  effort. 

Reference  9  calls  for  use  of  a  CG  computer  program  to  compute 
TDs  when  using  the  hybrid  mode  of  LORAN  operation.  The  TDs  are  mathematically 
determined  for  the  desired  position  and  at  some  other  nearby,  easy  to  find, 
reference  point.  The  vessel  visits  the  reference  point,  measures  TDs,  com¬ 
putes  corrections  and  applies  them  to  the  predicted  TDs  for  the  desired 
location.  Provided  that  the  reference  point  is  easily  accessible  and  no 
serious  local  anomolies  exist,  this  procedure  can  be  sufficiently  accurate  for 
backup  use.  The  augmented  and  repeatable  modes  are  preferable,  if  at  all 
possible. 

In  cases  where  the  CG  Headquarters' computer  program  is  not 
readily  available,  an  0PM  LORAN  precomputation  capability  may  prove  useful. 

Of  course,  the  0PM  LORAN  capability  must  also  be  locally  calibratable  and 
usable  in  the  quasi -di fferential  or  repeatability  modes  with  procedures 
promulgated  in  the  Aids  to  Navigation  Positioning  Manual  (CG-222-5).  Require¬ 
ments  for  such  an  0PM  precomputation  capab i 1 i ty  include  the  geographic  coordi¬ 
nates  of  the  desired  position  and  some  preselected  reference  location  in  its 
vicinity,  and  prerecorded  data  on  geographic  location  and  signals  for  stations 
in  the  applicable  chain.  An  example  of  the  LORAN  data  required  is  provided  by 
table  4-1.  The  table  includes  the  geographic  latitude  and  longitude  of  each 
station  in  the  Northeast  LORAN-C  chain  and  the  respective  programed  time 
delay  between  the  master  and  secondary  signal.  The  programmed  time  delays 
found  in  the  table  are  slightly  different  than  other  published  time  delays  but 
have  been  developed  through  research  associated  with  the  CAP  effort  of  ANPAR. 

The  most  efficient  use  of  LORAN-C  signals  for  positioning  is 
to  calibrate  locally  and  not  even  consider  the  programmed  time  delays. 
Calibrating  in  the  local  area,  in  accordance  with  the  ATON  Manual,  allows 
determination  of  the  TDs  necessary  at  the  desired  location. 

The  gradient  vector  magnitude  of  a  hyperbolic  LOP  is  provided 
by,  (reference  10) 
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where  c  is  the  velocity  of  electromagnetic  waves  over  sea  water  and  -  (3S 
is  the  difference  in  forward  azimuths  to  the  two  stations,  master  and  secon¬ 
dary.  The  LOP  direction,  9  ,  is  given  by:  (figure  3-3) 
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The  positive  gradient  direction  is  calculated  by, 

+A  -  90°  /^m^s 

- 2 - 

^  m  +A  +  90°  m  >/^s 

2 

/^vn  =/^s 

(depends  on  position  relative  to  two  stations) 


Table  4-1 

LORAN-C  Data  Northeast  Chain  (NAD  27) 


LATITUDE 

LONGITUDE 

PROGRAMMED  DELAY 

MASTER 

42°.  71401944 

76°. 82623333 

- 

W 

46°. 80773889 

67°. 92754444 

13797.20 

X 

41°. 25332778 

79°. 97791944 

26969.93 

Y 

34°. 06266944 

77°. 91311111 

42221.65 

Z 

39°. 85207778 

87°. 48653055 

57162.06 
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4.2.4  Circle  of  Position  by  Horizontal  Angle  Measurement 

The  expected  horizontal  angle  measurement  is  the  difference 
between  the  forward  azimuths  of  the  geodesics  connecting  the  desired  geograph¬ 
ic  location  and  the  landmarks  to  be  observed.  The  angle,  6,  with  which  the 
LOP  crosses  the  meridian,  passing  through  the  desired  location,  is  determined 
by:  (reference  9)  (figure  3-3) 


9  a  ^1  +  02  “  03  ±  1800  (4.6) 

0°  <  0  <  360° 

where  0i  is  the  forward  or  back  azimuth  of  the  geodesic  between  the  desired 
geographic  location  and  the  left  most  landmark,  02  is  the  same  measure  but 
for  the  other  landmark,  and  £3  is  the  forward  or  back  azimuth  of  the  geo¬ 
desic  between  the  landmarks.  It  is  noted  that  the  derivation  of  equation 
(4-6)  in  reference  9  is  an  approximation.  The  LOP  angle,  9,  is  found  on  a 
plane  which  approximates  the  spheroid  in  the  region  of  interest.  The  result 
of  this  approximation  is  an  incorrect  LOP  angle.  The  inaccuracy,  however,  has 
a  negligible  effect  on  later  derivations. 

The  gradient  vector  magnitude  is  determined  from  the  linear¬ 
ized  formula:  (reference  9) 


G  /roeters\  =  *  ab  (4-7) 

\minutey  (180  x  60‘ )p 

where  a  is  the  chord  length  approximation  of  the  geodesic  between  the  desired 
position  and  the  left -most  landmark,  b  is  the  same  approximation  to  the  other 
landmark,  and  p  is  the  same  approximation  between  the  two  landmarks.  The 
computational  errors  will  not  significantly  affect  the  results  of  further 
derivations.  The  positive  gradient  direction,  Y,  of  the  line  of  position  is 
determined  from  the  LOP  angle,  0,  and  the  azimuths  of  the  two  observed  land¬ 
marks  by 


9  +  90° 
0  +  90° 


01  *  02 
01  =02 


0°  —  #  <  360°  (4-8) 
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5.0  POSITIONING 

This  section  is  divided  into  two  parts:  (I)  systematic  error  and  (2) 
accuracy  and  precision.  The  first  part  concerns  error  sources  for  which  com¬ 
pensation  procedures  are  either  provided  here  or  readily  available  in  other 
references.  Once  again,  the  mathematics  are  left  to  the  appendices.  The 
second  portion  puts  accuracy  and  precision  into  perspective  for  use  in  this 
report. 

5.1  Systematic  Errors 

Compensation  for  systematic  errors  is  a  necessary  part  in  accurate 
positioning  and  can  take  place  before  or  after  computations  to  determine  a 
position,  depending  on  the  nature  of  the  error.  Compensation  for  systematic 
observation  errors  take  place  before  computations,  whereas  compensation  for 
eccentrics  (translational  offsets),  such  as  that  induced  when  the  measurement 
observers  are  not  located  at  the  chain  stopper,  can  be  performed  after  the 
computations  of  the  position  of  the  angle  takers.  Systematic  observation 
errors  such  as  sextant  error,  radar  error,  gyro  error,  personal  error  and 
inclined  angle  error  are  compensated  for  through  procedures  readily  available 
in  references  4  and  9. 


5.1.1  Observer  to  Chain  Stopper  Vector 

Observers  are  prohibited  from  standing  near  the  sinker  drop 
point  during  positioning  operations.  The  displacement  from  the  reference 
position  of  angle  takers,  AT,  to  the  sinker  drop  point  may  be  compensated  for 
after  determination  of  AT.  8oth  the  direction  and  distance  from  the  sinker 
release  point  to  AT  are  required.  The  mathematical  equations  for  this  trans¬ 
lation  are  found  in  section  C.3.  The  sinker  release  point  is  of  prime  impor¬ 
tance  as  it  represents  the  MPP.  The  AP-to-MPP  distance  is  determined  from  the 
AP-to-AT  distance  by  applying  this  eccentric  correction  either  with  the 
mathematical  equations  or  by  placing  a  scaled  model  of  the  positioning  unit  on 
the  gradient  diagram. 

5.1.2  Lack  of  Observer  Coincidence 

It  is  desirable  for  observers  to  stand  at  the  same  point  when 
making  measurements.  Of  course,  this  is  difficult  due  to  intervisibility 
conditions  and  often  a  systematic  error  from  observer  lack-of-coincidence 
exists.  The  distance  and  direction  of  each  observer  from  the  reference 
position  of  angle  taker,  AT,  are  used  to  compensate  for  this  sytematic  error. 
The  mathematical  equations  for  compensation  are  found  in  section  C.2. 

5.2  Accuracy  and  Precision  in  Positioning 

When  a  position  is  determined,  the  accuracy  is  the  difference 
between  the  actual  position  and  the  determined  position.  High  accuracy  is 
present  when  the  difference  is  small.  It  is  desirable  to  describe  how  much 
the  difference  is  so  that  a  measure  of  accuracy  can  be  assigned  to  the  posi¬ 
tion.  A  more  accurate  means  of  positioning  is  required  to  make  this  determin¬ 
ation;  but  such  a  means  is  seldom  available  so  accuracy  is,  in  general, 
inmeasurable  for  usage  in  aids  to  navigation  positioning.  The  need  remains  to 
describe,  as  best  as  possible,  the  accuracy  in  a  determined  position.  The 
best  sidestep  of  this  predicament  is  to  employ  statistical  techniques  to 
describe  accuracy  in  a  manner  which  best  suits  the  needs  of  this  report.  For 
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those  needs,  accuracy  is  described  harmoniously  with  precision.  The  accuracy 
statistic  derived  and  computed  is  supplemented  by  an  equally  important  preci¬ 
sion  statistic.  Mistakes  and  uncompensated  for  systematic  error  in  the 
observations  cannot  be  ignored  in  computing  accuracy  and  precision  statistics 
and  must  be  detected  as  often  as  possible. 

The  planning  methods  described  earlier  are  not  sufficient  for 
computing  the  accuracy  and  precision  statistics  of  a  determined  position. 

They  provide  expectation  probabilities  rather  than  actual  positioning  statis¬ 
tics.  The  probable  measurement  error,  while  indicative  of  a  large  sample  of 
error  statistics,  is  not  respresentati ve  of  a  small  unique  set  of  actual 
errors.  It  ignores  the  possibility  that  mistakes  might  have  been  made.  The 
methods  used  in  planning,  at  best,  provide  some  expectation  of  the  precision 
of  a  position  but  do  not  predict  the  accuracy  of  a  position. 

The  statistical  approach  to  positioning  is  the  obvious  extension  of 
the  probabilistic  methods  used  in  planning.  The  statistical  approach  is  based 
on  the  evaluation  of  each  particular  position  as  a  separate  entity.  The 
statistical  evaluation  is  sensitive  to  mistakes  and  uncompensated  systematic 
errors.  The  statistical  approach  provides  measures  of  both  accuracy  and 
precision  in  positioning.  A  statistical  model  and  procedures  that  provide  the 
accuracy  and  precision  statistics  are  provided  by  Rosenblatt  (reference  12). 
The  Rosenblatt  model  and  the  least  squares  principle  used  in  the  Error  Sensi¬ 
tivity  Model  (ESM)  (reference  2),  have  been  extended  in  appendix  C  to  meet  the 
requirement  of  this  report. 

Various  computational  procedures  are  available  but  all  can  be 
''escribed  in  general  terms  as  follows.  With  respect  to  a  specified  "assumed 
position"  (in  practice  the  desired  location),  each  observation  defines  a  line 
of  position.  If  there  were  no  measurement  errors,  each  LOP  would  pass  through 
the  assumed  position  with  a  known  direction,  i.e.,  the  LOPs  would  be  the 
precomputed  LOPs  which  define  the  point.  In  practice,  observed  measurements 
contain  some  uncompensated  error  and  the  LOPs  representing  the  observations 
will  not  all  pass  through  the  assumed  position.  In  the  standard  linearized 
approach  for  small  relative  errors,  the  actual  LOPs  (based  on  measurements) 
are  taken  to  be  parallel  to  the  precomputed  LOPs  and  the  distance  between 
them,  d,  is  computed  by  (reference  12), 

d  =  GAm  (5-1) 

where,  G,  represents  the  gradient  magnitude  of  the  LOP  and  Am  represents  the 
small  error  in  measurement. 

The  degree  to  which  the  measurements  are  consistent  with  each  other 
and  with  the  precomputed  measurements  form  a  sample  for  analysis  of  the 
precision  and  accuracy  of  the  determined  position.  With  the  relative  weight 
of  each  measurement  provided  by  the  type  of  measurement,  the  sample  can  be 
used  to  statistically  find  a  position  which  best  represents  the  actual  posi¬ 
tion  and  also  provides  a  source  for  determining  the  confidence  in  that  best 
representation.  The  position  calculated  is  called  the  reference  position  of 
the  angle  takers  (AT).  The  displacement  vector  V,  from  the  assumed  position 
to  AT,  is  derived  in  appendix  C.  The  precision  of  AT,  is  derived  statistical¬ 
ly  by  study  of  the  measurement  residuals.  An  unbiased  estimate,  sS  of  the 
reference  measurement  variance, <rft,  (references.  13  and  14)  is  made  and 
general ized-T^  statistics  (reference  15)  are  used  to  define  the  confidence 
regi on. 
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6.0  STANDARDS  -  MEASURING  SUCCESS  IN  POSITIONING 
6.1  Standards  -  General 


Quantitative  standards  should  prescribe,  in  numerical  terms,  the 
acceptable  degree  of  uncertainty  associated  with  the  positions  of  aids  to 
navigation.  The  objective  of  this  section  is  to  discuss  the  quantitative 
(numerical)  measures  of  accuracy  and/or  precision  with  which  the  positions  of 
aids  to  navigation  could  be  established  and  maintained.  Standards  must  be 
derived  from  the  quantitative  measures  so  that  they  will  give  clear  guidance 
to  operating  personnel  and  reflect  the  needs  of  the  mariner.  A  standard  is 
defined  as  an  authoritative  value,  Sm,  for  comparison  with  a  position  error 
measure,  S,  to  determine  if  S  is  acceptable  (which  indicates  that  the  posi¬ 
tioning  evolution  is  complete).  Sm  may  depend  on  the  specific  environment 
at  the  time  of  positioning,  the  availability  of  measurements,  the  importance 
or  criticality  of  the  aid,  and  the  difficulty  in  achieving  Sm  at  the  time  of 
positioning.  In  reality,  both  S  and  Sm  may  be  a  set  of  position  error 
measures  and  standards.  The  specific  combination  depends  on  the  aid.  Suc¬ 
cessful  aid  positioning  is  determined  by  both  accuracy  and  precision;  which 
means  that  the  set  of  S^  values  should  contain,  at  least,  standards  for 
these  two  qualities. 

Measures  of  accuracy  and  precision  need  not  be  restricted  to  use  in 
geometry  selection  (section  4.1.1)  and  first  time  positioning  evolutions. 
Audit  of  aid  to  navigation  positions  also  requires  some  appropriate  quantita¬ 
tive  measure. 


Practically  all  aid  positioning  units  within  the  Coast  Guard  have 
used  point  plotting  with  a  three-arm  protractor  and  chart  as  a  prime  method  to 
obtain  a  measure  of  success  with  established  confidence.  This  method  is 
gradually  being  replaced  by  more  sophisticated  hydrodetic  procedures,  some  of 
which  are  described  in  Aids  to  Navigation  Manual  -  Positioning  (COMDTINST  Ml 6500 . 1 ) 
The  equations  and  procedures  within  this  report  allow  for  creation  of  more 
definitive  standards  suggesting  use  of  grids  and/or  calculators. 

An  assortment  of  position  error  measures  are  discussed  in  tr.e 
following  paragraphs  with  detailed  mathematics  left  to  appendix  D.  Frequency 
histograms  of  the  position  error  measures  actually  found  through  research  of 
historical  data  are  also  included  in  the  appendix. 


6.2  Standards  -  Specific 


6.2.1 


A  Posteriori 


Estimate  of  Reference  Measurement  Variance 

f 


In  section  5.2,  accuracy  and  precision  were  placed  in  per¬ 
spective  and  defined  harmoniously.  A  statistical  approach  to  positioning 
alluded  to  the  reference  measurement  variance  (abbreviated  to  reference 
variance  or  a  £7!  The  reference  variance  is  an  arbitrary  constant  with 
arbitrary  dimensions.  Its  square  root  is  referred  to  as  the  reference  mea¬ 
surement  standard  deviation  or  just  reference  standard  devi at  ion,  cr0.  The 
precision  with  which  the  AT  can  be  defined  is  governed  by  the  stochastic 
nature  of  the  measurements  from  which  the  AT  was  calculated.  The  stochastic 
nature  (random  error)  of  the  instruments  causes  inconsistency  in  a  set  of 
measurements.  The  inconsistency  is  evident  in  the  residuals  of  a  set  of 
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measurements  which,  in  turn,  are  used  to  make  an  unbiased  estimate  of  the 
reference  variance.  The  estimate  is  called  the  "A  Posteriori  Estimate  of  the 
Reference  Variance"  and  is  a  possible  measure  for  which  standards  can  be  set 
for  precision  associated  with  AT.  The  reference  variance  estimate,  s  ,  can 
in  turn,  be  used  as  a  part  in  defining  a  confidence  region  for  AT  (described 
in  detail  later). 


The  reference  variance  might  be  considered  a  benchmark  to 
which  measurement  consistency  can  be  compared.  It  is  determined  through 
analysis  of  historical  data.  Its  units  are  considered  arbitrary  because  it 
can  be  chosen  to  represent  the  measurement  variance  of  any  instrument  to  be 
used  in  positioning.  Only  the  relative  measurement  variances  are  of  impor¬ 
tance  and  they  are  so  only  for  weighting  different  measurement  types.  In  the 
case  of  horizontal  angles  with  the  sextant,  its  units  are  minutes  squared. 
Appendix  D  provides  some  reference  variance  estimates  that  were  observed 
aboard  the  REDWOOD  (WLM-685)  during  pilot  analytical  positioning  efforts. 

6.2.2  A  Posteriori  Estimate  of  Reference  LOP  Variance 


A  procedure  similar  to  that  presented  in  the  last  section  can 
be  used  to  make  an  unbiased  estimate  of  reference  LOP  varian:e.  LOP  variance 
represents  the  random  error  associated  with  the  location  of  a  line  of  posi¬ 
tion.  The  LOP  variance  is  related  to  the  measurement  variance  through  the 
gradient  of  the  line  of  position, 
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where is  the  measurement  variance. 


The  estimate  of  LOP  variance  for  particular  measurement  sets 
are  computed  through  use  of  the  distance  residuals  which  are  the  measurement 
residuals  converted  to  distance  by  the  gradient.  The  reference  LOP  variance 
being  estimated  by  the  distance  residuals  is  an  average  of  the  LOP  variances 
of  the  set  of  LOPs  in  the  geometry.  Similar  to  reference  measurement  vari¬ 
ance,  analysis  of  historical  data  is  needed  to  establish  reasonable  standards 
for  the  LOP  variance  estimate.  Some  values  for  this  estimate  are  provided  in 
appendix  D  from  analysis  of  pilot  analytical  positioning  efforts  aboard  the 
REDWOOD  (WLM-685). 


6.2.3  Confidence  Ellipse  Parameters 

The  stochastic  nature  of  the  set  of  measurements  used  to 
determine  a  position  is  converted  to  a  two-dimensional  probability  distribu¬ 
tion  about  the  AT  via  a  mathematical  transformation.  The  two-dimensional 
distribution  is  called  a  bivariate  probability  distribution.  An  important 
consideration  here  is  that  the  variances  of  the  distribution  in  both  dimen¬ 
sions  are  unknown.  One  of  the  most  important  groups  of  problems  in  statistics 
relates  to  questions  concerning  the  mean  (the  MPP  here)  when  the  variance  is 
unknown.  Anderson  (reference  15)  thoroughly  exhausts  the  subject. 
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From  that  work,  the  general ized-T^  statistic  was  taken  for 
use  in  determining  confidence  regions  with  unknown  distribution  variances. 
The  regions  are  bounded  by  ellipses  and  therefore  are  named  confidence 
ellipses. 


The  statistical  analysis  needed  to  define  confidence  regions 
is  found  in  appendix  C.  In  that  appendix,  equations  for  computation  of  the 
following  confidence  ellipse  parameters  are  presented: 

a.  major  semi-axis 

b.  area 

c.  semi -diameter  in  specified  direction 

d.  circle  of  confidence 

All  four  of  the  confidence  ellipse  measures  are  functions  of 
the  estimate  of  the  reference  variance,  the  level  of  confidence  desired  in  the 
MPP,  the  number  of  measurements  used  to  determine  the  MPP  and  the  weight  and 
orientation  of  each  line  of  position.  Each  confidence  ellipse  parameter  is  a 
potential  position  error  measure  for  which  standards  can  be  established. 
Confidence  limits  of  9016  have  been  specified  for  field  use  (reference  9). 

6.2.4  AP-to-MPP  Vector 


The  AP-to-MPP  vector,  Vg,  discussed  in  section  5.1  is  a 
potential  measure  of  success  in  positioning  in  tandem  with  other  position 
error  measures.  By  itself,  the  vector,  like  any  vector,  has  two  parts:  a 
magnitude  and  direction. 

The  magnitude  is  a  measure  of  the  accuracy  in  the  position 
determined.  The  measure  of  accuracy  is  not  absolute  because  it  is  not  pos¬ 
sible  to  check  its  correctness  (in  geodetic  land  si.-vey  this  can  be  performed 
by  sidechecks,  final  closure,  higher  order  endpoints,  baseline  measurement, 
etc.).  The  measure  is  of  accuracy  only  in  the  sense  that  it  defines  how  close 
the  positioning  team  thinks  the  sinker  drop  point  is  to  the  point  they  are 
marking.  Together  with  a  measure  of  precision,  the  vector  is  the  best  measure 
of  accuracy  available  in  hydrodetic  survey.  Some  values  for  the  magnitude  of 
the  AP-to-AT  displacement  are  provided  in  appendix  D  from  analysis  of  pilot 
analytical  positioning  efforts  aboard  the  REDWOOD  (WLM-685).  (NOTE:  The 
AP-to-AT  has  not  been  corrected  for  observer  displacement  from  the  sinker  drop 
poi nt. ) 


The  direction  of  the  AP-to-MPP  vector  is  useful  in  situations 
where  accuracy  is  important  only  in  specified  directions.  The  component  of 
the  AP-to-MPP  vector  in  that  direction  is  the  desired  position  error  measure. 


6.2.5  P-in-R 


There  are  only  a  few  measures  that  combine  both  the  accuracy 
and  precision  of  the  MPP.  One  such  measure  is  the  probability  that  the 
determined  position  is  somewhere  within  a  designated  circular  region  of  radius 
R  centered  on  the  desired  position,  P-in-R.  Calculation  of  P-in-R  is  outlined 
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in  appendix  C  of  reference  2  and  requires  two-dimensional  numerical  integra¬ 
tion,  which  is  very  time  consuming.  Numerical  errors  related  to  its  calcula¬ 
tion  and  their  causes  are  presented  in  appendix  D.5.  The  time  required  to 
calculate  P-in-R  on  scene  using  programmable  calculators  and  equations  defin¬ 
ing  regions  of  applicability  are  also  provided  in  the  appendix. 

The  use  of  P-in-R  as  a  measure  of  success  is  restricted  to 
well  defined  probability  distribution  functions.  This  means  that  the  random 
error  associated  with  each  measurement  must  be  assumed  prior  to  its  computa¬ 
tion.  In  a  way,  this  is  a  step  backwards  from  the  statistical  methods  pro¬ 
moted  by  this  report.  However,  its  strength  as  a  measure  of  success  could 
overshadow  this  drawback. 

The  bivariate  probability  density  function  defined  about  the 
MPP  using  the  general i zed-T^  technique  is  not  well  defined  and  further 
analysis  would  be  necessary  to  define  a  measure  similar  to  P-in-R.  Although 
the  P-in-R  position  error  measure  may  be  excellent  in  planning  and  assessing 
the  success  that  is  thought  to  have  been  achieved,  it  has  limited  importance 
in  actual  positioning  efforts. 

Values  of  P-in-R  are  provided  in  appendix  D  from  analysis  of 
pilot  analytical  positioning  efforts  (observer  measurement  standard  deviation 
of  5  minutes)  aboard  the  REDWOOD  (WLM-685).  The  radius,  R,  was  set  at  10,  15 
and  20  meters  for  the  calculations. 

6.2.6  R-for-P 


The  radius  of  a  circular  region,  centered  on  the  desired 
position,  that  contains  at  least  some  specified  fraction  of  the  probability 
mass  is  called  R-for-P.  Accurate  calculation  of  this  measure  numerically  is 
beyond  the  scope  of  this  report.  However,  a  crude  but  very  conservative 
approximation  that  never  exceeds  10%  error  (the  10%  refers  to  error  in  the 
probability  mass  enclosed  by  the  radius,  R-for-P)  is  presented  in  appendix 
0.6.  The  R-for-P  approximation  can  be  made  regardless  of  whether  or  not  the 
probability  density  function  is  well  defined.  REDWOOD  data  was  also  used  to 
provide  some  representative  values  of  this  position  error  measure. 

6.2.7  Difference  Between  Actual  and  Computed  Measurements 


6.2. 7.1  All  Measurements 


After  numerous  operations  at  a  given  station,  high 
confidence  will  be  attained  in  the  expected  or  precomputed  measurements.  In 
the  cases  where  high  confidence  exists,  the  measurements  made  in  positioning 
can  and  should  be  compared  with  those  that  are  expected.  The  agreement  in 
this  comparison  can  be  measured  statistically  by  comparing  the  differences  to 
a  priori  assumed  measurement  variances.  The  standard  in  this  case  is  a 
dimensionless  number  taken  from  statistical  tables. 

This  procedure  depends  on  known  probability  density 
function  variances  and  analysis  of  historical  data  is  required  to  provide  a 
reasonable  assumed  measurement  variance.  This  position  error  measure  is 
called  the  sun  of  the  squared,  weighted  differences,  swd,  and  is  defined 
mathematically  in  appendix  D.7.  This  measure  of  success  combines  position 
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accuracy  and  precision.  Also  defined  in  appendix  D.7  is  the  gradient  weighted 
differences,  gwd.  This  measure  of  success  describes  how  far  the  LOPs  are  from 
the  AP.  It  has  no  advantage  over  swd  except  that  the  differences  are  distance 
measures.  Data  on  measurement  differences  and  gradient  weighted  measurement 
differences  are  presented  in  appendix  0.  The  differences  are  not  weighted  by 
assumed  variances. 


6.2. 7.2  One  Measurement 


One  of  the  most  frequently  employed  procedures  in 
present  Coast  Guard  aid  positioning  is  called  the  fixed  glass  procedure.  The 
fixed  glass  procedure  is  performed  by  setting  two  sextants  on  precomputed 
angles  and  maneuvering  the  ship  until  both  measurements  agree  exactly  with  the 
prescribed  angles.  Obviously,  this  procedure  is  limited  to  only  two  measure¬ 
ments  and  must  be  extended  to  satisfy  the  requirement  for  at  least  three 
measurements.  To  do  this,  an  observer  quickly  measures  a  third  angle,  each 
time  a  measurement  set  is  desired,  thus  allowing  a  check  on  the  measurement 
set.  The  problem  is  to  define  limits  within  which  the  measurement  set 
indicates  an  acceptable  measurement  set.  The  derivation  of  the  relationship 
between  the  third  measurement  difference  and  resulting  position  error  is  found 
in  appendix  D.7.  The  results  show  the  linear  relationship  that  exists  between 
the  measurement  difference  and  the  major  semi-axis  of  the  confidence  ellipse. 
The  slope  of  the  linear  function  is  dependent  on  the  confidence  level  desired, 
the  weight  of  the  third  measurement  relative  to  the  first  two  measurements, 
and  the  geometry  of  the  fix.  An  example  of  this  procedure  is  incluoed  with 
the  derivation.  Similar  calculations  can  be  performed  for  other  position 
error  measures.  The  procedure  is  easily  adaptable  to  the  presently  used  grid 
method  (reference  9)  with  or  without  using  the  fixed  glass  procedure. 

6.3  Summary-Combined  Standards 

Success  in  positioning  must  be  indicated  by  measures  of  both  the 
precision  and  the  accuracy  of  the  resulting  position.  The  position  error 
measures  defined  previously  are  listed  in  table  6-1.  The  table  indicates 
whether  the  measures  are  of  accuracy,  of  precision,  or  both. 

It  is  not  the  purpose  of  this  report  to  designate  the  measures  that 
are  most  appropriate  for  setting  standards;  however,  the  following  considera¬ 
tions  are  important: 

a.  The  numerical  measure  should  provide  immediate  feedback  to 
the  positioning  team  to  remove  reliance  on  graphical  measures 
of  success. 

b.  The  physical  significance  of  the  measure  must  be  easily 
understood  by  senior  members  of  the  positioning  team. 

c.  The  standards  set  at  any  specific  location  must  take  into 
account  the  peculiarities  of  that  location;  this  implies  that 
all  aids  are  not  of  equal  importance  to  the  mariner. 

d.  The  measure  should  be  easy  to  advertise  to  the  mariner  and 
easily  defended  in  court;  the  measure  should  not  be  too 
restrictive  on  the  positioning  team. 


POTENTIAL  POSITION  ERROR  MEASURES  FOR  USE  IN 
SETTING  STANDARDS 


measurement  difference. 


7.0  DETECTION  OF  MEASUREMENT  ERRORS  AND  REJECTION  OF  MEASUREMENTS 


Under  the  broad  concept  of  error  properties  of  observations,  the  conven¬ 
tional  theory  of  errors  includes  blunders  or  mistakes  in  addition  to  random 
errors  and  systematic  errors.  Blunders  may  be  caused  by  numerous  failures  to 
follow  prescribed  procedures.  From  a  statistical  point  of  view,  blunders  are 
observations  that  cannot  be  considered  as  belonging  to  the  distribution  in 
question.  They  should  not  be  used  with  the  sample.  Consequently,  measure¬ 
ments  should  be  planned  and  observational  procedures  designed  to  allow  for 
blunder  detection  and  rejection.  In  practice,  there  are  a  variety  of  ways  to 
detect  blunders: 


a.  Taking  multiple  measurements  and  checking  for  consistency 

b.  Careful  checking  of  all  reading  and  recording 

c.  Checking  and  verifying  performance  of  instruments 

d.  Vary  procedures  to  get  same  desired  result 

e.  Applying  simple  geometric  and  algebraic  checks 

f.  Note  which  measurements  deviate  from  the  norm  by  a  signifi¬ 
cant  amount 

Despite  precautions,  some  blunders  may  still  remain.  Their  detection  and 
rejection  should  be  carried  out  according  to  principles  of  statistical  test¬ 
ing.  These  principles  require  a  priori  knowledge  of  the  distributions  of  the 
random  variables  involved.  It  is  normally  assumed  that  the  observations  are 
normally  distributed  and  the  variance  of  each  observation  type  is  known.  The  a 
posteriori  reference  variance  estimate  discussed  in  section  6.2.1  provides  an 
excellent  vehicle  for  testing  a  particular  measurement  set  for  existence  of  a 
blunder.  Assuming  that  a  priori  knowledge  of  the  expected  variance  can  be 
established,  the  reference  variance  estimate  can  be  tested  via  the dis¬ 
tribution  (reference  14).  The  test  is  outlined  in  appendix  E  in  three  differ¬ 
ent  subsections: 


a.  A  Posterior  Estimate  of  Reference  Measurement  Variance 

b.  Measure  Combinations 

c.  Difference  Between  Actual  and  Precomputed  Measurements 

The  first  section  discusses  how  the  estimate  can  be  compared  to  the  a 
priori  known  measurement  variance  determined  from  historical  data.  A  statis¬ 
tically  significant  difference  leads  to  suspicion  of  a  blunder. 

The  second  section  discusses  ways  of  finding  which  measurement  of  a  set 
is  the  most  likely  to  contain  the  blunder. 

The  final  section  considers  another  approach;  that  of  statistical  analy¬ 
sis  of  the  actual  and  computed  measurement  differences.  Both  blunder  detec¬ 
tion  and  rejection  are  considered  using  this  procedure. 
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8.0  OPERATIONAL  PROCEDURES  FOR  POSITIONING 

The  hydrodetic  procedures  presented  can  be  adopted  for  Coast  Guard  use  in 
different  ways,  achieving  desirable  results.  Two  ways  are  explained  in  the 
following  sections:  graphical  procedures  and  calculator  based  procedures.  No 
claim  is  made  that  one  is  better  than  the  other  and,  in  fact,  a  combination  of 
desirable  qualities  of  each  may  be  the  best  track  to  follow. 

Each  section  has  been  divided  into  six  parts:  description,  accuracy, 
adaptability,  ease  of  use,  feedback,  and  recording.  The  description  section 
is  self-explanatory.  The  accuracy  section  is  concerned  with  the  potential 
accuracy  and  precision  of  positions  determined  using  these  procedures.  The 
adaptability  section  considers  how  compatible  the  procedures  are  with  present 
Coast  Guard  procedures  and  training  levels.  The  "ease  of  use"  section  con¬ 
cerns  itself  with  the  technical  level  required  to  employ  such  a  procedure  and 
the  cumbersomeness  of  the  actual  equipment  required.  The  feedback  section 
discusses  the  time  requirement  between  measurements  and  comparison  of  position 
error  to  a  standard.  The  recording  section  involves  storage  of  data  accumu¬ 
lated  in  the  positioning  effort.  These  include  the  measurements,  the  position 
error  measures,  and  other  data  of  interest. 

8.1  Graphical  Procedures  -  Grid  Diagrams 

8.1.1  Description 


Geographic  coordinates  of  landmarks  in  the  area  and  the 
designated  position  are  found  through  use  of  procedures  presently  being 
employed  by  the  Coast  Guard  (reference  9).  The  alternative  fix  geometries  are 
examined  using  a  chosen  position  error  measure  and  a  priority  listing  is 
created  for  forwarding  to  the  unit.  The  first  priorities  are  used  to  develop 
grid  diagrams  as  is  present  Coast  Guard  procedure.  On  the  grid  plot  are 
established  limits  on  each  line  of  position  which  represent  the  acceptable 
limits  on  the  measurements  used  to  define  the  LOP.  The  lines  of  position  are 
labeled  in  measurement  units.  The  limits  correspond  to  the  selected  set  of 
standards  (remember,  third  measurement  differences  can  be  directly  related  to 
most  position  error  measures). 

The  conning  officer  safely  maneuvers  his  vessel  to  the  loca¬ 
tion  near  where  the  buoy  is  to  be  placed.  Measurements  are  read  continuously 
to  the  conning  officer  for  check  against  the  standard  and  the  grid  diagram. 

The  measurements  are  sounded  in  angular  units  or  in  fractions  of  a  glass  de¬ 
pending  on  the  preferred  procedure.  The  grid  diagram  is  drawn  so  that  each 
line  corresponds  to  the  units  sounded.  The  position  is  approximated  on  the 
grid  diagram  and  corrected  for  eccentric  errors  by  use  of  a  scaled  model  of 
the  ship  on  the  grid  chart.  When  it  is  evident  that  the  sinker  will  drop 
close  enough  to  the  desired  location,  it  is  dropped  and  the  final  measurements 
are  taken.  See  appendix  A  for  other  hydrodetic  procedures  to  be  followed 
before  and  during  positioning  efforts. 

8.1.2  Accuracy 

This  method  is  quite  accurate.  The  eccentric  errors  are 
corrected  using  scaled  models  of  the  positioning  unit. 


‘M 


34 


8.1.3  Adaptabi 1 ity 

Grid  diagrams  are  being  prepared  for  buoys  throughout  the 
Coast  Guard.  It  is  not  difficult  to  calculate  the  maximum  acceptable  measure¬ 
ment  tolerances  for  the  lines  of  position  on  the  diagram  using  equations  in 
the  appendix  of  this  report.  Little  additional  training  is  needed  to  imple¬ 
ment  this  procedure. 

8.1.4  Ease  of  Use 


The  grid  diagram  procedure  requires  very  little  technical 
expertise.  The  computations  are  performed  at  operational  levels  higher  than 
the  positioning  unit  (reference  9).  The  systematic  eccentricity  errors  do, 
however,  involve  a  knowledge  of  both  navigation  and  the  grid  diagram  for 
proper  compensation.  The  preparation  of  grid  diagrams  is  a  one-time  effort 
and  the  paperwork  is  very  manageable. 

8.1.5  Feedback 


The  time  between  the  measurement  taking  and  position  error 
determination  is  dependent  only  on  how  quickly  LOPs  can  be  plotted  on  the  grid 
and  how  quickly  eccentricity  errors  can  be  compensated  for.  For  typical 
situations  involving  three  angles,  feedback  is  on  the  order  of  30  to  90 
seconds. 


8.1.6  Recording 

All  recording  of  data  is  done  manually. 
8.2  Calculator  Based  Procedures 


8.2.1  Description 


Geographic  coordinates  of  the  landmarks  in  the  area  and  the 
designated  position  are  found  throt!  >  use  of  procedures  presently  employed  by 
the  Coast  Guard.  The  alternative  fi  geometries  are  examined  at  either  the 
field  or  district  level  using  a  position  error  measure  chosen  at  either  level 
and  a  priority  listing  is  created.  The  first  priorities  are  used  in  position¬ 
ing  the  aid  if  possible.  If  not  possible,  the  unit  simply  employs  a  fix 
geometry  that  is  possible. 

The  ship  is  safely  maneuvered  to  the  area  around  the  desired 
location  where  measurements  are  taken  and  entered  into  a  programmed  calcula¬ 
tor.  Immediate  output  of  position  error  and  guidance  assist  is  provided  by 
the  calculator.  The  eccentricity  error  compensation  is  a  part  of  the  calcula¬ 
tion.  The  position  error  desired  (or  set  of  position  errors)  is  compared 
against  the  standards  and  the  decision  whether  or  not  to  drop  the  sinker  is 
made.  The  sinker  is  dropped,  final  measurements  (not  restricted  to  three)  are 
made,  and  the  important  information  is  recorded  on  magnetic  tapes  and  paper 
output.  As  part  of  the  positioning  calculations,  an  outlier  detection  and 
rejection  routine  makes  blunder  detection  realizable.  Routine  hydrodetic  pro¬ 
cedures  such  as  discussed  in  appendix  A  can  be  programmed  and  employed  con¬ 
tinuously  or  at  will . 
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8.2.2  Accuracy 

The  calculator  assisted  procedure  allows  use  of  many  measure¬ 
ments  simultaneously,  compensates  for  systematic  eccentricity  errors  and 
detects  obvious  measurement  inconsistencies.  It  is  very  accurate. 

8.2.3  Adaptability 

Very  few  Coast  Guard  personnel  have  any  experience  with 
programmable  calculators.  However,  programs  created  from  equations  in  the 
appendix  require  little  computer-user  interface  once  the  control  data  has  been 
entered  into  the  calculator.  Training  requirements  for  handling  the  control 
data,  caring  for  the  equipment,  and  general  calculator  operation  are  far 
greater  than  those  of  graphical  procedures.  The  method  is  adaptable  for  use 
by  personnel  on  any  size  aid  to  navigation  unit. 

8.2.4  Ease  of  Use 


State-of-the-art  programmable  calculators  have  been  built 
with  user  definable  keys  as  if  designed  for  this  particular  application.  The 
technical  level  required  of  the  operator  is  higher  than  with  other  methods  but 
surely  not  beyond  the  capabilities  of  a  high  school  graduate.  Simple  operat¬ 
ing  instructions  would  make  operation  no  more  difficult  than  normal  ship 
navigation  procedures.  The  hardware  required  to  perform  the  task  must  be 
light,  portable  and  durable.  Oata  loading  and  recording  must  be  kept  to  a 
minimum. 


8.2.5  Feedback 


The  time  between  the  measurement  taking  and  position  error 
determination  is  dependent  on  the  selected  measure  of  position  error.  It  can 
vary  from  seconds  to  many  minutes.  The  most  useful  measures  are  retrievable 
in  seconds  (<7  seconds  on  the  HP-41C  by  Hewlett  Packard). 

8.2.6  Recording 

Recording  of  data  can  be  performed  immediately  upon  comple¬ 
tion  of  the  operation  automatically  in  hard  copy  form. 
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9.0  CONCLUSIONS 


a.  Analysis  of  aid  positioning  data  indicates  that  there  is  room  for 
improvement  in  present  aid  positioning  procedures. 

b.  Present  aid  positioning  efforts  can  be  supplemented  by  analytical 
procedures  employing  either  calculators  or  extended  grid  techniques 

c.  Classifications,  standards  of  accuracy,  and  specifications  for 
geodetic  survey  (reference  25)  are  not  applicable  to  positioning 
floating  aids  but  procedures  are  available  to  achieve  the  same 
objectives  with  a  lower  degree  of  accuracy. 

d.  Planning  and  positioning  models  are  used  for  distinctly  different 
reasons  but  both  are  necessary  to  minimize  possibility  of  error 
either  systematic  or  stochastic  in  nature. 

e.  Exact  geodetic  computations  can  be  approximated  with  acceptable 
accuracy  for  Coast  Guard  use  in  analytical  positioning. 

f.  Mixing  of  measurement  types  poses  no  significant  problem  when  using 
analytical  procedures. 

g.  Detection  of  measurement  outliers  is  feasible  through  use  of  analyt 
ical  procedures. 


The  "A  posteriori  estimate  of  the  reference  measurement  variance"  i 
a  useful  statistic  easily  employed  to  test  fix  quality  through 
analytical  procedures. 


i.  Many  different  mathematical  measures  of  position  error  exist  for 
planning  and  positioning. 
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10.0  RECOMMENDATIONS 

The  following  planning  should  be  performed  prior  to  each  positioning 
evolution. 

1.  Selection  of  the  priority  fix  geometries  and  position  error 
measures  which  provide  the  best  indicators  of  the  accuracy  and 
precision  requirements  for  each  aid  location. 

2.  Determination  of  reasonable  standards  for  the  selected  posi¬ 
tion  error  measures  by  use  of  the  planning  model. 

3.  Specifications  should  be  stated  for  the  number  of  measurements 
needed,  the  check  measurements  required,  and  the  angle  clo¬ 
sures  to  be  performed. 

4.  Note  where  problems  in  positioning  might  be  expected,  such  as 
landmark  descriptions,  landmark  orientations,  and  poor  geome¬ 
tries. 

5.  Updates  of  control  data  on  each  prospective  landmark  must  be 
performed. 

6.  Gradient  diagrams  should  be  prepared  with  all  standards 
readily  visible  in  measurement  units.  Diagrams  should  specify 
scale  model  size  to  be  used  as  a  representation  of  the  posi¬ 
tioning  unit  on  the  grid  diagram.  Diagrams  should  be  prepared 
for  all  top  priority  geometries. 

The  following  hydrodetic  procedures  should  be  considered  as  part  of  the 
normal  positioning  routine  aboard  all  positioning  units  (regardless  of  whether 
or  not  presently  used  analytical  procedures  are  extended). 

1.  Horizon  closure  before  and  after  movements  to  ensure  instru¬ 
ment  accuracy. 

2.  Angle  sum  (difference)  measurement  simultaneous  with  angle 
measurement  when  a  landmark  is  used  for  more  than  one  measure¬ 
ment. 

3.  Observe  only  geodetically  controlled  landmarks  (however,  the 
classification  is  likely  to  be  of  little  importance). 

4.  As  prescribed  in  reference  9,  routine  determination  of  index 
error  and  instrunent/observer  random  error  should  be  performed. 

5.  Scaled  positioning  units  (on  grid  diagrams  for  angle  taker  to 
chain  stopper  compensation)  should  be  used. 

6.  Remove  any  reliance  on  the  three-arm  protractor  and  chart  for 
aid  positioning. 

It  is  recommended  that  the  following  approach  be  taken  with  regards  to 
establishment  and  implementation  of  standards  based  on  analytical  procedures. 
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1.  Where  grid  diagrams  are  being  used,  they  should  be  extended  to 
include  the  third  measurement  difference  standards  discussed 
in  this  report. 

2.  From  information  provided  here  and  in  error  modelling,  stan¬ 
dards  on  the  magnitude  of  the  AP-MPP  vector  magnitude  should 
be  established. 

3.  When  parts  1.  and  2.  are  in  force,  classification  of  aid 
positions  should  be  made  by  use  of  the  R-for-P  position  error 
measure. 

In  addition  to  the  grid  diagram  procedures,  calculator  based  procedures 
should  be  explored  for  practical  areas  of  application.  The  following  approach 
is  reconmended: 

Program  all  applicable  positioning  equations  on  state-of-the-art 
programmable  calculators.  Reinstate  CAP-II  aboard  the  REDWOOD 
(WLM-685)  to  demonstrate  that  calculators  can  be  used  to  ease  the 
burden  of  aid  positioning  as  well  as  make  the  position  determined 
more  accurate  and  precise. 


The  fear  of  mathematics  and  computer  science  will  not  be 
removed  unless  we  subject  personnel  at  all  levels  to  their  use.  Even  if  the 
analytical  procedures  are  placed  at  levels  higher  than  the  field  unit,  the 
field  unit  personnel  should  have  the  computer  capacity  and  computing  knowledge 
to  understand  and  check  the  standards  and  procedures  supplied  by  higher  levels 
of  administration.  Many  field  units  have  personnel  with  the  mathematical  and 
programming  capability  to  expand  on  the  ideas  presented  in  this  report. 

The  results  of  the  following  reports  should  be  reviewed,  understood  and 
integrated  to  compare  with  original  ANPAR  objectives  and  to  form  recommenda¬ 
tions  for  additional  project  efforts. 

1.  Positioning/Error  Model-First  Interim  Report 

2.  Positioning/Error  Model-Error  Sensitivity  Model-Second  Interim 

Report 

3.  Positioning/Error  Model  -Analytical  Positioning  of  Aids  to 
Navigation  (This  Report) 

4.  Sinker  Orop  Error  Analysis 

5.  Off  Station  Buoy  Analysis 

6.  Watch  Circle  Analysis  (unpublished) 
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7.  Position  Accuracy  Report  (unpublished) 

8.  Sextant  Evaluation  Laboratory  Report  (unpublished) 

Institute  a  class  on  "The  Positioning  of  Aids  to  Navigation"  at  the 
Coast  Guard  Training  Center,  Governors  Island.  Direct  the  class  toward  the 
junior  officer,  senior  enlisted  man,  and  cadet  destined  for  service  on  a 
positioning  unit.  Structure  the  course  to  include: 

1.  References 

a.  ATON  Manual  -  Positioning  .COMDTINST  1116500.1,1978 

b.  Hydrographic  Manual,  Umbach,  M.J.,US  Dept  of  Commerce, 1976 

c.  Handbook  on  Grid  Usage  (to  be  published  for  course) 

d.  Chart  No.  1  Nautical  Charts  Symbols  and  Abbreviations 

e.  Pocket  Calculator  Use  in  Positioning  (to  be  published  for 
course) 

f.  Geodesy  for  the  Layman,  US  Dept  of  Commerce.  NOAA,  NOS. 
1977 

g.  Sextant  Adjustment  Manual,  CG  ATON  School,  Governor's 
Island 

h.  Classifications,  Standards  of  Accuracy,  and  General 
Specifications  for  Geodetic  Control  Surveys,  FGCC,  US  Dept 
of  Commerce,  1976. 

i .  Specifications  to  Support  Classifications,  Standards  of 
Accuracy,  and  General  Specifications  for  Geodetic  Control 
Surveys,  FGCC,  US  Dept  of  Commerce,  1975. 

2.  Course  Content 

a.  Grid  Usage 

b.  Calculator  Usage — •*  Analytical  Procedures 

1.  Positioning 

2.  Blunder  Detection  Methods 

3.  Methods  for  correcting  poor 
situations 

c.  Survey  and  Geodesy  Basics 

d.  Sextant  Usage 

e.  Discussion  of  Accuracy  and  Precision 
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GLOSSARY 


Accuracy  -  Degree  of  conformity  with  a  standard.  Accuracy  relates  to  the 

quality  of  a  result.  The  accuracy  attained  in  a  positioning  effort  is  a 
product  of  the  procedures  to  be  followed  in  executing  the  work  and  the 
precision  with  which  those  instructions  are  followed. 

Adjustment  -  The  determination  and  application  of  corrections  corresponding  to 
the  errors  affecting  the  observations,  making  the  observations  consistent 
among  themselves,  and  coordinating  and  correlating  the  derived  data. 
Adjustment  is  comnonly  performed  by  the  method  of  least  squares. 

Administrative  Procedures  -  Accumulation  and  processing  of  positioning  data 
for  ver i f i cat i on ,  1 ega 1  and  planning  purposes. 

Analytical  Planning  -  Employment  of  numerical  computations  in  a  systematic 
procedure  for  deriving  the  expected  results  of  a  positioning  effort  and 
setting  standards  for  the  effort. 

Astronomical  Sighting  -  Observation  of  the  azimuth  of  a  celestial  body. 

Audi t  -  Independent  check  on  the  position  of  aids  to  navigation. 

Azimuth  -  The  horizontal  direction  reckoned  clockwise  from  the  meridian 

plane.  In  the  basic  control  survey,  azimuths  are  measured  clockwise  from 
south.  The  common  procedure  for  navigation  is  to  measure  clockwise  from 
north.  North  is  used  in  this  report. 

Azimuth?  Astro  -  The  data  obtained  through  an  astronomical  sighting.  At  the 
point  of  the  observation,  the  angle  measured  from  the  vertical  plane 
through  the  celestial  pole  to  the  vertical  plane  through  the  observed 
object. 

Azimuth,  Forward  -  For  a  geodesic  line  from  A  to  B,  the  angle  between  the 
tangent  to  the  meridian  at  A  and  the  tangent  to  the  geodesic  line  at  A. 

Azimuth,  Back  -  For  a  geodesic  line  from  A  to  B,  the  angle  between  the  tangent 
to  the  meridian  at  B  and  the  tangent  to  the  geodesic  line  at  B. 

Basel i ne  -  The  side  of  one  of  a  series  of  connected  triangles,  the  length  of 
which  is  measured  with  perscribed  accuracy  and  precision,  and  from  which 
the  lengths  of  the  other  triangle  sides  are  obtained  by  computation. 

Base  Measurement  -  Determination  of  the  length  of  the  baseline  classified 
according  to  the  character  of  the  work  they  are  intended  to  control. 
Probable  error  in  measurement  not  to  be  exceeded  for  the  various  classes 
is  prescribed. 

Bearing  -  The  horizontal  direction  of  one  terrestrial  point  from  another.  It 
is  measured  from  north  clockwise  through  360°. 

Bivariate  Statistics  -  Statistics  of,  relating  to,  or  involving  two  variables. 
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Blunder  -  A  gross  error  or  mistake  resulting  usually  from  stupidity, 
igio  ranee,  or  carelessness. 

Chord  -  A  straight  line  segment  joining  two  points  on  a  curve  ( e . g . ,  circle, 
sphere,  ellipsoid). 

Clarke  Spheroid  of  1866  -  See  spheroid.  The  Clarke  Spheroid  has  a  major 
semi-axis  of  6,378,206.4  meters  and  a  minor  semi-axis  of  6,356,583.8 
meters. 

Class  (survey)  -  The  division  or  rating  of  surveys  based  on  the  overall 
acc  uracy  and  precision  required  of  the  survey. 

Closure,  Error  of  -  The  amount  by  which  a  value  of  a  quantity  obtained  by 
surveying  operations  fails  to  agree  with  another  value  of  the  same 
quantity  held  fixed  from  earlier  determinations  or  with  a  theoretical 
value  of  the  quantity.  The  quantity  may  be  an  angular  measure,  distance 
measure  or  spacial  coordinates. 

Conditions  -  An  equation  which  expresses  exactly  certain  relationships  that 
must  exist  among  related  quantities,  which  are  not  independent  of  one 
another,  exist  a  priori,  and  are  separate  from  relationships  demanded  by 
observati on. 

Confidence  Ellipse  -  One  of  a  family  of  contours  of  equal  probability  density 
of  the  bivariate  probability  density  function  of  a  position  determined  by 
measurements.  Two  types  of  ellipses  are  considered:  (1)  those  which  may 
be  expected,  based  on  the  laws  of  chance  and  the  assumed  values  of 
measurement  random  error  as  used  in  planning,  and  (2)  those  which  result 
by  statistical  inference  from  analysis  of  a  particular  fix. 

Convergence  -  The  angular  drawing  together  of  the  geographic  meridians  in 
passing  from  the  equator  to  the  poles.  For  a  geodesic,  the  azimuth  at 
one  end  differs  from  the  azimuth  at  the  other  end  by  180°  plus  or  minus 
the  amount  of  the  convergence  of  the  meridians  at  the  end  points. 

Conversion  Angle  -  The  angle  between  the  rhumb  line  and  the  great  circle 
between  two  poi nts . 

Designated  (desired,  computed)  Position  -  The  point  specified  by  geographic 
coordinates  that  authorities  have  decided  an  aid  is  to  be  placed. 

Direct  and  Reverse  Sightings  -  A  method  to  reduce  systematic  reading  errors  by 
which  the  telescope  is  rotated  and  turned  such  that  the  readings  are  made 
opposite  on  the  horizontal  plate  (i.e.,  180°  different). 

Pi  recti  on  -  The  position  of  one  point  in  space  relative  to  another  without 

reference  to  the  distance  between  the  points  and  provided  in  terms  of  the 
angular  difference  from  the  reference  direction  (normally  the  reference 
is  north). 

Direction  (survey)  -  Horizontal  angles  at  a  triangulation  station  are  reduced 
to  a  common  initial  and  termed  horizontal  directions. 
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Distance  Angles  -  An  angle  in  a  triangle  opposite  a  side  used  as  a  base  in  the 
solution  of  the  triangle,  or  a  side  whose  length  is  to  be  determined. 

Distance  Resi dual s  -  The  computed  distance  along  a  line  perpendicular  to  a 
line  of  position  between  the  most  probably  position  of  the  observers  and 
the  line  of  position. 

Eccentricity  -  the  ratio  of  the  distance  between  the  center  and  the  focus  of  a 
ellipse  to  the  length  of  its  major  semi-axis. 

Ellipse  -  A  closed  plane  curve  generated  by  a  point  moving  in  such  a  way  that 
the  sums  of  its  distances  from  two  fixed  points  is  a  constant. 

El  1  ipso  id  -  The  three-dimensional  figure  formed  by  rotating  an  ellipse  about 
its  major  or  minor  axes. 

Error,  Modeling  -  A  system  of  postulates,  data,  and  inferences  presented  as  a 
mathematical  description  of  errors  found  in  a  process. 

Error.  Propagation  of  -  involves  obtaining  the  stochastic  characteristics  of 
(functionally)  dependent  variables  given  the  characteristics  of  the 
independent  variables  and  the  functional  relationships  relating  the  two 
sets  of  variables. 

Error,  Random  -  A  chance  error,  unpredictable  in  magnitude  or  sign. 

Error,  Residual  (Residuals)  -  The  difference  between  any  value  of  a  quantity  in 
a  series  of  observations,  corrected  for  known  systematic  errors,  and  the 
value  of  the  quantity  obtained  from  the  combination  or  adjustment  of  that 
series . 

Error,  Standard  -  Sample  standard  deviation  of  the  mean. 

Error,  Systematic  -  An  error  that  is  not  determined  by  chance  but  whose  sign 
and,  to  some  extent,  magnitude  bear  a  fixed  relation  to  some  condition  or 
set  of  conditions. 

Geodetic  Coordinates  -  Quantities  which  define  the  horizontal  position  of  a 
point  on  the  spheroid  of  reference  with  respect  to  the  planes  of  the 
geodetic  equator  (latitude)  and  of  a  selected  geodetic  meridian  (longi¬ 
tude)  . 

Geodetic  Datum  -  Numerical  or  geometrical  quantity  or  set  such  quantities 
which  may  serve  as  a  reference  or  base  for  other  quantities.  Geodetic 
datum  consists  of  5  quantities:  latitude,  longitude,  azimuth  of  some 
geodesic  at  the  point  of  concern  and  the  two  constants  needed  to  define 
the  terrestrial  spheroid. 

Geodetic  Survey  Standards  -  Quantities  specified  by  the  Federal  Geodetic 
Control  Commission  and  published  by  the  U.S.  Department  of  Commerce  to 
ensure  high  precision  and  good  accuracy  in  horizontal  and  vertical 
control  surveys. 


Geodesy  -  The  science  which  treats  mathematically  the  figure  and  size  of  the 
earth. 

Geodesic  (Geodesic  Line)  -  A  line  of  shortest  distance  between  any  two  points 
on  any  mathematically  defined  surface. 

Geoid  -  The  figure  of  the  earth  considered  as  a  mean  sea-level  surface 
extended  continuously  through  the  continents. 

Geometry  (of  a  fix)  -  Consists  of  all  the  lines  of  position  which  represent  the 
measurements  that  compose  the  fix,  their  respective  measurement  standard 
deviations,  and  their  orientation  with  respect  to  each  other. 

Global  Survey  Network  -  Horizontal  control  survey  net:  arcs  of  triangulation, 
sometimes  with  lines  of  traverse,  connected  together  to  form  a  system  of 
loops  or  circuits  extended  over  an  rea.  Global  indicates  that  all 
available  tri angulation  nets  are  considered. 

Gradient  Diagram  (Grid)  -  A  graphical  representation  of  the  computed  lines  of 
position  that  correspond  to  a  specified  geometry  at  an  aid  location. 

Lines  of  position  are  spaced  at  regular  intervals  in  measurement  units 
parallel  to  the  desired  lines  of  position,  thus  forming  a  grid. 

Gradient  Vector  -  Describes  quantitively  by  magnitude  and  direction  the 

transverse  movement  of  a  line  of  position  for  a  change  of  one  measurement 
unit. 

Graphical  Plotting  Procedure  -  In  resection,  the  use  of  the  three-arm 
protractor  and  chart  in  position  determination. 

Horizontal  Control  Data  -  Data  established  from  surveys  which  are  used  with 
measurements  to  define  the  position  represented  by  the  values  of  the 
measurements . 

Horizontal  Directions  -  See  Directions. 

Hydrodetic  Procedures  -  Procedures  based  on  as  many  geodetic  survey  standards 
and  methods  as  possible  but  which  are  performed  in  the  marine  environment. 

Instrument  Specification  -  Indicates  the  instrument  resolution  required  for 
conduct  of  an  acceptable  survey  of  a  given  order  and  class. 

Inverse  -  The  computation  of  the  length  and  forward  and  back  azimuths  of  a 

geodesic  by  computation  based  on  the  known  geodetic  positions  of  the  ends 
of  the  line. 

Least  Squares  Principle  -  A  mathematical  method  of  determining  the  most 
probable  values  of  a  series  of  quantities  from  a  set  of  observations 
greater  in  number  than  are  necessary  to  determine  those  quantities. 

Line  of  Position  -  A  line  on  some  point  of  which  an  observer  may  be  presumed 
to  be  located,  as  a  result  of  observation  or  measurement. 
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LOP  Variance  -  The 
coordinate  of 
square  of  the 
determine  the 


square  of  the  standard  deviation  of  the  transverse  spacial 
a  line  of  position.  The  LOP  variance  is  the  product  of  the 
gradient  magnitude  and  the  measurement  variance  used  to 
line  of  position. 


Lower  Order  Stations  -  In  a  survey,  the  stations  which  are  not  principal 
stations  (less  than  3rd  order  stations). 

Major  Semi -Ax is  -  The  distance  from  the  center  of  an  ellipse  to  the  ellipse 
along  the  longer  axis  of  the  ellipse. 

Marine  Geodetic  Standards  (Hydrodetic  Standards)  -  Standards  applicable  to 
survey  work  at  sea. 

Marine  Geodesy  -  Geodesy  applied  to  the  environment  of  the  sea  with 
intentions  of  extending  geodetic  control  to  that  environment. 

Measurement  Combination  -  A  subset  of  the  measurements  made  in  determining  a 
position. 


Measurement  Variance  -  The  square  of  the  standard  deviation  of  measurement. 

Meri di an  -  A  north-south  line  from  which  longitudes  and  azimuths  are  reckoned; 
or  a  plane,  normal  to  the  geoid  or  spheroid  defining  such  a  line. 

Minor  Semi -Axis  -  The  distance  from  the  center  of  an  ellipse  to  the  ellipse 
along  the  shorter  axis  of  the  ellipse. 

Misclosure  -  See  Closure,  error  of. 

NAD  27  -  North  American  Datum  of  1927.  The  geodetic  datum  which  is  defined  by 
the  geographic  position  of  tri angulation  station  Meades  Ranch  and  the 
azimuth  from  that  station  to  station  Waldo,  on  the  Clarke  Spheroid  of 
1866. 


Number  of  Positions  -  Specifies  the  number  of  locations  on  the  horizontal 
circle  of  a  direction  theodolite  to  be  used  for  the  observation  on  the 
initial  station  of  a  series  of  stations  which  are  to  be  observed  on. 

Operational  Procedures  -  Those  procedures  used  in  the  practical  application  of 
~~  the  principles  and  processes  of  positioning. 

Order  (Survey)  -  A  category  describing  the  quality  of  a  survey. 

Outlier  -  A  measurement  that  is  far  from  the  mean  due  to  chance  or  blunder. 

P-in-R  -  The  probability  mass,  P,  contained  by  a  circular  region  of  radius,  R, 
centered  on  the  desired  position. 

PI anninq  Model  -  A  mathematical  model  of  the  positioning  process  used  to 
assess  the  expected  results  of  specific  positioning  efforts,  and  to 
establish  standards  for  positioning,  to  study  problems  incurred  during 
positioning  efforts. 


47 


Position  -  A  location  on  the  horizontal  circle  of  a  direction  theodolite  to  be 
used  for  the  observation  of  the  initial  station  of  a  series  of  stations 
which  are  to  be  observed  on. 

Position  Error  Measure  -  A  quantity  that  represent  the  accuracy,  precision  or 
both  accuracy  and  precision  of  determined  position;  to  be  used  in  setting 
standards  and  as  a  measure  of  positioning  success. 

Precision  -  The  degree  of  refinement  with  which  an  operation  is  performed  or  a 
measurement  stated.  Usually  represented  by  the  standard  deviation  of  a 
set  of  measurements  of  the  same  quantity. 

Prime  Vertical  -  Vertical  circle  perpendicular  to  the  plane  of  the  celestial 
meridian.  The  plane  of  the  prime  vertical  cuts  the  horizon  in  the  east 
and  west  points. 

Principal  Station  -  A  station  through  which  basic  data  are  carried  in  the 
extension  of  a  survey  system.  The  principal  station  is  a  higher  order 
station  relative  to  those  stations  whose  purpose  is  limited  to  the 
control  of  local  surveys. 

Probability  Expectation  Approach  -  Any  approach  which  uses  or  assumes  a  known 
probabil ity  distribution  for  each  step. 

Propagation  of  Error  -  See  Error,  propagation  of. 

Propagation  Velocity  -  The  speed  at  which  electromagnetic  radiation  passes 
through  a  medium.  Usually  termed  the  speed  of  light. 

Published  Data  -  The  geodetic  datum  of  all  control  stations,  formalized  for  use 
as  a  reference  base  in  future  surveys. 

Radius  of  Curvature  -  The  reciprocal  of  the  curvature  of  a  line. 

Recommended  Spacing  -  The  specified  distance  between  principal  stations  in  a 
survey  network. 

Reconnaissance  -  Preliminary  survey  to  gain  information  and  to  plan  a  geodetic 
survey. 

Recording  -  Record  keeping  of  data  for  use  in  planning,  research  and  for  legal 
purposes.  Recording  includes,  but  is  not  limited  to,  magnetic  storage  of 
data. 

Reference  Measurement  Variance  -  An  arbitrary  constant  with  arbitrary 

dimensions,  usually  selected  as  the  best  known  estimate  of  measurement 
variance  of  the  observer-instrument  combination. 

Reference  Position  of  the  Angle  Takers  (AT)  -  In  many  cases,  measurements  can 
not  be  taken  from  the  exact  same  location  due  to  intervisibility  condi¬ 
tions.  The  point  selected  in  the  region  where  measurements  are  being 
taken  to  define  displacement  vectors  for  each  observer's  location  is  the 
reference  position  of  the  angle  takers.  When  the  least  squares  principle 
is  employed,  the  actual  position  determined  is  AT,  and  a  translation  to 
the  point  of  sinker  drop  is  needed  to  define  the  MPP. 
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Rejection  Limit  -  A  specification  on  the  difference  between  a  measurement  and 
the  mean  of  the  set  of  measurements  of  which  it  is  an  element. 

Residual  Error  Averaging  -  The  mathematical  operation  of  spreading  a  residual 
error  to  the  measurements  from  which  the  residuals  are  calculated  so  that 
the  measurement  set  is  consistent  with  theory. 

Resecti on  -  The  determination  of  the  horizontal  position  of  a  survey  station 
by  observed  directions  from  the  station  to  points  of  known  position. 

Residuals  -  See  error,  residual. 

Resolution,  Instrument  -  The  smallest  increment  of  change  that  an  instrument 
and  observer  can  measure. 

R-for-P  -  The  radius,  R,  of  a  circle  centered  on  the  desired  position  that  is 
required  to  contain,  at  least,  a  desired  probability  mass,  P. 

Side  Check  -  A  specification  that  requires  agreement  in  triangle  side  lengths 
as  computed  in  various  chains. 

Sphere  -  The  three-dimensional  solid  figure  which  consists  of  the  set  of  all 
points  equidistant  from  a  point  constituting  its  center. 

Spherical  Excess  -  The  amount  by  which  the  sum  of  the  three  angles  of  a 
triangle  on  a  sphere  exceeds  180°  (also  assumed  for  spheroids). 

Spheroi d  -  Any  figure  differing  but  little  from  a  sphere.  In  geodesy,  a 

mathematical  figure  closely  approximating  the  geoid  in  form  and  size,  and 
used  as  a  surface  of  reference  for  geodetic  surveys  (normally  an 
ellipsoid  of  revolution). 

Standards  -  An  exact  value,  or  concept  thereof  established  by  authority, 

custom,  or  common  consent,  to  serve  as  a  model  or  rule  in  the  measurement 
of  quantity,  or  in  the  establishment  of  a  practice  or  procedure. 

Standards  of  Accuracy  -  Standards  established  on  the  degree  of  perfection 
obtained  in  a  survey.  They  serve  as  benchmarks  for  determining  the 
quality  of  a  result. 

Station  -  A  definite  point  on  the  earth,  whose  location  has  been  determined  by 
surveying  methods. 

Strength  of  Figure  -  Expresses  the  comparative  precision  of  computed  lengths 
in  a  tri angulation  net  as  determined  by  the  size  of  the  angles,  the 
number  of  conditions  to  be  satisfied,  the  distribution  of  base  lines,  and 
lengths  determined  in  previous  adjustments.  An  expression  of  relative 
strength. 

Systematic  Error  Tendency  -  A  vector  quantity  describing  how  much  displacement 
the  position  determined  by  a  set  of  measurements  changes  when  all  mea¬ 
surements  are  changed  one  unit. 

Tangent  -  Touching  a  curved  surface  at  only  one  point. 
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Traverse  -  A  sequence  of  lengths  and  directions  of  lines  between  points  on  the 
earth,  obtained  by  or  from  field  measurements  and  used  in  determining 
positions  of  the  points. 

Triangle  Closure  -  The  amount  by  which  the  sum  of  the  three  observed  angles  of 
a  triangle  fails  to  equal  exactly  180°  plus  the  spherical  excess  of  the 
triangle. 

Tri angulation  -  A  method  of  surveying  in  which  the  stations  are  points  on  the 
ground  at  the  vertices  of  a  chain  or  network  of  triangles,  whose  angles 
are  observed  instrumentally  and  whose  sides  are  determined  by  computation 
from  selected  triangle  sides  called  base  lines,  the  lengths  of  which  are 
obtained  from  direct  measurements  on  the  ground. 

Tri angulation  Net  -  See  Global  survey  network 

Tri later at  ion  -  A  method  of  extending  horizontal  control  where  the  sides  of 
triangles  are  measured  rather  than  the  angles  as  in  tri angul ation. 

Truth  -  Quantity  accepted  to  be  perfect.  Measurement  results  are  compared  to 
truth  to  determine  accuracy.  Truth  in  survey  is  provided  by  higher  order 
stations. 

Tentative  Truth  -  Quantities  may  be  determined  by  higher  order  procedures  and 
specifications  but  until  they  are  adjusted  to  fit  the  triangulation  net, 
they  are  only  tentative  truth. 

Verification  -  Process  of  checking  computations  to  ensure  compliance  with 
applicable  standards  and  freedom  of  correctable  error. 
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APPENDIX  A 


MATHEMATICS  OF  GEODETIC  SURVEY 


A.l  SURVEY  PROCEDURES 

Geodetic  survey  specifications  and  standards  of  accuracy  are  discussed  in 
this  appendix.  Each  section  is  titled  by  an  appropriate  geodetic  survey  term, 
those  introduced  in  section  3.2. 

A. 1.1  Strength  of  Figure  (reference  6) 

The  strength  of  figure  is  derived  from  that  portion  of  the  formula 
for  probable  error  of  a  triangle  side  which  is  independent  of  the  accuracy  of 
the  observations,  as  follows: 


~N'd— NdNC  Xl  [§A  +  SA  SB+  Sb_  (A-1) 

in  which  Nj  and  Nc  are  the  numbers  of  directions  observed  and  of  condi - 
tions  to  be  satisfied,  and  S  ^  and  8g  are  the  rates  of  change  in  the 
sines  of  the  distance  angles  A  and  B,  usually  expressed  by  the  differences  of 
the  logarithms  of  the  sines  for  a  difference  of  1  second  in  the  angles,  the 
sixth  decimal  being  the  unit  place.  By  summing  up  the  values  obtained  by 
formula  for  the  simple  figures  composing  a  triangulation  net,  the  strength  of 
figure  for  the  net  can  be  obtained.  As  a  triangulation  net  is  usually  com¬ 
posed  of  several  different  systems  of  simple  figures,  comparable  values  of 
different  systems  are  obtained,  and  the  strongest  route  can  then  be  selected 
through  which  to  carry  a  computation  of  length.  Reconnaissance  for  a  triangu¬ 
lation  net  is  usually  executed  under  instructions  which  specify  limiting 
values  for  the  strength  of  figure  for  the  best  and  second-best  chains  of  tri¬ 
angles  between  adjacent  base  lines,  the  sites  for  stations  and  for  baselines 
being  selected  accordingly. 

Strength  of  figure  is  not  applicable  to  buoy  positioning  in  tha  Coast 
Guard.  It  is,  however,  a  necessary  element  of  fixed  aid  position  computations. 

For  planning  purposes,  other  figures  of  merit  must  be  used  prior  to  buoy  positioning. 
The  figures  of  merit  suggested  in  this  report  are  described  in  Appendix  B. 

A. 1.2  Recommended  Spacing  of  Principal  Stations 

Principal  stations  are  the  link  between  local  geodetic  surveys  and 
the  global  survey  network.  In  geodetic  survey  the  link  is  made  a  specified 
number  of  times  per  distance  surveyed  to  ensure  accuracy  and  consistency  of 
the  system. 


Recommended  spacing  is  not  applicable  to  the  CG  except  *hen 
surveying  fixed  aids  to  navigation.  For  floating  aids,  the  ’ ■ n*  to  the  geo¬ 
detic  survey  network  is  made  through  the  landmarks  used  to  post-on  *■:. 
Third-Order  Class  II  or  better  stations  are  required  for  s-qnt  nq  •>*  -i.v  - 
tal  angles.  Error  modeling  in  reference  2  provided  pre’immary  ev  •  :e"  >  “u* 
large  errors  (greater  than  target  region  dimension:)  in  *  ’an.iman*  •  •  "it**-. 
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cause  significant  position  error.  In  cases  where  insufficient  Third-Order 
Class  II  or  better  landmarks  are  available  for  sighting,  lower  order  landmarks 
are  acceptable  for  sighting  (with  small  chance  of  causing  position  error) 
until  Third-Order  Class  II  or  better  landmarks  are  made  available.  In  no  way 
does  sighting  Third-Order  Class  II  landmarks  allow  classification  of  the 
surveyed  position  to  be  Third-Order  Class  II. 

A. 1.3  Base  Measurement 


The  methods  of  triangulation  are  quite  demanding.  They  require  a 
considerable  number  of  sightings  with  a  minimum  amount  of  distance  measure¬ 
ment.  The  triangles  are  developed  into  a  net  of  interconnected  figures,  and 
certain  lines,  called  base  lines,  must  be  measured  in  order  to  compute  the 
lengths  of  other  sides  in  the  net.  The  base  lines  must  be  measured  with 
extreme  precision,  since  errors  propagated  through  the  system  originate  with 
those  measured  lines. 

The  precision  with  which  a  base  line  must  be  measured  is  specified 
for  the  various  classes  of  survey  by  a  ratio  of  standard  error  of  the  mean  of 
base  line  measurements  to  the  length  of  the  base  line.  The  standard  error  of 
the  mean,  cr  m,  is  computed  by: 


"m  =  ^  n"(?-lT  (A"2) 

where  v  is  a  difference  between  a  measured  length  and  the  mean  of  all  measured 
lengths,  and  n  is  the  number  of  measurements.  For  Third-Order  Class  II  tri- 
angulation,  the  base  measurements  must  be  precise  to  1  part  in  250,000. 

Base  measurement  is  not  easily  adaptable  to  buoy  positioning  but 
it  is,  of  course,  appropriate  for  fixed-aid  positioning.  Precise  base  mea¬ 
surement  decreases  errors  at  their  origin  before  they  are  propagated  through 
the  system.  A  similar  goal  in  resection  methods  is  to  decrease  the  number  and 
size  of  systematic  errors  in  horizontal  angle  measurement  prior  to  angle  mea¬ 
surement.  One  way  to  accomplish  this  is  to  close  the  horizon  before  and  after 
ship  movements.  Assuming  a  priori  the  random  error  of  a  sextant-observer 
pair  acceptable  ranges  can  be  computed  for  horizon  closure.  The  standard 
deviation  of  the  sum  of  n  angles  of  equal  precision  that  close  the  horizon  is 
V^- times  the  standard  deviation  of  each  angular  measurement. 

Ninety-five  percent  of  the  sums  of  n  such  measurements  are 
expected  to  fall  within  two  standard  deviations  of  360°  if  all  measurements 
are  free  of  systematic  errors.  It  is  reasonably  safe  to  assume  a  systematic 
error  or  mistake  exists  if  the  sum  falls  outside  of  these  limits. 

A.  1.4  Horizontal  Directions 


The  instrument  used,  the  number  of  positions,  and  the  rejection 
limit  from  the  mean  are  considered  under  horizontal  directions. 


Resolution  is  the  smallest  increment  of  change  that  an  instrument 
and  observer  can  detect.  For  geodetic  survey  instruments,  the  resolution  is 


specified  in  tenths  of  seconds  for  First  and  Second  Order  surveys  and  as  one 
second  for  Third  Order  surveys.  Normally,  sextants  are  graduated  in  tenths  of 
minutes.  Both  systematic  and  random  error  of  the  average  sextant  are  on  the 
order  of  minutes  (reference  30),  therefore,  specifications  for  resolution  are 
not  critical  for  sextants  and  the  normal  tenth-of-a-minute  reso-lution  is 
adequate.  However,  it  is  important  to  specify  acceptable  random  error  in 
angular  measurement.  This  is  already  prescribed  in  reference  9. 

The  number  of  positions  specification  requires  the  number  of 
sighting  groups  required  for  measurement  of  each  angle.  This  specification 
requires  that  the  instrument  be  stationary  between  measurements  which  is,  of 
course,  impossible  when  using  the  sextant  aboard  a  floating  unit.  Each  sight¬ 
ing  group  requires  four  separate  sightings.  The  angle  must  be  turned  twice 
and  the  telescope  must  be  dumped  twice  for  direct  and  reverse  sightings.  For 
Third  Order  Class  II  stations  two  positions  require  a  total  of  8  sightings  . 

Similar  redundancy  requirements  for  resection  methods  can  be  made 
by  designating  the  number  of  simultaneous  (or  nearly  so)  measurements  required 
to  determine  a  position.  Current  requirements  (reference  9)  are  for  a  minimum 
of  three  horizontal  angle  measurements. 

Three  measurements  provides  some  redundancy  in  the  position  deter¬ 
mination  but  by  no  means  does  it  allow  Third-Order  Class  II  specifications  to 
be  achieved.  To  satisfy  Third-Order  Class  II  specifications,  eight  simulta¬ 
neous  measurements  of  each  of  three  different  horizontal  angles  would  be 
required  by  six  different  12-man  positioning  teams  (six  measurements  needed  to 
reach  1-second  resolution  from  6-second  instrument,  24  men  to  assign  four  per 
angle).  Thus,  144  simultaneous  measurements  would  be  taken.  The  folly  of 
this  process  confirms  that  Coast  Guard  specifications  and  standards  of 
accuracy  are  required;  the  geodetic  survey  standards  are  not  achievable  from 
floating  platforms. 

The  rejection  limit  specification  for  geodetic  survey  applies  only 
when  many  measurements  of  the  same  quantity  are  made  at  different  positions  on 
the  horizontal  circle  and  the  mean  of  the  measurements  at  each  of  the  posi¬ 
tions  is  calculated.  In  turn,  the  mean  of  the  resulting  position  means  is 
calculated.  Each  position  mean  is  compared  to  the  overall  mean  and  if  any 
differs  by  more  than  a  specified  amount  from  that  mean,  the  set  at  that 
position  is  rejected  and  repeated.  To  explain  what  this  accomplishes  mathe¬ 
matically,  let  p  represent  the  number  of  positions,  let  r  be  the  rejection 
limit,  let  m-j  be  the  mean  of  measurements  at  the  i*h  position,  let  m  be 
the  overall  mean.  The  requirement  is  that: 


|  m-j  -  m |  £  r  for  i  from  1  to  p  (A-3) 

This,  in  effect,  is  placing  an  upper  limit  on  the  standard  error  of  the  mean. 
Using  equation  (A-2)  with  the  maximum  acceptable  values  for  differences  from 
the  mean: 
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m  3 


r 

\ZFT 
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which  equals  5  seconds  for  Third-Order  Class  II  survey.  Detection  and  rejec¬ 
tion  methods  applicable  to  this  report  are  found  in  appendix  E. 

A. 1.5  Triangle  Closure 


In  geodetic  survey,  whenever  a  triangle  is  formed  within  the 
triangulation  net,  the  sum  of  the  three  angles  of  the  triangle  must  be  within 
a  specified  range  of  180°  plus  spherical  excess.  The  triangle  closure  test 
is  the  simplest  available  in  the  field  to  ascertain  the  accuracy  of  triangula¬ 
tion  observations.  There  are  specifications  for  both  the  average  and  "maximum 
seldom  to  exceed"  triangle  closures  for  a  tri angul ation  net.  The  effect  of 
this  test  is  the  same  as  that  of  closing  the  horizon  with  horizontal  angle 
measurements. 

The  triangle  closure  specification  is  not  adaptable  to  resection 
methods.  A  viable  alternative  is  to  measure  angles  which  are  geometrically 
the  sum  or  difference  of  other  simultaneously  measured  horizontal  angles. 

With  an  a  priori  accepted  value  of  measurement  standard  deviation,  the  dif¬ 
ference  between  the  sum  (difference)  of  two  measurements  and  the  measured 
horizontal  angle  can  be  compared  statistically  as  a  check  for  measurement 
accuracy.  The  mathematics  of  this  comparison  are  similar  to  that  of  horizon 
closure.  The  checking  difference  between  the  check  measurement  and  the  sum 
(difference)  of  the  checked  measurements  should  be  zero.  The  standard  devia¬ 
tion  of  the  checking  difference  is  the  square  root  of  the  sum  of  the  variances 
of  the  component  measurements.  The  acceptable  difference  in  the  check  angle 
should  be  limited  by  a  multiple  of  the  standard  deviation  of  the  checking 
difference  of  the  component  measurements .  For  example,  two  measurements  are 
made  using  a  common  landmark,  each  with  a  standard  deviation  of  five  minutes. 
The  check  measurement  is  made  with  the  same  precision  on  the  sum  (difference) 
of  the  two  measurements.  The  checking  difference  between  the  sum  (difference) 
of  the  two  measurements  and  the  check  measurement  should  be  within  17  minutes 
(2x5\/T"minutes)  of  zero  95  percent  of  the  time.  It  is  safe  to  question  the 
accuracy  of  all  measurements  if  this  is  not  the  case.  Residual  error  aver¬ 
aging  can  be  employed  if  the  measured  angles  differ  by  less  than  a  prescribed 
standard. 


A. 1.6  Side  Checks 


In  geodetic  survey  procedures,  side  checks  are  made  periodically 
in  the  triangulation  net  as  the  best  check  on  accuracy  in  the  field.  Side 
checks  are  comparisons  of  common  sides  of  triangles  determined  through  dif¬ 
ferent  chains  in  the  tri angul ati on  net. 

The  side  check  is  not  adaptable  to  resection  methods.  A  simple, 
but  not  strong  (reference  2)  substitute  is  to  specify  limits  on  angles,  bear¬ 
ings,  ranges,  and  time  difference  measurements  as  independent  checks  on  the 
position.  The  specified  acceptable  difference  in  a  measurement  from  that 
expected  should  be  in  a  distance  dimension.  This  requires  a  conversion  to 
measurement  units  for  the  instrument  used  (divide  by  gradient). 

The  specified  limits  are  determined  by  the  random  error  associated 

with  the  instrunent. 
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A. 1.7  Closure 


The  accuracy  of  a  fix  is  the  final  determination  for  classifica¬ 
tion.  If  all  of  the  standards  and  specifications  are  complied  with,  the  final 
measure  of  accuracy  should  be  readily  attained.  The  final  measure  in  geodetic 
survey  is  the  ratio  of  the  error  in  the  final  position  to  the  length  of  a  side 
determined  by  previously  performed  higher  order  surveys.  The  acceptable 
ratios  have  been  determined  empirically. 

For  resection  on  floating  platforms  the  same  check  with  higher 
order  stations  is  not  available  (i.e.,  there  is  no  previously  validated  sta¬ 
tion  to  compare  computation  results  with).  An  alternative  is  to  compare  the 
determined  position  to  the  designated  position.  The  final  classification  is 
made  by  consideration  of  both  the  accuracy  and  precision  (defined  harmonious¬ 
ly)  statistics  of  the  determined  position  (those  discussed  in  section  6.0, 
STANDARDS). 
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A. 2  GEODESICS 


A. 2.1  Length  of  Geodesics 

The  equation  of  the  reference  ellipsoid  in  a  Cartesian  coordinate 

system  is: 


x2  +  y2  +  =  a2  { A- 5 ) 

l-e^ 

where  a  is  the  semi -major  axis  (center  to  equator)  and  e^  is  the  square  of 
the  eccentricity.  Let  the  positive  z-axis  pass  through  the  north  poles,  the 
positive  x-axis  pass  through  0°  longitude,  and  the  positive  y-axis  pass 
through  90°E  longitude.  For  the  Clarke  Spheroid  of  1866,  a  =  6,378,206.4 
meters  and  e  =  8.2271854  x  10*2. 

A  geodesic  is  the  shortest  line  connecting  two  points  on  the 
spheroid.  Define  two  points  Pq(<?!>i,  A  i)  and  P2  2  *  A  2 )  on  The 
surface  of  the  spheroids.  <f>  and  A  are  the  conventional  latitudes  and  longi¬ 
tudes  of  the  points  respectively,  as  they  are  published  on  nautical  charts  and 
as  geodetic  control  data.  They  are  also  called  geodetic  or  geographic  coordi¬ 
nates  .  West  longitudes  are  negative.  The  azimuth  of  P2  from  P]  is  called 
the  forward  azimuth  and  is  determined  clockwise  from  geodetic  north  in  opposi¬ 
tion  to  the  conventional  geodesy  reference  clockwise  from  geodetic  south.  The 
length  of  the  chord  connecting  P]  and  P2  is  found  by  the  Pythagorean 
relation  in  the  Cartesian  coordinate  system.  Pi(^>q,  a  -j)  and  P2(<>2»  x2)  are 
converted  to  X,  Y,  Z  coordinates  by: 


X  =  a  cos  <t>  p  cos  A 
Y  =  a  cos  4>p  sin  A 
Z  a  sin  <f>p 


(A-6) 


where  <t>p  is  the  parametric  latitude.  Parametric  and  geodetic  latitude 
differ,  as  follows:  the  parametric  latitude  is  the  angle  at  the  center  of  the 
ellipsoid  between  the  radius  vector  to  the  point  of  interest  and  the  equa¬ 
torial  plane.  The  geodetic  latitude  is  the  angle  between  the  normal  to  the 
surface  at  the  point  and  the  equitorial  plane.  The  two  latitudes  are  related 
by: 


tan 


tan  4> 


( A- 7 ) 


The  points  Pi  and  P2  are  now  in 
length  of  the  chord  between  the 


defined  by  P  i(  X^ ,  Y] ,  Zi )  and  P  2  ( x2 » Y2  •*  z2 ) 
two  points  is:  ^ 


The 


D  =  ((  AX)2  +  (  a Y)2  +  (  aZ)2)V2 


(A -8) 
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K  -M 


where 


AX  =  X2  *  X] 

AY  =  Y2  -  Y}  (A-9) 

A Z  a  Z2  Z! 


The  length  D  is  the  geodesic  approximation  GIC  of  section  3. 3. 1.3. 

The  radius  of  curvature  of  the  geodesic  is  used  to  define  a  sphere 
that  approximates  the  spheroid  at  P}  and/or  P2.  The  radius  of  curvature 
is  found  at  one  end  of  the  geodesic  for  approximation  GL1  of  section  3. 3. 1.2 
and  averaged  for  both  ends  of  the  geodesic  for  approximation  GL2  of  section 
3. 3. 1.1. 

The  radius  of  curvature  of  the  geodesic  is  found  at  some  point  on 
the  geodesic  by  the  calculus  ( reference  7).  The  curvature  depends  on  the 
azimuth  of  the  geodesic  at  the  point  and  is  found  tr igonometr ically  by  combin¬ 
ing  the  radius  of  curvature  of  the  meridian  at  the  point  of  interest  and  the 
radius  of  curvature  of  prime  vertical  at  the  point.  The  radius  of  curvature 
of  the  meridian,  M,  depends  on  the  geodetic  latitude  as: 


M  = 


a(l-e  )? 

( l-e2sin2</>)  ^ 


(A- 10) 


The  radius  of  curvature  of  the  prime  vertical,  N,  is  also  a  function  of  geo¬ 
detic  latitude: 


N= — 7  -  -;-)/?  (A-n) 

(l-e^sin2^)  ' 

The  radius  of  curvature  of  the  geodesic  at  any  point  on  the  geodesic  is  found 
by: 


N  cos +  M  sin< 


(A- 12) 


where  a  is  the  azimuth  of  the  geodesic  at  the  point  of  interest.  For  Gil, 
Ra  is  found  using  the  forward  azimuth  at  P].  For  GL2,  R„  is  found  using 
the  forward  azimuth  at  P}  and  the  back  azimuth  at  P2  and  arithmetically 
averaging  the  two. 

The  radius  of  curvature  of  the  geodesic  (for  Gil  or  GL2)  is  now 
used  to  define  a  sphere  of  radius  R  equal  to  the  radius  of  curvature.  The 
length  of  the  chord,  0,  between  the  two  points  Pi  and  P2  and  the  radius, 

R,  are  used  to  approximate  the  true  length  of  the  geodesic  on  the  reference 
ellipsoid.  From  figure  A- 1 : 
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where  the  inverse  sine  is  calculated  in  radians. 


(A- 13) 


FIGURE  A- 1 


A. 2. 2  Azimuth  of  Geodesics 

The  vector  from  P]  to  p£  is: 

D  =  Jl  X?  +  ^  Yj  +  AZk 


(A- 14) 


where  i,  j  and  k  are  unit  vectors  in  the  x,  y,  and  z  directions,  respective- 
ly.  To  find  the  azimuth,  p  ,  of  the  geodesic  at  P ■] ,  it  is  desired  to  find 
the  angle  between  a  unit  vector  tangent  to  the  surface  at  Pi  in  the  direc¬ 
tion  of  geodetic  north  and  a  vector  defined  by  projecting  D  onto  a  plane 
tangent  to  the  spheroid  at  P].  Define  the  following  unit  vectors  at  Pr 


ns  = 


+  Yj 

*1 

fJi_) 

U-e2/ 

A 

k 

(xf  ♦  Y*  ♦  | 

f  z,  \2> 

U-e2/  ) 

r 

(A- 15) 


which  is  normal  to  spheroid  at  X],  Y],  Z] . 

?  ~Y i i  +  X]j 

'  (x* :  v? , 1/2 

which  is  in  the  direction  of  increasing  longitude. 


(A- 16) 
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<t>  =  ns  x  a 


(A- 17) 


which  is  in  the  direction  of  increasing  latitude. 

The  vector  D  is  projected  onto  the  plane  defined  by  A  and  The  projected 
vector,  IJ|j,  has  components  in  the  %  and  <£  direction  as  follows: 


1)£  =  (0  '  A)  a  (A- 18) 

DJ-  (0  •  $  )  £  (A- 19) 

and 

0p=D^  +  (A- 20) 

The  angle,  «  ,  between*Dp  and  £  is  used  to  find  the  desired  azimuth,  a  is 
found  by. 


( t 


and  the  forward  azimuth  is 


( A— 2 1 ) 


j  360  - 

I 


A  2  <  A  -j 
A 2  ^  *1 


and  the  back  azimuth  (calculated  similarly  as  /Ip)  at  P2  is 


(A-22) 
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B  = 


360  -  a 

a 


<  A 


<  A 


(A- 23) 
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APPENDIX  B 


MATHEMATICS  OF  PLANNING 


The  least-squares  principle  is  employed  to  define  the  mathematical  planning 
model  from  which  expected  results  can  be  obtained.  The  principle  is  discussed 
extensively  elsewhere  (references  2,  9,  13,  14)  and  leads  to  the  following 
matrix  equation: 

X  =  (At  W  A)"1  At  W  L  (B-l ) 


where 


— 

- 

" 

dm] 

3  m] 

Am-j 

7T 

0y 

a  m2 

* 

L  = 

* 

3mn 

3mn 

Amn 

dx 

0y 

- 

±x  =  x  component  (east)  of  vector  from  AP  to  AT 
Ay  =  y  component  (north)  of  vector  from  AP  to  AT 
Am1  =  aQi  -  aci 

"q2  =  arbitrary  reference  variance 

=  ith  measurement  (function  of  x  and  y  coordinates) 

“oi  =  ^  observed  measurement  corrected  for  known  systematic  errors 

aci  =  ith  computed  angle  at  desired  position 

ly,  =  coveriance  matrix  of  observations  (diagonal) 

This  matrix  equation  is  equivalent  to  a  simpler  set  of  equations  for  plan¬ 
ning.  The  following  information  is  needed  for  each  line  of  position  and  is 
providied  by  equations  in  section  4.2: 
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a.  The  positive  gradient  magnitudes, 

b.  The  positive  gradient  directions, 

c.  The  assumed  measurement  random  errors,  <j  \ ,  for  each  instrument 
type. 

The  simpler  set  of  equations  is  used  to  define  position  error  measures  needed 
(1)  for  fix  geometry  selection  criteria,  (2)  to  determine  expected  results,  and 
(3)  for  standard  setting. 

The  simpler  set  of  equations  is  derived  as  follows: 

The  origin  of  the  coordinate  system  is  the  desired  position  of  the  aid;  the 
+y-axis=North,  and  +x-axis=East.  Passing  through  the  origin  with  positive 
gradient  directions,  y  are  the  available  lines  of  position.  The  angles 
y  are  w.r.t.  geodetic  north.  Select  any  n  element  subset  (n>l)  of  lines 
of  position  from  those  available  (of  course  the  y -j  must  be  different  for 
n-2). 


The  matrix  A  is  equivalent  (figure  8-1)  to: 
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The  least  squares  equation  now  can  be  expressed  as:  (derivation  for  two 
measurement  case  for  clarity). 
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with 


Now 


or 


wi  i 


and  <rg2  arbitrarily  in  min2 
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The  AP-to-AT  vector  is  defined  by  ax  and  Ay. 


V  =  ax  T  +  Ay  j' 


(B-6) 


A  A 

where  i  and  j  are  unit  vectors  in  the  position  x  and  y  directions,  respective¬ 
ly. 

The  direction  /3At,  w.r.t  north,  of  V  is, 

ax  >  0  ax  <  0 


Ay  >  0 


AT  = 


Ay  <  0 


t--'  -if 

3600  +  tan-l  1 

180°+  tan-1  — 
Ay 

1800  +  tan-1 

Ay 

0°  <  £AT  <  3600  (B-6) 


In  planning,  the  ^  matrix  is  initially  defined  with  all  zero  elements  because 
all  lines  of  position  are  defined  to  pass  through  the  desired  position.  Thus 
both  ax  and  Ay  are  zero.  The  usefulness  of  this  mathematical  relation  in 
planning  is  to  answer  the  following  question: 

"If  the  measurements  are  in  error^how  serious  is  the  resulting 

position  error?" 

This  question  can  be  answered  (demonstrated  here  for  sextant  case)  by  entering 
a  known  systematic  error  into  the  L  matrix;  for  example,  all  elements  of  L  are 
1  minute.  The  summation  notation  Tor  this  is  provided  by  G  and  H  of  equaTions 
(8-3).  The  resulting  vector  is  called  the  systematic  error  tendency  (SET)  and 
is  defined  by: 


AXset  = 
•^set  = 


HE-CG 

AE-C2 

AG-HC 

AE-C2 


(8-8) 


The  obvious  use  of  this  is  to  examine  all  possible  fix  geometries  to  determine 
which  ones  are  least  susceptible  to  inaccuracy. 

The  component  of  V  in  a  predesignated  direction  is  discussed  in  appendix  D  as 
a  possible  measure  for  use  in  setting  standards. 

The  derivation  of  selection  and  prediction  equations  now  turns  to  the  stochas¬ 
tic  nature  of  the  observations.  The  random  measurement  error  associated  with 
each  line  of  position  and  the  arbitrary  reference  variance  have  been  assumed 
for  planning  purposes.  These  assumptions  allow  determination  of  the  confi¬ 
dence  in  the  determined  position.  The  confidence  is  described  in  terms  of  a 
two-dimensional  probability  distribution  centered  on  the  desired  position  (no 
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systematic  error  involved).  Various  descripters  of  this  distribution  are 
defined  as  position  error  measures.  The  derivation  is  as  follows. 

The  matrix  of  assumed  measurement  variances  is: 


(B-9) 


By  a  Jacobian  transformation  (reference  13),  the  matrix  of  coordinate 
variances  is 


(AT  W  A)"1  = 


For  two  measurements  (for  clarity), 


is  reduced  to 


( B— 10 ) 


all2wll+a212w22  alla12wll+a21a22w22 
a12al lwll+a22a21w22  a122wl l+a222w22 


<7 


2 

0 

(B-ll) 


Normally,  a  correlation  exists  between  coordinates  and  C^O.  To  remove  this 
correlation  (reference  13)  the  axes  are  rotated  clockwise  ai  angle  /?ce 
according  to 
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ft CE  =  V2  tan-1  (B-12) 

With  the  following  conditions,  /3q£  represents  the  orientation  of  the  major 
semi-axis  w.r.t.  north: 

I<  0  (45° 

=  0  then  /3  c£  =  |  0° 

>  0  (  135° 


if  A <  E  then  the  minor  semi-axis  has  been  found  w.r.t.  north 
and  90°  must  be  added  to  find  the  orientation  of  the  major  semi-axis. 
Finally, 


if  /J  CE  <  0  add  180° 

The  new  uncorrelated  (u,v)  coordinate  system  allows  the  covariance  matrix  to 
be  defined  as 


with 


(B-13) 


A+E-((A-E)2+4C2)  1/2 
A+E+((A-E)2+4C2)  1/2 


(B- 14) 


A]  and  B^  are  called  the  geometry  factors. 

The  confidence  ellipse  is  defined  on  th^  xi  v. uncorrelated  coordinates 
(reference  12)  by 


/u  -  *u\  2  +  /v  -  ^  v\  2 
\  ^maj  /  \  a min  / 


(B- 15) 


where  «  is  the  confidence  level  desired  and  and  ? v  are  the  Cartesian 
coordinates  of  the  AT.  In  the  planning  case,  and  fiv,  are  zero  unless 
systematic  error  is  simulated. 
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The  expected  position  error  is  given  in  terms  of  the  bivariate  normal  proba¬ 
bility  distribution  defined  above.  The  major  and  minor  semi-axes  of  a  desired 
confidence  ellipse  are  determined  by  a  .  The  semi-axes  are  multiples  of 
°  maj  and  <rmin  where  the  multiple,  M  „  ,  is  determined  by 


( B - 16 ) 


(B-17) 


Table  B-3  provides  a  selected  set  of  multiples  and  the  corresponding  confi¬ 
dence  that  it  can  be  expected  that  the  true  position  will  lie  within  the 
specified  ellipse  centered  on  the  AT  (reference  17). 


TABLE  B-3 


MULTIPLE 


V4.605  =  2.15 
V5T99T  =  2.45 
V372T0  =  3.03 


PROBABILITY  LEVEL 

90% 

95% 

99% 


The  area  of  the  one  standard  deviation  ellipse  is  determined  as  follows  (det 
represents  taking  the  determinant  of  the  matrix): 


Area(T=  *  (det  _£u) 1/2  =  *  (det  (  °  g  p-  1(AT  W  A)-]P)) 1/2  (B- 18) 

where  P  rotates  ^  u  to  uncorrelate  coordinates.  Continuing, 

AreaCT=  *  (  <r  J  I  P'^A1  W  A)  “  ^1 ) 1/2  =  I  (AT  W  A)'1!  I  P  I  1/2 

(B-20) 

2  2 
T  a  0  _  T  °  Q 

I  At  W  A  I  1/2  ( AE-C2) 1/2 

The  area  of  other  confidence  ellipses  are  determined  by: 


Area, 


2 

0 


(AE-C2) 


T72 


(B-21 ) 
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where  is  determined  at  some  confidence  level,  a  . 

Ii  some  situations  the  confidence  ellipse  semi-diameter  in  a  specified  direc¬ 
tion  is  important. 

Three  previously  defined  quantities  are  needed  to  calculate  the  semi -diameter 
of  a  particular  confidence  ellipse  in  a  direction,  s  ,  clockwise  from  geodetic 
2  2 

north.  They  are  ^  ma  j .  a  m-jn>  and  the  ellipse  orientation  angle,  /?  ce  •  The 
major  semi -diameter  oriented  in  a  direction  $  is 


( B - 22 ) 


Another  position  error  measure  useful  in  planning  is  the  probability  mass  (of 
the  p.d.f.  centered  about  the  AT)  contained  within  a  designated  circular 
region  of  radius  R  (centered  on  the  AP),  called  P-in-R.  Calculation  of  this 
measure  is  outlined  in  appendix  G  of  reference  2.  P-in-R  is  a  function  com¬ 
bining  both  the  accuracy  and  precision  in  positioning.  Further  discussion  of 
P-in-R  is  contained  in  appendix  D  herein  which  discusses  standards. 

The  final  expected  position  error  measure  is  the  2-drms  value.  It  is  found  by 


(B-23) 


A  circle  with  a  radius  of  2-drms  centered  on  the  center  of  the  ellipse 
encloses  at  least  95%  of  the  probability  mass. 
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APPENDIX  C 


MATHEMATICS  OF  POSITIONING 


C.l  COORDINATE  SYSTEM  FOR  ERROR  ANALYSIS  AND  NOTATION 

The  origin  of  the  coordinate  system  is  the  assumed  position.  The+x 
direction  is  to  the  east  and  the+y  direction  is  to  the  north.  Each  fix  geome¬ 
try  of  n  measurements  consists  of  n: 

a.  Lines  of  position  with  their  gradient  magnitude  and  direction, 

G-j  and  'f  and  relative  weights,  w^-j. 

b.  Precomputed  measurements,  *c-j. 

c.  Observed  measurements  (corrected  for  systematic  errors),  a0-j. 

C. 2  DETERMINATION  OF  THE  ANGLE  TAKERS  POSITION  (AT) 

The  equations  to  calculate  the  AP-to-AT  vector  in  planning  are  used 
iteratively  in  calculating  the  AP-to-AT  vector  when  positioning.  Each  itera¬ 
tion  differs  only  by  the  choice  of  the  assumed  position.  The  desired  position 
is  used  for  the  first  iteration  and  successively  determined  AT's  are  used  as 
the  AP  for  following  iterations.  This  process  accounts  for  error  due  to  the 
linearization  of  the  system.  The  final  values  of  ax  and  Ay  are  the  sums  of 
all  incremented  changes  in  x  and  y.  The  sums  are  used  to  convert  back  to 
latitude  and  longitude  with  equations  (C-4). 

From  appendix  B  the  components  of  the  AP-to-AT  vector  are: 


BE -CD 

Ax  =  — £  (summed  over  all  iterations) 

Ml  —  v 


(C-l) 


AD-BC 


2  (summed  over  all  iterations) 


The  magnitude  of  the  AP-to-AT  vector,  V,  is: 


(  ax)2  +  (  A y)2 


The  direction,  /fyf,  of  V,  w.r.t.  geodetic  north,  is: 


(C-2) 


ax  >  0 


Ay>0  tan-1 


AX  <  0 


360°  +  tan-1 

•a  Jr 


Ay  <  0  180  o+  tan-1  —  180°  +  tan*1  -4^ 

Ay  Ay 


0°  <  flAT  <  3600  ( c-3) 


C-l 


The  conversion  of  the  ax  and  Ay  (which  are  in  meters)  to  ax  and  a<£ 
is  accomplished  through  (referencel7) : 

A<£  *  (Diff.  per.  sec.)  Ay 

(converts  a  y  in  meters  to  a£  in  seconds) 

ax  *  H  ax  (C-4) 

(converts  a  x  in  meters  to  ax  in  seconds) 


where  from  USC4GS  S.P#241  (other  references  exist  that  provide  these  equations 
in  slightly  different  forms), 

Diff.  per.  sec.  =  1/(111132.09  -  566 .05cos2^p  +  1 .2cos4<*>Ap)  (C-5) 

H  =  -Diff.  per.  sec./cos4>Ap  (C-6) 

The  geographic  latitude  and  longitude  of  the  AT  are  found  by: 
aAT  *  XAP  +  AA 

.  .  CC-7) 

‘PAT  s  ^AP  +  A* 


If  the  angle  takers  are  displaced  from  the  chain  stopper,  the  vectors  V  and 
Dj  (appendix  C.3)  should  be  added  before  finding  geographic  coordinates. 

C.3  SYSTEMATIC  ERRORS 

C.3.1  Lack  of  Observer  Coincidence 


It  is  desirable  for  observers  to  stand  at  the  same  point  when 
making  measurements.  Of  course,  this  is  difficult  due  to  inter-visibility 
conditions  and  therefore  a  systematic  error  of  observer  lack-of -coincidence 
exists  and  should  be  compensated  for  before  computations  are  performed  on  the 
measurement  set.  Reference  2  explains  the  significance  of  the  lack  of 
coincidence.  It  was  found  to  be  small  in  most  cases  but  the  equations  to 
compensate  for  the  error  are  presented  here  for  completeness. 

Define  the  displacement  vector,  0-j,  to  each  observer  from  some 
reference  point  in  the  region  where  they  make  measurements.  The  corrected 
observation  is  aoi,  which  represents  the  measurement  that  would  have  been 
made  if  the  observer  were  standing  at  the  reference  point.  The  corrected 
observation  is  found  as  follows.^  The  gradient  vector  corresponding  to  each 
LOP  is  G-j .  The  vectors  G -j  and  are  used  to  find  <»0^  from  the 
uncorrected  measurement,  a  j. 


®oi  *  n  i 


I  Oi  I 

“i  +  .■==”  cos(/?h  +  .doi  -  *i) 

I  U-j  I 
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where  the  direction  of  0^  is  determined  by  the  sum  of  the  heading  direction, 
and  the  direction  0-j  makes  the  clockwise  from  which  is  called 

8  oi  • 


C.3.2  Angle  Taker  to  Sinker  Drop  Point  Vector 

Observers  are  not  able  to  stand  near  the  sinker  drop  point  during 
positioning  operations.  The  displacement  vector,  Tfs ,  from  the  observation 
point  to  the  sinker  drop  point  may  be  compensated  for  after  preliminary  deter¬ 
mination  of  the  position  of  the  angle  takers  (AT).  The  vector  from  the 
desired  position,  AP,  to  the  calculated  position  of  the  observers  (AT)  is~V. 
Let  ,3  h  represent  the  true  heading  of  the  positioning  platform  (clockwise 
from  north)  and  let  ps  represent  the  direction  of  Ds  clockwise  from 
/3h.  The  AP-to-MPP  vector,  7C,  is  computed  by 


V  +  D< 


(C-8) 


_  The  direction  of  0S  is  the  sum  of  8 ^  and  ps.  The  direc¬ 
tion  of  V  was  discussed  in  appendix  C.2. 

The  angle  taker  to  sinker  drop  point  vector  can  also  be 
compensated  for  by  placing  a  scaled  model  of  the  positioning  unit  on  the 
gradient  diagram. 

C.4  DETERMINATION  OF  PRECISION  IN  POSITIONING 

The  residuals  of  a  fix  can  be  used  to  determine  the  precision  associated 
with  the  determined  position.  A  residual  is  the  difference  between  an 
observed  measurement  and  an  adjusted  measurement.  The  least  squares  procedure 
minimizes  the  sum  of  the  squares  of  weighted  residuals.  In  matrix  notation, 
the  residuals  are  (using  the  matrix  notation  of  appendix  B). 


R  =  AX  -  L 


(C-9) 


The  sum  of  the  weighted  squared  residuals  is: 


RTVJR  =  (AX  -  L)T  W(AX  -  L) 

=  xWwAX  -  XTATWL  -  LtWAX  +  JWL 

When  R~HwR  is  minimized  (reference  13,31): 

X  =  (AtWA)'1AtWL 

XtAtWAX  =  XTAtWL  (C-10) 

RTWR  =  LTWL  -  L^WAX 


A 


3 


1 


An  unbiased,  most  likely  a  posteriori  estimate,  s?,  of  the  reference 
2 

variance,  a  ,  is  provided  by  (references  13,31): 

0 

RtWR  LtWL  -  JWAX 

c2  =  —  —  =  -  — _ — — r-  ir 


where  n-2  is  the  number  of  degrees  of  freedom  of  the  estimate.  Two  degrees  of 
freedom  were  lost  to  determining  the  most  probable  position. 

The  estimator,  s^,  can  be  reduced  to  summation  notation  (see  appendix 
B)  as  follows  (for  2  measurement  case;  higher  numbers  follow  readily). 


2 

wn  0  1 

Ami 

'wn 

0  ‘ 

"4 11 

412" 

AX 

(n-2)s  =  [Am]  Am2] 

0  W22 

a  m2 

-  [Am]  a  m2] 

0 

w22_ 

.421 

422. 

Ay 

(n-2)  s^  =  Ami^wii  +  Am2^W22  -  Am]W]]a]]  ax  -  Am2W22a21  Ax  “ 

Am-|W]iai2Ay  -  ^m2w22a22  ( C-  12 ) 

and  finally 


s2  =  (F  -  B  ax  -  0  Ay)/n-2  (C- 13) 

(units  of  reference  variance) 

The  variate  (n-2)  s^/  <r  *  is  x  n  ,,  distributed  and  is  used  to  calculate 
the  confidence  ellipse  dimensions. 

From  equation  ( B - 1 5 ) ,  the  confidence  ellipse  on  the  axis  of  uncorrelated 
coordinates  is  defined  as. 


u  -MU" 

2 

r 

> 

a. 

1 

> 

1 _ 

"maj 

+ 

‘'min 

2,« 

2  2  2 

The  reference  variance,  *  is  related  to  <t  „  .  and  .  by  equations 

u  maj  min 

(B-14).  Using  equations  (B-13)  in  (B-14), 

2 


( C- 14) 


|u  -  M  u 

2 

-f 

V  -  My 

r 

Q 

O 

1 _ 

B]  <^0 

2,« 


(C- 15) 


By  dividing  both  sides  by  (n-2)s2/<rjj 
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J 


,1 


u  -  M  2  /v_-j*v\  2  2 

A1  0  /  V  B1  0  /  _  ‘2,  a 

(n-2)  s2  (n-2)  s2  2 

2  - 2  *  n-2, a 

"  0  17  0 


an  F -distributed  variate  arises  using  General ized-T^  statistics  (references 
12  and  15). 


2  F2,n-2,« 


( C- 16) 


/“  -  2  +  h  -  ^v\  2 

\  Sm3j  /  \  smin  J 

where  Ms  =  y  2F2tn-2,a- 

Ms  is  found  for  different  values  of  n  and  <*.  Table  C-l  provides  a  set 
of  selected  Ms  values  for  various  n  and  a  (reference  17).  The  major  and 
minor  semi-axes  of  the  confidence  ellipse  are: 


(C-  17) 


As  =  Ms  smaj  -  Ms  A]  s 
Bs  =  Ms  Smaj  =  Ms  s 


TABLE  C-l 

CONFIDENCE  ELLIPSE  MULTIPLIERS 


n 


Confidence 

2 

3 

4 

90% 

00 

9.94 

4.24 

3.30 

95% 

00 

19.97 

6.16 

4.37 

99% 

30 

100.00 

14.07 

7.85 

Equation  (C-16)  demonstrates  that  only  the  relative  weights  of  the 
measurements  are  important  to  confidence  ellipse  determination.  The  relative 
weights  are  used  in  calculation  of  A-| ,  B],  and  s.  If  only  one  measurement 
type  is  used,  all  weights  are  equal  and  the  reference  variance  is  that  of  the 
measurement  type  used. 

The  orientation  of  the  major  semi-axis  w.r.t.  north  is  provided  by 
equation  (8-12)  as  it  was  for  planning. 
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The  semi -diameter,  As§  ,  of  the  confidence  ellipse  of  confidence  level 
in  some  direction,  $  ,  clockwise  from  geodetic  north,  can  be  determined  by: 


AsS  = 


Ms  s 


iin2(  S  -/?c£)  cos 2 (  8  -/?CE)  \ 


'1  U1 
The  area  of  a  confidence  ellipse  of  confidence  level  a  is. 


(C- 19) 


?  *  s2 

Areas  =  AiSis2  =  - —TTT  (C-20) 

S  (AE-C2)1/2 

where  Ms  is  determined  by  a  and  the  number  of  measurements. 

A  measure  of  position  error  which  combines  the  major  and  minor  semi-axis 
into  one  number  is  the  circle  of  confidence,  CO C.  The  circle  of  confidence 
encloses  at  least  1  -<*  of  the  probability  mass  and  is  found  by 


COC  =  Ms  s  ^  A2  +  B2  (C-21) 

The  circle  of  confidence  is  centered  on  the  most  probable  position. 

C.5  REPLICATION  -  TIME  AVERAGING  LOPS 

Suppose  that  the  ability  to  replicate  each  of  the  n  lines  of  position  in 
a  fix  geometry  m  times  exists  (e.g.,  LORAN  receivers  do  this).  The  effect  of 
time  averaging  "n  times  m"  measurements  on  fix  precision  is  discussed  in  this 
appendix . 

All  position  error  measures  of  precision  are  functions  of  the  a 
posteriori  estimate  of  the  reference  variance,  the  confidence  level  multiple, 
Ms  (appendix  C)  and  the  geometry  factors.  A]  and  B].  If  each  of  the  n 
lines  is  measured  with  equal  precision,  the  geometry  factors  are  not  dependent 
on  the  number  of  replications.  The  main  influence  on  precision  is  through 
Ms  and  s2. 

For  n  lines  of  position: 


s 


2  = 


(C-22) 


Let  each  residual  consist  of  a  systematic  component,  rs,  and  a  random 
component,  rr.  The  residual  matrix  is  now  written  as: 

=  R_s  +  Rj. 
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The  estimate  is: 


S2  = 


(R$  +  R r)T  W  {R3  +  Rj,) 

rT2 


(C-23) 


Replicate  each  line  of  position  m  times  randomly  and  average  the 
residuals.  The  random  components  will  average  close  to  zero  if  m  is  large 
(>100)  and  only  the  systematic  components  remain.  Thus, 


<2  3  «sT  M  h 
n-2 


(C-24) 


If  no  systematic  error  exists,  the 
variance  would  be  near  zero. 

The  confidence  level  multiplier  Ms 


estimate  of  the  reference  measurement 
is  given  by  equation  (C-16), 


Ms  * 


(C-25) 


As  the  number  of  measurements  gets  large  (with  “=.!),  F£  n_ ?  ^ 

approaches  2.3  and  Ms  approaches  2.15.  This  is(not  coincidentally)  the 

same  as  x  |  1  • 

The  geometry  factor  A]  is  not  dependent  on  the  number  of  replications 
made  of  each  line  of  position  as  it  is  computed  using  the  orientations  and 
gradients  of  each  of  the  n  lines  of  position. 

The  effect  of  m  replications  is  displayed  by  the  equation  (C-13)  for  the 
major  semi -axis. 


VrJ  w  Rs 

„.2-  (C-26) 

The  result  of  this  derivation  indicates  that  little  is  to  be  gained  by 
time  averaging  of  lines  of  position  if  systematic  error  exists.  There  is  an 
upper  limit  on  the  precision  with  which  a  position  can  be  determined.  The 
results  of  this  section  dispel  the  common  belief  by  many  mariners  that  time 
replication  always  provides  a  very  accurate  fix  with  high  precision. 
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APPENDIX  D 


STANDARDS-POSITION  ERROR  MEASURES 


The  purpose  of  this  appendix  is  to  explore  various  alternatives  for  a  position 
error  measure,  S,  when  considering  standards  for  aid  positioning.  A  standard 
is  defined  as  an  authoritative  measure,  for  comparison  with  S  to  deter¬ 
mine  if  S  is  acceptable,  which  indicates  that  the  positioning  evolution  has 
been  acceptable.  Sn,  may  depend  on  the  specific  environment  at  the  time  of 
positioning,  the  availability  of  signals,  the  quality  of  the  signals,  the 
importance  or  criticality  of  the  aid,  and  the  difficulty  in  achieving 
(i.e.,  limit  on  number  of  attempts  at  reaching  at  the  time  of  position¬ 
ing.  More  than  one  Sm  value  may  be  applicable  to  any  one  positioning  evolu¬ 
tion.  The  specific  combination  depends  on  the  requirements  placed  upon  the 
aid  at  the  time.  Successful  aid  positioning  is  determined  by  both  accuracy 
and  precision;  this  means  the  applicable  set  of  Sm  values  must  at  least 
contain  tests  for  these  two  qualities. 

Frequency  histograms  of  position  error  measures  found  through  extensive 
research  of  historical  data  received  from  the  REDWOOD  (WLM-685)  are  pro¬ 
vided  at  the  end  of  this  appendix.  The  data  for  each  position  error  measure 
are  separated  only  by  chart  scale.  All  fixes  are  3-line  fixes  using  Third 
Order  or  better  reference  objects. 

D.l  A  PORTER IOR I  ESTIMATE  OF  REFERENCE  VARIANCE 

The  reference  variance  estimate,  s^,  of  equation  (C-ll)  can  be  used  as 
a  measure  of  the  precision  of  a  fix  (reference  13).  The  estimate  is  unbiased 

and  can  be  tested  against  some  a  priori  value,  *q,  using  the  x2  test.  With 
n  measurements. 


x 


2 

n-2 


-  (n-2)  s? 
2 

"0 


(D-l) 


has  a  x  distribution  with  n-2  degrees  of  freedom.  There  are  two  possible 
hypothesis  tests.  First, 


HO 


versus 


Hi 


and  reject  Hq  when 


(0-2) 


2  2 

x  <  x  or 

n-2  n-2>l-  «/2 


2  2 

x  >  x 
n-2  n-2,  <*/2 


(0-3) 
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x 


Second, 


o  2 

Hq  :  <r c  -  a  g  versus  H] 
2  2 

and  reject  H0  when  x  „  „  >  x  ~  0 

u  n-2  n-2,« 

0PM. 

The  standard  is  selected  at  some  a 

~  ^2  >  w^ere  ^m  = 

*0 


:  0  2  >  <7  g  (D-4) 

The  second  is  more  applicable  to  the 

such  that 


x 


2 

n-2,  “ 


(0-5) 


2 

indicates  an  adequate  positioning  evolution.  <?  ^  need  not  be  constant  or  even 

representative  of  measurement  random  error.  It  can  represent  some  critical 
value  based  on  other  considerations  beyond  which  rejection  occurs.  It  is  most 

2 

advisable,  however,  to  establish  a  g  through  analysis  of  historical  data. 

Figure  0-1  is  a  cumulative  frequency  histogram  of  s  for  fixes  performed 
by  the  crew  of  the  REDWOOD  (WLM-685).  The  data  has  been  separated  into 
two  groups  by  chart  scale.  The  figures  indicate  that  positioning  efforts  on 
aids  found  on  charts  of  scale  1:20,000  or  less  result  in  s  values  of  on  the 
average  7.2  minutes.  Positioning  of  aids  located  on  larger  scale  charts  leads 
to  an  average  s  of  4.4  minutes.  Remember  s  is  an  unbiased  estimates  of  <rg. 


D.2  A  POSTERIORI  ESTIMATE  OF  LOP  VARIANCE 

The  distance  residuals  are  defined  as  the  distances  in  meters  between  the 
AT  and  the  LOPs  along  a  line  perpendicul ar  to  the  LOPs.  In  terms  of  the 
measurement  residuals,  the  matrix  of  distance  residuals  is  given  by 


D=GR=G(AX-L) 


( D-6 ) 


where  G  is  an  nxn  diagonal  matrix  of  LOP  gradient  magnitudes. 

The  distance  residuals  are  weighted  by  the  usual  weighting  matrix  defined 
in  equation  (8-1),  squared  and  summed. 

The  sum  is  differentiated  and  set  equal  to  zero  which  leads  to 


x  =  (aVwga)-Wwgl 


(0-7) 


0-2 


This  result  is  identical  with  that  presented  in  appendix  B  with  W 
redefined  as  G^WG.  In  this  case,  the  reference  variance  is  in  units  of  the 
LOP  variance  Tmeter2).  One  way  to  calculate  the  LOP  reference  variance  is 


a  Z  .  1 
lopo  n 


n 


Gi 


G2  <t  2 


(D-3) 


2 

where  the  <*  .  are  the  assumed  measurement  variances. 

! 

The  calculation  of  the  reference  variance  estimate  can  be  altered  to 
provide  an  LOP  variance  estimate.  The  estimate  is  found  as  follows: 


DTWD  =  LtGtWGL  -  LtGtWGAX 


An  unbiased  most  likely  a  posteriori  estimate  of  the  LOP  variance, 
s^opt  is  provided  by 


2  DtWD  LtGtWGL  -  LtGtWGAX 

4P  ■  =s?  ■  — vr""  <D-9> 

where  n-2  is  the  number  of  degrees  of  freedom  in  the  estimate.  The  possible 
hypothesis  are  similar  to  those  in  the  previous  section. 

The  confidence  level,  a  ,  can  be  selected  to  account  for  the  conditions 
at  the  location  of  the  aid.  The  standard  is  selected  at  some  <*  such  that 


(n-2)  s 


lop 


lopo 


Sm  where  Sm  « 


2 

x 

n-2,  a 


(0-10) 


indicates  an  acceptable  positioning  evolution. 

Figure  D-2  is  a  cumulative  frequency  histogram  of  s op  for  fixes  per¬ 
formed  by  the  crew  of  the  REDWOOD.  For  small-scale  charts,  s -j op  averages 
6.8  meters  and  for  large-scale  charts,  siop  averages  9.9  meters.  Remember, 
Slop  an  unbiased  estimate  of  <*  lopg* 

D.3  CONFIDENCE  ELLIPSE  PARAMETERS 

The  dimensions  and  orientation  of  the  confidence  ellipse  provide  four 
position  error  measures,  S,  at  any  confidence  level,  <*  .  They  are: 

a.  The  major  semi-axis. 

b.  The  area  of  the  confidence  ellipse. 

c.  Semi-diameter  in  specified  direction  S  . 

d.  The  square  root  of  the  sum  of  the  squares  of  the  major  and  minor 
semi-axes.  Table  D-l  directs  the  reader  to  the  equations  needed 
to  calculate  the  position  error  measures. 
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a* 


TABLE  D-l 


r 

! 

{ 


POSITION  ERROR  MEASURES  FOR  CONFIDENCE  ELLIPSE 


MEASURE  EQUATION 


C- 18 

Ass 

C- 19 

Areas 

C-20 

COC 

C-2 1 

Figures  D-3,  0-4,  and  0-5  are  cumulative  frequency  histograms  of  confi¬ 
dence  ellipse  parameters  for  CGC  REDWOOD  fixes.  Average  results  are: 

SCALE  As ( 0 . 90 )  Areas(0.90)  COC  (0.90) 

<  20  39  meters  4,200  m2  43.3  meters 

>  20  66  meters  12,400  m2  75.8  meters 

0.4  AP-to-MPP 

The  magnitude  and  direction  of  the  AP-to-MPP  vector,  Vc,  is  calculated 
using  equations  (C-l),  (C-2),  (C-3),  and  (C-8).  Setting  a  standard  on  the 
magnitude  of  Vc  is  equivalent  to  establishing  the  level  _of  acceptable 
accuracy;  acceptable  precision  is  not  established  using  Tc.  In  many  situa¬ 
tions  the  component  of  Tc  in  a  specific  direction  may  be  important.  For 
example,  in  marking  a  narrow  channel,  the  component  along  the  channel  is  of, 
little  importance  relative  to  the  transverse  component.  The  component  of  Vc 
in  any  direction  <P ,  clockwise  from  geodetic  north,  is  found  by 

m  ffc  *  t)t  =lYcl  cos  (/?Mpp  -  *  )t  (0-11) 

A 

where  t  is  a  unit  vector  at  the  AP  in  the  <p  direction.  The  standard, 
is  established  such  that. 


1  Vc1  <  5m 

or 

nrC|i  <  Sr* 


indicates  an  acceptable  positioning  evolution. 

Figure  D-6  is  a  cumulative  frequency  histogram  of  the  magnitude  of  the 
AP-to-AT  vector  for  REDWOOD  fixes.  The  average  vector  magnitude  for  the 
small-scale  chart  case  is  13.1  meters.  The  large-scale  case  gives  a  magnitude 
of  23.4  meters. 


(0-12) 

(0-13) 


0.5  P-in-R 


The  probability  mass  (of  the  p.d.f.  centered  about  the  MPP)  contained 
within  a  designated  circular  region  of  radius  R  (centered  on  the  AP)  is  called 
P-in-R.  Calculation  of  this  success  measure  is  outlined  in  appendix  G  of 
reference  2.  P-in-R  is  a  function  combining  both  the  precision  and  the 
accuracy  in  positioning.  Calculation  of  P-in-R  requires  two-dimensional 
numerical  integration  (which  is  very  time  consuming);  numerical  error  in 
calculating  P-in-R  is  largely  dependent  upon: 

a.  The  number  of  points  at  which  the  integrand  is  evaluated. 

b.  The  limits  of  the  integration. 

c.  The  quadrature  method  employed. 

d.  The  higher  order  derivatives  of  the  integrand. 

Table  D-2  displays  the  relative  error  in  calculating  P-in-R  using  various 
orders  or  Gaussian  Quadrature  (reference  18)  for  an  assortment  of  R/amaj 
and  "maj/ "min  ratios.  The  relative  errors  are  presented  as  the  percent 
difference  between  the  calculated  values  and  the  true  P-in-R  for  each  case. 

The  time  required  to  calculate  P-in-R  is  a  function  of  the  order  of  integra¬ 
tion.  From  the  table. 


Time  (secs)  =  5  (#  pts)3/2  (on  HP-41C)  (D-14) 

where  (#  pts)  is  the  order  of  Gaussian  Quadrature. 

The  regions  of  the  table  where  the  relative  error  is  less  than  10%  are 
boxed  in  to  indicate  the  applicability  of  P-in-R  as  a  measure  of  success  in 
positioning.  The  calculations  indicate  that: 

a.  P-in-R  cannot  be  calculated  accurately  with  low-order  ( <  8) 
quadrature  methods  unless  the 


"maj 

ffmin 


<  3  and 


R 

"  maj 


<  3 


(D-15) 


b.  If  P-in-R  is  to  be  calculated  for  accurate  inmediate  feedback  when 
positioning  (<1  min),  then. 


Isfil  +  _R  , 

"min  "maj 


5 


( D  - 1 6 ) 


c.  Accurate  feedback  in  less  than  15  seconds  is  restricted  to 


+  JL_  <  3 

amin  "maj 


(0-17) 


D-5 


-  — '-'A- 


TABLE  D-2 


PERCENT  ERRORS  IN  CALCULATING  P-IN-R  AND  COMPUTATION  TIME  (HP-4 1C) 


"  maj 

R 

ORDER  OF  INTEGRATION  (#PTS//TIME  IN  MINUTES) 

"min 

"min 

4770767m 

6// 1 ,25m 

8//2.00m 

16//5.3m 

1 

!  o 

0 

0 

0 

0 

2 

21 

0 

0 

0 

0 

1 

3 

61 

\l _ 

0 

0 

0 

4 

90 

28 

3 

0 

0 

5 

99 

59 

14 

1  2 

0 

6 

100 

82 

32 

1  8 

0 

1 

4 

0 

0 

0 

o, 

2 

62 

11 

- Lo. 

0 

0 

3 

96 

54 

15 

1  3 

0 

4 

100 

87 

48 

19  j 

1 0 

5 

100 

92 

77 

44 

ll 

6 

100 

100 

93 

68  | 

Li 

i 

25 

1 

0 

0 

0 

2 

94 

47 

13 

9 

0 

3 

100 

91 

59 

28 

1 0 

4 

100 

100 

91 

65 

Ul 

5 

100 

100 

99 

89 

17 

6 

100 

100 

100 

97 

36 

■>'rr*>rniirr 

■1 

1 

55 

[_9 

0 

0 

0 

2 

100 

76 

38 

14 

0 

3 

100 

98 

87 

60 

lil 

4 

100 

100 

92 

21 

5 

100 

100 

100 

99 

48 

6 

100 

100 

100 

62 

1 

80 

23 

U _ 

0 

0 

2 

100 

90 

63 

33 

f  01 

3 

100 

100 

97 

84  | 

m 

4 

100 

100 

100 

99 

46 

5 

100 

100 

100 

100 

75 

6 

100 

100 

100 
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91 

1 

93 

40 

10 

2 

0 

2 

100 

95 

80 

53 

m 

3 

100 

100 

99 

100 

31 

4 

100 

100 

100 

100 

69 

5 

100 

100 

100 

100 

91 

6 

100 

100 

100 

100 

98 

Time  *  5  ( #pt s ) 3/2seconds 
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d.  P-in-R  can  be  calculated  accurately  if  time  is  not  an  important 
restriction  (>5m)  for. 


gmaj 
a  min 


R 

^maj 


<  3 


(0-18) 


In  calculation  of  table  3-5,  the  AP-to-AT  vector  magnitude  was  assumed  to 
be  zero.  Vc_does  affect  the  relative  error  but  does  so  in  a  positive  way. 

In  fact,  if  Vc  exceeds  R,  the  accuracy  of  the  integration  is  significantly 
improved.  The  reason  for  this  being  that  the  region  of  integration  is  where 
the  p.d.f.  changes  slowly  with  position. 

Figures  0-7,  0-8,  and  0-9  are  cumulative  frequency  histograms  of  the 
P-in-R  values  calculated  for  REDWOOD  fixes.  The  average  P-in-Rs  are  as 
follows: 


RADIUS  IN  METERS 

SCALE 

R=Tu 

R=T5 

<20 

0.47 

0.68 

0.80 

>20 

0.10 

0.19 

0.34 

0.6  R-for-P 

The  radius  of  a  circular  region  (centered  on  the  AP)  that  contains  at 
least  the  probability  mass  P  (of  the  p.d.f.  centered  on  the  MPP)  is  called 
R-for-P. 


R-for-P  is  defined  as  the  sum, 
R-for-P  =  As  +  IV 


(0-20) 


which  is  just  the  sum  of  the  major  semi -ax is  of  the  confidence  ellipse  and  the 
magnitude  of  the  AP-to-MPP  vector.  The  standard,  Sm,  is  established  at  some 
level  a  such  that 


R-for-P  <  Sm 


(0-21 ) 


indicates  an  acceptable  positioning  evolution  or  classifies  the  position 
determined . 

Although  this  approximation  is  rather  crude  at  some  R,  the  P  is  in  error 
by  no  more  than  10  percent  and  any  numerical  methods  for  calculating 
R-for-P  more  accurately  are  extremely  time  consuming. 

Figure  0-10  is  a  cumulative  frequency  histogram  of  R-for-P  (for  »  =  o.l) 
from  data  taken  from  the  REDWOOD.  For  small-scale  charts,  the  average  R-for-P 
is  51.6  meters.  Fixes  associated  with  aids  on  large-scale  charts  average  89.2 
meters . 
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0.7  DIFFERENCE  BETWEEN  0BSERVE0  MEASUREMENTS  AND  COMPUTED  MEASUREMENTS 


0.7.1  All  Measurements  -  Statistical  Test 


In  section  4.2,  procedures  for  precomputing  measurements  (expected 
to  be  observed  at  the  assigned  position)  were  discussed.  After  numerous 
operations  at  a  given  station,  the  positioning  team  should  develop  high  confi¬ 
dence  in  the  expected  measurements.  In  cases  where  high  confidence  exists, 
measurements  can  and  should  be  compared  to  those  that  are  expected.  The 
expected  measurements  have  no  estimated  parameters,  which  causes  the  statisti¬ 
cal  comparison  to  be  one  level  less  complicated  than  the  consistency  checks 
described  in  the  previous  sections. 

Assuming  a  priori  measurement  variances  (which  should  also  be 
known  quite  well  after  frequent  visits  to  a  station),  the  ratios  Am-j/wi 
are  squared  and  summed.  This  test  is  much  more  sensitive  than  the  test  on  the 
residuals  in  that  it  takes  into  account  more  than  measurement  inconsistency; 
in  addition,  it  tests  measurement  accuracy. 

The  test  is  a  x  2  test  as  follows: 


swd  = 


"  Ami2  . 
f-1  <r  2 
i 


* 


2 

n 


where  swd  is  pronounced  "the  sum  of  the  squared,  weighted  differences." 

The  positioning  evolution  is  acceptable  if  the  calculated  swd  does 
not  exceed  the  standard,  S^,  which  is  found  in  x  2  tables  at  various 
confidence  levels,  a  . 

In  the  same  way  s^op  is  related  to  s2,  the  sum  of  the 
squares  of  the  gradient  weighted  differences,  gwd,  is  related  to  swd.  The 
formula  for  gwd  is, 


gwd 


n  G?  Am? 

V  -  1 

1-1  G2  *  2 


2  2 
x  where  Sm  =  x 

n  n ,  a 


Figures  D-ll  and  D-12  provide  cumulative  frequency  histograms  for  Am  and 
Gam  from  all  REDWOOD  fixes  studied.  The  average  values  are: 


SCALE 

G  a  m 

a  m 

<20 

13.8  meters 

16.3  minutes 

>20 

20.1  meters 

13.9  minutes 
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0.7.2  Special  Case  for  Difference  Between  Observed  Measurements  and 

Computed  Measurements;  Example  Use  of  Standards 

One  of  the  most  frequently  employed  procedures  used  by  the  Coast 
Guard  in  positioning  aids  is  called  the  fixed  glass  procedure.  The  fixed 
glass  procedure  is  performed  by  setting  two  sextants  at  prescribed  angles  and 
maneuvering  the  ship  until  both  measurements  agree  with  the  prescribed  angles. 
Presently,  this  procedure  often  includes  only  two  measurements  and  should  be 
extended  to  satisfy  the  requirement  for  at  least  three  measurements.  To  do 
this,  the  positioning  team  continually  makes  a  third  measurement,  thus  allow¬ 
ing  a  continuous  check  on  the  consistency  of  the  measurement  set.  The  problem 
is  to  define  the  limits  within  which  the  third  measurement  must  lie  so  that 
the  measurement  set,  as  a  whole,  indicates  an  adequate  positioning  effort. 

The  derivation  of  the  relationship  that  exists  between  the  third  measurement 
difference,  am3,  and  resulting  position  errors  is  as  follows. 

The  fixed  glass  method  required  that  two  observed  measurements 
agree  with  the  corresponding  computed  measurements.  Let  the  third  observed 
measurement  differ  from  the  corresponding  computed  measurement  by  some  value 
a  m3  and  still  indicate  an  adequate  positioning  effort.  That  is,  some 
measure  of  success  S  is  a  function  of  Am3  and  the  standard  Sm,  corre¬ 
sponds  to  a  limit  on  a  m3.  As  an  example,  the  major  semi -axis,  As,  is 
found  as  a  function  of  Ann,  start  with  equations  ( C- 18 )  and  ( C- 13)  with 
(n=3) . 

As  =  MSA1S 

,  LtWL  -  lTwax  t  t 

S2  =  r- -=r-r  =  iJWL  -  LTWAX 

Because  two  measurements  are  on  the  mark,  only  one  a  m-j  remains.  Which  one 
depends  on  which  of  the  three  is  the  third  measurement,  s'f  depends  on 
A  m-j  as  follows: 

s<t  =  Am-jw-jj  Amn-  -  AtriiWii  [aii  ax  +  a^  Ay] 

Using  the  summation  notation  of  equation  (B-3),  ax  and  a y  are: 


BE-CD 

AE-C2 

AD-BC 

AE-C^ 


si  =  AmjWjj 


-  AmiWii 


ajitBE-CD) 

AE-C2 


a ,2( AD-BCi 


AD-A107  all 
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For  this  case,  B  and  D  can  be  reduced  to 


and. 


B  =  Airmail  0  *  a ,2 


A 


Arn^wii 


A  m-jw-j  ^ 
AE-C2 


jja-j !  ( a -j i  A  m-jE-Ca-j2  An,i )  +  ai2(A  Amia^-Caii  Ami  )1 


But, 


sin  ri 

- 51 — 


cos  Y  i 


ai2 - 


and  finally. 


2  2 
Si  a  A  m -w  -j  i 


The  major  semi -axis  is 


/Acos2  y  i  +  Esin2  r-j-  2Csin  Yicos  yA* 
\  G2.  (AE-C2)  / 


As  *  MsAi  Ami(wii)1/2(i_ki)l/2 

where 


Acos2^,  +  Esin2  y i  -  2Csin  Yicos  Yi 
G2  (AE-C2) 


The  minor  semi-axis  is 

Bs  *  MSB-|  Ami(w1-j)V2(i-ki) V2 
and  the  confidence  ellipse  area  is 


2  2 

Areas  =  M^Bi  Am^w-j i ( 1-k ^ )  = 


Anfojjd-kj) 
(AE-C2) 1/2 
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The  approximating  circle  of  confidence,  COC,  is 


COC 


=  Ms  ^  +  8^  AmiCwii)1/2  (!-ki)V2 


and  finally  the  semi-diameter  of  the  confidence  ellipse  in  some  direction  S, 
clockwise  from  geodetic  north  is  given  by 


As 


Ms  Ami(wii)1/2  (1-kj)172 
4in2{  S  -#ce)  .  cos2(  8  -Pqi) 


'i 


B 


1 


where  Pqe  is  the  orientation  of  the  confidence  ellipse  clockwise  from 
geodetic  north. 

The  results  of  this  derivation  shows  the  linear  relationship  that, 
exists  between  the  measurement  difference,  a  m3,  and  the  major  semi -axis. 

The  slope  of  the  linear  function  is  dependent  on  the  confidence  level  desired, 
the  weight  of  the  third  measurement  relative  to  the  two  other  measurements  and 
the  geometry  of  the  fix. 

In  the  following  example,  the  major  semi-axis  is  calculated  as  a 
function  of  the  third  measurement  difference  for  a  geometry  which  consists  of 
three  sextant  measurements  of  equal  weight  and  which  determine  LOPs  with 
y  i=(0,60,120)  G-j=l  .0.  The  90%  confidence  level  ellipse  major 

semi-axes  is  calculated  producing  the  following  results: 


As  =  MsAt  Ami(wii)l/2(i-ki)l/2 
Wii  =  W22  =  W33  =  1.0 
A  =  E  =  1.5  k]  =  k  2  =  k  3  =  0.67 

C  =  0  A]  =  0.67 

Ms  =  V2F  2,1,0.10  =  9-94 


Using  Ms  in  the  As  expression,  we  have 


As  =  (9. 94)(0. 67)(l-0. 67)V2  Ami  =  (3.33  5*1)  Ami 


The  geometry  is  symmetric  about  the  designated  crossing  point 
which  means  the  major  semi-axis  is  not  dependent  on  which  of  the  three 
measurements  is  chosen  to  be  the  check  angle.  The  maximum  acceptable  Am^ 
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which  still  indicates  an  adequate  fix  is  determined  by  the  chosen  semi -major 
axis  standard  and  the  proceding  equation. 


Am, 


S 

m 


(3.83 


met. 

mTrr 


For  example,  if  a  major  semi -axis  standard  of  twenty  meters  is 
established  for  the  example  geometry  studied,  the  first  two  measurements  are 
caused  to  coincide  by  maneuvering  the  ship  and  the  consistency  is  adequate  if, 


a  m3 


< 


20  meters 


3.84 


met 

min 


*  5.2  minutes 


In  other  words,  the  angular  measurement  must  be  within  5.2  minutes  of  the 
computed  measurement.  Similar  calculations  can  be  performed  using  other 
standards;  such  as,  confidence  ellipse  area,  directional  semi -diameters. 


NOTE:  Once  the  third  measurement  differences  have  been  calculated  using  one 
of  the  above  equations,  they  can  be  used  anywhere  on  the  grid  diagram 
to  check  measurement  consistency.  This  procedure  is  as  follows: 

On  the  grid  diagram,  use  dividers  to  find  the  distance  in  minutes 
from  the  point  where  two  lines  of  position  intersect  to  the  third 
line  of  position  (along  a  perpendicular  line).  The  result  must  be 
less  than  the  standard  set  for  the  third  measurement  difference. 
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APPENDIX  E 


DETECTION  OF  MEASUREMENT  OUTLIERS 


E . 1  A  POSTERIORI  REFERENCE  VARIANCE  ESTIMATE 

The  use  of  s2  in  detecting  outliers  is  similar  to  its  use  as  a  position 
error  measure.  In  appendix  0.1,  the  statistical  tests  are  outlined  for  com¬ 
paring  s2  to  expectations.  Rejection  of  the  hypothesis  that  <j  2=  <r  2 
should  lead  to  serious  review  of  the  data  for  blunders.  ^ 

E.2  MEASUREMENT  COMBINATIONS 

In  any  case  where  more  than  three  measurements  are  made  in  positioning,  a 
statistical  analysis  of  various  combinations  of  the  measurements  can  be  per¬ 
formed  with  the  goal  of  detecting  a  measurement  which  is  an  outlier  at  some 
level  of  confidence.  This  method  was  alluded  to  by  Rosenblatt  (reference  12) 
and  outlined  in  reference  2,  appendix  F. 

The  analysis  begins  with  a  set  of  n  inconsistent  measurements  at  a  pre¬ 
designated  confidence  level  «  .  This  is  indicated  by  the  test 


1  n-2Ls.2 
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The  problem  is  to  identify  which  of  the  n  measurements  are  causing  the  incon¬ 
sistency.  The  inconsistency  may  not  be  due  to  only  one  poor  measurement;  how¬ 
ever,  it  is  logical  to  test  for  this  possibility  first.  The  approach  is  to 
temporarily  remove  one  measurement  from  the  set  of  n  measurements  (Q)  times 
thereby  creating  n  different  reduced  fixes.  The  reduced  a  posteriori  estimate 
of  the  reference  variance  from  each  of  the  (")  reduced  fixes  has  n-3  degrees 
of  freedom.  One  degree  of  freedom  lies  with  the  measurement  residual  of  the 
removed  measurement.  This  residual  is  calculated  by  comparing  the  actual 
measurement  with  that  measurement  which  would  have  been  made  if  the  observer 
were  at  the  reduced  MPP.  There  are  now  (0)  reduced  estimates  of  the  reference 
variance  and  an  equal  number  of  removed  residuals.  The  removed  residual  is 
weighted  by  the  corresponding  element  of  the  weighting  matrix,  W,  which 
allows  the  following  ratio  to  be  calculated  (Ij1)  times. 


(E-2) 


The  numerator  has  one  degree  of  freedom,  the  denominator  has  n-3  degrees 
of  freedom.  The  largest  of  the  (9)  ratios  is  tested  against  the  F  statistic 
at  some  level  <*  .  Thus 


>  Fl,n-3,a  (E-3) 

max 

The  hypothesis  (that  the  maximum  ratio  is  not  excessive)  is  rejected  if  the 
calculated  ratio  exceeds  the  critical  F]  n_3  a  value.  In  any  case  where 
the  hypothesis  is  rejected,  the  corresponding  measurement  (that  was  tempor¬ 
arily  removed  for  testing)  can  be  rejected.  The  remaining  n-1  measurements 
can  be  further  tested  for  compliance  with  any  applicable  standards,  such  as 
against  the  newx 2  value  (n)  now  being  reduced  by  one  from  that  in 
equation  (E-l). 

As  mentioned  earlier,  there  is  no  guarantee  that  only  one  measurement  of 
the  n  measurements  is  unacceptable.  In  fact  all  measurements  can  cumulatively 
contribute  to  an  unacceptable  inconsistency  in  the  measurements.  The  proba¬ 
bility  that  all  measurements  are  in  large  error  is  small;  but  one  can  easily 
envision  scenarios  where  two  measurements  contribute  to  the  inconsistency. 

For  example,  suppose  that  a  rebuilt  but  not  resurveyed  landmark  is  sighted  in 
two  separate  sextant  angular  measurements.  Though  the  measurements  are 
statistically  independent,  there  is  a  high  probability  that  they  are  mutually 
inaccurate.  The  problem  is  to  statistically  identify  two  inaccurate  measure¬ 
ments  and  (assuming  they  cannot  be  corrected)  identify  them  from  the  measure¬ 
ment  set.  It  is  obvious  at  this  point  that  the  number  of  measurements 
required  to  perform  statistical  outlier  detection  is  quite  large.  If  one  were 
to  attempt  to  identify  two  statistical  outliers  from  a  set  of  only  four 
measurements,  only  two-measurement  subsets  remain  and  the  tested  measurements 
are  of  the  same  strength  as  the  reference  measurements;  this  is  an  obvious 
predicament.  Therefore,  it  is  justified  to  eliminate  the  four  measurement 
cases  from  the  attempts  to  identify  two  inaccurate  measurements. 

Where  five  or  more  measurements  are  made  (be  they  any  type),  statistical 
detection  of  paired  inaccurate  measurements  is  possible  but  the  method  becomes 
quite  cumbersome.  Often  graphical  representation  of  the  lines  of  positions 
and  the  corresponding  landmarks  would  be  more  fruitful  as  performed  in  refer¬ 
ence  28.  The  measurements  are  now  grouped  into  (j))  different  (n-2)-measure- 
ment  subsets  with  the  corresponding  (£)  removed  two-measurement  groups.  The 
(n-2) -measurement  subsets  can  be  used  to  calculate  a  posteriori  reference 
variance  estimate,  s^,  with  n-4  degrees  of  freedom.  The  removed  two-measure¬ 
ment  groups  can  be  used  to  calculate  two  residuals  which  are  each  weighted  by 
the  corresponding  element  of  the  weighting  matrix,  W.  This  allows  the 
following  ratio  to  be  calculated  (£)  times. 
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The  numerator  has  two  degrees  of  freedom,  the  denominator  has  n-4  degrees 

of  freedom.  The  maximum  of  the  (?)  ratios  is  tested  against  the  F  statistic 
at  some  level  a  .  Thus 


>  F2,n-4, a  (E-5) 

max 

The  hypothesis  (that  the  maximum  ratio  is  not  excessive)  is  rejected  if  the 
calculated  ratio  exceeds  the  critical  n-4  value.  In  any  case  where 
the  hypothesis  is  rejected,  the  corresponding  two-measurement  subset,  that  was 
temporarily  removed  for  testing,  can  be  rejected.  The  remaining  n-2  measure¬ 
ments  can  be  further  tested  for  compliance  with  any  applicable  standard,  such 
as  against  the  new  x  £_2  value  (n  now  being  reduced  by  two  from  that  in 
equation  (E-l)).  Simultaneous  inaccuracy  in  more  than  one  measurement  can  be 
partially  avoided  if  no  landmarks  are  used  for  more  than  one  measurement. 


Measurement  combinations  increase  rapidly  after  the  case  (^)  so  further 

statistical  tests  are  superfluous.  In  fact  (?)  leads  to  ten  different  mea¬ 
surement  checks  and  might  be  superfluous  itself.  It  is  mentioned  here  only 
because  the  common  landmark  misplacement  case  can  contribute  sign  if icantly  to 
position  error  (reference  2).  The  search  for  two  outliers  should  be  performec 
subsequent  to  the  search  for  one  outlier. 

E.3  DIFFERENCE  BETWEEN  OBSERVED  MEASUREMENTS  AND  COMPUTED  MEASUREMENTS 

In  section  6.2.7,  swd  was  discussed  as  a  position  error  measure.  Swd  can 
also  be  used  to  detect  large  measurement  errors.  A  test  for  outliers  can  be 
made  on  any  subset  of  m  measurements  of  the  n  measurement  set.  This  test  is 
as  follows. 
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at  some  confidence  level  a  for  each  m  measurement  subset.  The  number  m  would 

logically  progress  from  one  to  two  with  the  tests  on  many  combinations  (")  and 

(~),  of  measurements.  Given  that  any  test  indicates  an  outlier,  the  outlier 
can  be  rejected,  and  the  remaining  measurements  can  be  used  to  compute  the 
position. 


APPENDIX  F 


HP-41C  AND  OFFSHORE  CALCULATOR -ASSISTED  RESECTION 


F.l  INTRODUCTION 

Calculator  routines  were  prepared  to  demonstrate  analytical  positioning 
using  the  HP-41C  Programmable  Calculator  System.  The  system  includes: 

HP-41CV  Programmable  Calculator 

HP-82153  Optical  Wand 

HP-82143A  Thermal  Printer 

The  report  sections  listed  below  were  chosen  for  use  in  the  demonstration 

3. 3. 1.3  4.2.1  5.1.1  6.2.3  8.2  C.l  D.l 

3.3.2  4.2.2  5.2  6.2.4  C.2  D.2 

4.2.4  C. 3.2  D.3 

C.4  D.4 

Three  sections  follow  in  this  appendix: 

1.  Section  F.2  is  a  listing  of  the  positioning  program  for  the 
HP-41C  System 

2.  Section  F.3  is  a  brief  flow  chart  to  instruct  the  user  how  to 
operate  the  routine.  The  program  and  operating  instructions 
are  provided  in  a  form  usable  by  those  already  familiar  with 
operation  of  the  calculator. 

3.  Section  F.4  is  a  bar-code  generation  routine  developed  for  use 
on  the  HP9825A  Calculator  and  the  HP-9872B  Plotter  (adapted 
from  Generating  Barcode  in  the  Hewlett-Packard  Format.  McNeal , 
Thomas,  BYTE,  January  1981).  The  routine  can  be  used  to 
generate  optical  bar-codes  for  the  HP-41C  programs  and  for  the 
geodetic  control  data  for  landmarks  and  buoys.  The  routine 
requires  the  23K-byte  option  on  the  HP-9825A  and  requires  the 
following  Read-Only-Memories  to  be  inserted: 

Advanced  Programming  9872B  Plotter 

String  Variables  General  Input/Output 

Extended  Input/Output 

Once  again,  the  routine  has  been  provided  in  a  form  usable  by  those 
already  familiar  with  the  calculator  and  its  operation. 


F-l 


The  bar-code  generation  routine  is  useful  in  two  different  modes  of 
operation,  they  are  as  follows: 

Program  Mode  -  In  this  mode,  the  routine  is  capable  of  creating 
complete  pages  of  bar-code  that  represent  calculator  programs  that 
can  be  loaded  into  memory  by  use  of  the  wand.  Once  the  data,  pro¬ 
grams,  and  special  function  keys  have  been  loaded  into  the  HP-9825A 
Calculator,  the  routine  is  operated  by  a  list  of  commands  that  are 
assigned  to  the  special  function  keys.  The  command  and  input  format 
are  given  by  File  14  of  this  appendix. 

Geodetic  Oata  Mode  -  In  this  mode,  the  calculator  uses  only  sub¬ 
routine  DATA  to  generate  the  optical  bar-code.  The  bar  generated 
represents  the  landmark  number,  the  latitude  of  the  landmark,  and 
the  longitude  of  the  landmark.  The  routine  prompts  the  operator  for 
the  required  input. 


F-2 


F.2  HP  -  41C  SUBROUTINE  LISTINGS 
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F.4  HP  -  9825/  HP  -  9872  BAR-CODE  GENERATION  ROUTINE 


FILE 

LIST 

FUNCTION 

NAME 

File 

0 

INITIALIZATION  ROUTINE 

BAR-CODE 

File 

1 

DATA  REQUIRED  TO  RUN  PROGRAM 

DATA  FILE 

File 

2 

PROGAM  BAR-CODE  PLOTTING  ROUTINE 

PROG 

File 

3 

SPECIAL  FUNCTION  KEY  STORAGE 

SPEC.FUNCT 

File 

4 

CONTROL  OR  PROMPTING  ROUTINE 

PROMPT 

File 

5 

AUTOMATIC  NUMBERING  ROUTINE 

NUMBER 

File 

6 

PROGRAM  LISTING  ROUTINE 

LIST 

File 

7 

GEODETIC  DATA  BAR-CODE  PLOTTING  ROUTINE 

DATA 

File 

8 

COMPILATION  ROUTINE 

COMPILE 

File 

9 

PROGRAM  RENUMBERING  ROUTINE 

RENUMBER 

File 

10 

empty  file 

File 

11 

empty  file 

File 

12 

PROGRAM  STORAGE  ROUTINE 

SAVEPROG 

File 

13 

PROGRAM  RETRIEVAL  ROUTINE 

GETPROG 

File 

14 

COMMAND  DEFINITION  ROUTINE 

COMMANDS 

File 

15 

BAR-CODE  PLOTTING  ROUTINE  FOR  USER 

DEFINED 

BIT  PATTERNS 

BARPLOT 

SPECIAL  FUNCTION 

KEY  ASSIGNMENTS 

fO 

/YES 

f8  /NUMBER 

fl6 

/SECPROG 

fl 

*ldf4 ,0 ,0 

f9  /COMPILE 

fl7 

*ldpO ,0 

f  2 

/EXIT 

flO  /DATA 

f  18 

blank 

f3 

/SAVEPROG 

fll  DELETE 

fl9 

blank 

f4 

/PROG 

fl2  blank 

f  20 

/RENUMBER 

f  5 

/SCRATCH 

fl3  blank 

f  21 

/ABORT 

f  6 

/NO 

fl4  blank 

f  22 

/BARPLOT 

f7 

/LIST 

fl5  /GETPROG 
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KEYS 


‘-'.r.Mm*,*.. 


file  o 

Os  cv-:..j»  ; 

::  as.-,  A*:i5003.F*trfC'.::.«*c'5o.:].nc?5o.i3 
3  s  31"  ft  !2  ]«sni  3.B:  ;?a  3.Tst  ?0  TO  ].  ■■■■sc  301.  SSt  50  ].&*[  1 
j:  :  J. :«  IS-S.7MC  ;i4].C*Cie.l  :.Cti03 

-s:  ;at  3 

f:  -  :  -  is  '30Sfi;  •  -l  -fUC  I  ;-*»[  ;  ];n*  t  ; 

E:  '  "AS 
7:  0-T 

i  i  8  ♦  C  t  *  :* 

9:  '* 

u*:  .2*  :.us«  i»*  cc*3 

11 J  inr  ’’Cdnnand  Iltt  in  n  ioni‘‘  *  yf 

12:  :*  :  ,*  ‘N"  or  ca&< U* •*  NO’ ; lof  4.0 

!!•:  '2t  14.0 


FILE  2 

VS  ’Sv0f  *• 

::  wr.  ?03«-VS35“ 

::  ...r-  755.  ' IP  100*  100.8000. 10004“ 

2?  i  *  f  1  9S  *  j'to  2 

4.:  or*  “Cowell*  Pro*"«ldf  4*0 

f:  -?•••  "P-D?ri.-.  No Mf*  Llnii  30  we^t^i)'‘»T$ 

*£i  iCl  0*8.0*  10 
7:  :ii:  1*2*1.23 

i:  218-C  s  b**o »  wai  ?  4005  b*e»*  wo  1 t  400;t»®*p 
8:  *  •*»?  3b  *  “P  ip*  r  in  Plot*.  «r“»  1  lb 

•  ■  •'*  0  *  C  *  I » 2  «  C  •  C 1 C  *  C  •  C  *  3  3  •  C  •  I  •  I*3*C»C»C*C*C*C5stP 

1 : :  o*n#  1 

:::  - 3. 5*?. "5*1 

:■!  ibi  Tf 
14:  •  ■'  R  ' ?  > -» ; 1 

If:  :*  ftnodT'  0*  X+l+X 
l’  cit  3*9.5.  1 

l":  1  r  I  "Nunc*''  ot  R*9ii*frs  N*»d»d 

1  :  1 0  1  ! ! 


33!  ino  6 
33:  3*S" : 

34:  itfuncc*!  :>"scs*i : 

33-  3  -  1 -3 : 3"i-3»E-l"£ 

3 £ :  ;•  E#0ij.i»  3 

3’:  0-f 

33:  »'o  "Be" 

33:  E.QSjmP  3 

I :  ?i  o  “  Be  ~ 

33:  ■.!  :cr:«:g*i*e 

;•  i?  * '  t'UC  C  ]->  •  143!  ?t  o  -26“ 


:*  !•  *  ■  MIC  C  ]  '*25  JM  lU 'MS[C:>#30!  jnp  4 
:•  *  3  MS [ C" 1  33fioalo"l4£ 
r*!  -»f 

»i  9  "8;" 

it  i*r ■  ntic v<i6  cr  Kf in*tc:> -as; jnp  9 

C  "  1  4  I 

it  it  *  'MSC  I  ]'<  16  or  itf(HSCn>,28;ji.B  3 


42: 

i  *  i  "i 

4  n 

•>  o  -2 

44: 

I  -C- 1  -»E !  r  -*-1  "F 

4;  : 

:•  E*Q"?!.o  "8c" 

44,* 

C-F 

47: 

•to  "Sc" 

4-3: 

0-E-F 

4?: 

.to  "Be" 

30! 

11 38  “ : 

c  *  • 

:*  i  •  f  f  MFC  C  l‘<  >29': 

no 

"3B“ 

w  *  • 

;•  ltf  <M«CJwi92 

or 

u  ■  11SCC3 

e  -  . 

'♦!-F 

e  ±  . 

i*  C"2''ftSji?p  3 

re  . 

2-r 

e_- . 

•  to  "8c" 

r  — . 

;  •  •  ■  f’ft  02  JOi’toellE- 

2+c 

e  :  . 

•  to  '3c" 

f  ?: 

'--£ 

4.j: 

f 

4  *  : 

?!  0  "  8c  ■ 

4  2 : 

38 

4  ‘ : 

;•  i'  t  ■  ft*tc:.  .  240 : 

j  «*i© 

4 

c4 : 

3-E 

44 : 

••o  "6c" 

•*  “ : 

1  *  t 1  rif  C  C  3  -toiHE-S 

82: 

•*1*F 

c  r : 

'6c " : 

i  v  • 

it  3  IE  o nd  C  rt*  •* 

7  •* 

J. 
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FILE  2(cont) 

:  *  -h; e-** 

"2:  -i  3  3 

F'i  1  *5 f « ♦  G •*» o si  1 6*SC  2  1 
"•*:  *0 r  :*i  »c  D 
"T :  M*5C  I  3^0'3256-m 
"■?:  ;♦  M  *2S*iM«od25E*l*M 
”>*  x  I 
: : 

'  2o rw«  rf  ■  J 
cO 5  G»:-G 
;  1 :  *  *  ®  1  =  1  -o  D 

s:=  sc  ;»• 

£2:  tor  Y*2**I  to  3 ♦  r- *  I  —  1  b*.  -l 

S4;  pNrtodS^Bt  >'  3 
IT:  -EC"  l*  2-Y 

it*:  ntrct  • 

5":  rtfM  I 

■:•*:  *-BC  1  ]-eC2>BC8D-4  2 
:  ?■:  1  *»BC  3D*3  2 

40  •  jf  C 
4  *  !  3P^4*  J 

-2:  oit  $• 4, S2-. 5G«od!5» 1 
42:  *  <d  0 

4-:  :  c  i  -*P0LI  "  s  1  b  1  C*  !bl  Lilbl  “  -  “?lbl  Bflbl  "  •" 

4T:  c  i :.  ■?«  4- .  5GnocU 8’  1 
?c*:  tor  1*1  *  o  J 

-"*•  :♦  BCl>0SioIt  0..3>2Si»:t  0* - . 3» -1 J  .tnP  2 

4  :• :  ic-lt  O  • »  2  *  2  S  i  p  1 1  0*-.3»-l«iPlt  .  0 1 8  *  0  ♦  i  p  1 1  0  > .  3  •  2 »  : » 1 1  0  •  • .  3 »  - 1 
licit  .04.0.1 
1O0:  ni/*  I 

*01*  !♦  Gr»Odli*0«  CUP  'Crsan?*  Poo®  r  1  ‘  "  •  it u 

102:  ina  S ■ B 

103*  3-D5  B-»L 

104:  it  E*0iL*l*L 

*•!“:  .t  Z  (4;  sit  o  "  ITS  '1 

1 0 1 :  o»*t  'BarCode  Conclftc" 

10":  idt  4,0 

■*  o  .3 1  *0 

FILE  4 

6>  s  "Prompt": 

1  I  W  •  T* 

2:  po;  ■.  Tt  i  ”  "  >  ■>  I 
3 i  O-f’-rl  I  ct  ,  1 
4:  1-N 

5:  it  1*0!  ,to  "Command” 

T*;i,!-s  :-u* 

7:  it  US*"DEL£T£'‘  !  jno  4 
i :  TIC !*1, 1+4 3*U* 

9:  Ur.  .'Jf»l-I 
i Os  It  5  1 

it  I-l:4*crt  "#  *00  lar,*"i,to  "Prompt" 

13:  to-  J=  I  - 1  to  1  b  -1 

if  uiIJOJ'T  or  >J*C  J.  “>n  *  or  Cmd"i,to  "Prompt  " 

i  4 :  ♦  .  riun UICJ»J3^*46)]0TK•♦'v, 

IT : 

1  J 

;7s  i*  V  750!p-t  '*  *oc  lar«"!,,.o  "Prompt" 

1  : :  i *  *  1 ,1.1  *  jmo  3 
13:  tt;  '  - 1  >»P*C  V  3 
30:  no  "Prompt" 

3::  ■:;:*iwi 
33:  41iT« "  '  "-ft* 

33:  ft;  7-'-P«V: 

34 :  T*liOtT$'*l*T 
25:  v  o  " P romPt  " 

4o :  Command  : 

27:  if  T**"NO”!.to  "cronot" 

23:  it  T**"BhPPLOT"S ldf  15.0 
33:  it  Tf*”DflTrt" 1 ldf  7.0 
30:  if  T, ."NUMBER"! ldf  5-0 
31!  .•  tf* "RENUMBER"! Id*  9.0 
33:  it  Tf » " L I 3T " ! 1 df  o.O 
33:  :•  Tf.'PPOG’i ldf  3*0 
3-s  Ti« "COMPILE"!  ldf  5*0 
35:  i •  "I*"YE5"ildf  14,0 
33:  it  TJ«"3flVEPP0G"! ldf  12*0 
3“:  it  TS*"GETPROG" ! ldt  13,0 
33i  TS.'SCPflTCH" 1 jne  5 

39:  «r.t  "Do  r.all.  want  to  BCRmIGHV  ,  js 

40:  if  oacf  TS>«"NO"  or  cap<T* '»"N'  !  ,to  "Pronpf 
4;:  1*1  •;  750ifti  ’  -1 : -nt[  l  3-Plt  l  3!  net  I 

43:  0-rE»rs-»7i  "••#!*!  ,to  "PronBt  " 

43:  it  T**"E:-.:T" ! a*B  'Don«"!»tB 

44:  I*  TS.'SECPROG" !  »rt  "''■'Unr.c  ond'.fto  "Pronot” 

45:  3-N 

49:  Id,  3 • 0 

-31403 
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FILE  5 

IN  **..n  ctr'i 
::  ent  " :  ■  a  r»  "  .  V 

2  J  ’  Inc .  "  *  r 9 

2:  it  iOdd*c-t  #  *  :*o  i  ere* "Sid*  4*0 

41  »  .  .3  0 

t:  d£*  * 

c :  en»  “  “  ♦  T  S 

•:  i  -  t="t. \  -Jf  4,0 

t  :  Ml\ ’’SU  "-Mi 

?:  1 1  »  vr  -Pf[  ] 

i*>!  7*!*rv  Tf  •♦1-T 

115 

12?  :  ric  -  9 

•2419 


FILE  6 

0:  "L 1  it " : 

is  * wi  diets 
•  tor  I»l  to  750 

it  :  1 1  ■  f  it  I  ]  ><Oi  ,t  o  ♦€ 

•or  J*1  to  50 

:•  **C  1 1 1  ,P«  I  jt,  .p*(  i  2  ,«.j  j,»  i  ••;  JrtB  2 

ne  :t  J 

:  9i£  :•  *  <  P*C  IJ>*t  i  it  f  'p*c  I  ]  “J-t  l+Si 

ft:  ijrt 
?:  ne.-t  I 

:0s  ut  4.0 

“  1  2?*DCi 


LE  7 


2  ■  o  :  ;  n  o  &.Ki04rl0-rii-ri2+rl3»""+TJ 
i.ir-  '05.  "VS23" 

300*0 

*np  .  rl0 

err  "Lti  N(i,..*o(;50 

?->r  "Lett  if..jde0<DD.  MMJSist  >“•  r 1 1 ;  erridf  rll .  -6>*rl  1 
*-s  "Lonai '  ijd*°<DE .  MM3  Jsss  /  “  .  r  12.  ornd'  r  1 2t  -6  ;*r  1 2 
*r>t  "Loca-ion  of  Bar-Cod*"'"  t  rl4S  if  rl4V16; jnc  0 
ri4  J-l 

..irt  -05.  "IP  000 .7000 .7900. 5000" ii  to  "Plot"  2,5: 

:  ■  i  rt  705" " rP0C0i 6000. 7900. 7000" : *t  0  "Plot"  o-. 
u-t  -05.  "TP000I  5P00. '000.  €000"  i  ,t  o  "Plot"  j"--. 
urr  -05.  IP000. 4000.  7900. 5000  Soto  -Plot"  oo.l 
..:rt  -05.  *  IP000. 3000. 7900. 4000  "«*to  "P'.Of  on. 
w-  ”05. " IP0O0. 200C . 7900*  3000"  i  if  c  “Plot"  ot. 
or-  705.  "IP000.  1000.  -900. 2000-:  ro  "Plot  "  «, 

705.  "I  p000. 000.  7900.  1000  ■•Silt  O  "Plot"  j7. 
urt.  705.  "IP5100I  7000,  13000.3000"!  910  "Plot"  ~/t\ 
•»*■•  -05.  ’  IF 5 100. 4000.  1 3000.  7000" 1 it O  "Plot" 

«-•  705.  "IP5100, 5000. 13000.  €000"!  *  to  "Plot"  Yd\ 

-  -  *  -05. "  I P5 100. 4000.  1  30O0 1  5000" '■  it  o  Plot"  37. 
u-t  -05.  "  IF5100.  3000  ■  1 3000*  4C00" :  it  o  "Plot"  . 

1  r.  -05-  '  1P51C0. 3000.  13000,  3000"!  it o  “Plot"  So, 
..irt  -05.  IP5100.  1000.  13000.2000"!  ito  "P  lot"^. 
705.  '  IP5100. 000- 1  3000*  1000'  '  it  0  "Plot"  «•; 


J^ILE  7  (cont ) 


,'t  -05" 
705.' 
*05.' 
-.-e  -05.  ' 

1  r.  -05. 

« rt  -05. 
•■r-  '05. 
•Plot": 

155 


FILE  8 


Sl+S 

"  -S* 

»v'-  J*1  to  730 

i*  i  m  'PSC  L  3K0*  **  o  "next  J” 

for  R*1  to  30 

if  itf  <P«t  iPSC  !"1  2 

n«  xt  P 

**C  it*  <P*C  d3)*t*  it  f  -,PSC  J 1  +  T$ 

S-r9 

icin’*. 

if  ABORT" « ldf  4,0 
“--U**v* 

0*ri +r2 
-i-v 

jop  3 

*•;  0'.5,CS1 

:  +  l-5 

?  licti  2!cti  3!ct  »  4!cti  5!o.i  6 

•S*$« 

i*  TI*"£ND"  or  Tl»" .  END.  “ !  it  o  'End" 
i-  -fCltll,”"  ard  T*t  l ,  2  3, "ft'  ' !  too  14 

;■  Tft "! j.-ip  3 

7*C2  1-T*:ort  T $ 

1 1 1  i 

l«rt  T*  i  »L 
It  L  13!  je,p  3 


4o:  oen#  1 
27:  sol  0> ”t0. 1 
25:  -s:x  10.it. 125 
29:  “Lab*:"; 

30:  2 1 5-C ,  b**p!  wa 1 1  400;  oe* p!  ...a  i •  4 
3n  *n*  ib.  "facer  in  Plotter  .lit 
32:  yrt  0.C.C.C,C.C,-;.C,C.33tC.C.C. 
33:  nit  , 3>.7,i;ibl  "Ln  " : t . d  0!1l 
2a:  Pi*  . 3. . 55. 1 ! lbl  "Lat  "!*rd  7; \ 
35:  Pit  2.. 55.1! lbl  "Lon  “ : lbl  '12 
3€:  "  Number”: 

37:  in-  rl0/-100)-Kt2: 

33:  1 00f  rc f  r 1 0/ 1 00 ; -KC  2  1 
39:  "Lat Lon"; 

4p:  <nr  I»u  to  12 

41:  int '  rl  4»4f  l-l  l  ■  3 

42:  int  1 100f  ro<  r I )  ;*t  £  5*4< l-l 1 y 3 

43:  ir.t  (  l00f  rc  <  1 06r I  -K£  6*4 .  ; - 1 ;  ; 

44:  int  <  100f  ro'  10000rl  >  >»5£  7»4t  I-n 

45:  next  I 

4©:  "Chkium". 

47:  0-rl3 
43:  for  1*2  to  11 
49:  rl3*KCn»rl3 
50:  next  I 
5 1 :  r 1 dnod25S-KC 1 1 
52s  "Concert": 

53:  for  1*1  to  11 
54:  KC I  lex 

55:  tor  7*2+31  to  3+3'.  1-1.’  Ot  -1 

56:  Xmod2  +  6C  Y 3 

57:  .X-BCV1J  -2+X 

53:  next  7 

59:  next  I 

60!  0+BC  1  1-BC21+BI  92  3 
61 :  1-BC91  1 
62:  "Pit": 

63:  Pi •  . 3* . 2. 1 
64:  tor  !*1  t.o  92 

65:  it  3£I3*0!i»lt  0..3"2!ie!f  0.-. 
66:  ip  It  C't .  3>  2!  ipH  0.  -.  3!  iPit  .01 
€“:  ip  1*  ,03.0.1 
6: :  re  t  I 

69:  ol*  0.0. ; ! IPit  4.5>0.2!iPlt  0> 1 
70:  c*n#  O 

7;:  enp  More  LM'l""'.-i 

72:  it  C0P,T»'*  'N"  or  o«p.  T,  •■'NO'  ! 

7":  *•  o  "Eat  a" 


. .  9  .  t  *  L  .  .  • 
10! It! 
r  1  1 


.  3 :  i  p  1  •  0  -  - . 


■•It  -4. 5.0;  ic  It  .j.-: 
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FILE  8  (cont ) 


cr t  "l_in«  too  lonfi  »".J 
jnp  3 

If  TIC  Lt  L  3« . ijnB  4 

c-t  "Rlpna  Error  *  *J 
fib  "Error" 

«■!  "Scon" 
iff  2 

fto  “IS" 

tor  to  lenCTi' 

■ft  M  3-U* 

it  and  U*'»“9“ijnp  4 

If  or  ll*»"-"'  and  l«n>T*  < ;  1 1  one  3 

l*  US*'  "  or  Uf*"E”  or  Uf«"."<jne  2 

MP  4 

ni.rt  I 

t.  f  4  £ 

ftO  "IS" 

pos  <  T*  •  “  '">  +  r  1 

it  rllOfjnB  3 

iff  3 

fto  "IS" 

tic  r  1*1  :-ui 

T*C  It  rl-t  3-T* 

ien’.lJf't-L 

if  UfCI.n#"’"  or  UfCL.U#"’ “ijrtp  5 

if  L  '9 *  ft  o  -25 

UfC  2>  L-l 3-US 

iff  4 

ft o  "IS" 

post  Ijf  < "  "  1*r2 

if  r2=0>.mp  8 

JiC  1 • r2-l 3-V* 

if  VI*" I  HD " ! Jnp  4 

prt  "Operand  Error  #"tj 

pso  "Error" 

ft o  "Scan" 

itf  5 

'Jft  r2»l  3-Ui 
if  I tn I U* t < »2! ft o  “IS" 
or t  "Nun.  Op.  Error  #"tj 
fib  "Error" 
fto  "Scan" 

"IS": 

It  f ; s2»l i an b  26 

ft:  <240-l-23»m«S3 

if  MflMWti  <imt1*CS3'-l>*MICS3 

S-l-5 

1-.Y 

*?-Y 

flflllijnp  3 
*  >  i  ■'tSD-MlISJ 
3-1-S 

* o **  :*2  to  l-i 
'-EJ  v-C 

i r.-  '  •.B-C>'2)-M 
:  f  "'ll  I  -  !  IfCttMJ!  ;nn  5 
;•  TfC  I  •  I  3/CIC II  3>  ft- 1-8 
i<  TIC  I  >  I  3'  CIC  M  3!  M-l  -C 

It  6-*Cijnn  -4 
0-11 
n-; 


sr ; 

if  2i0t  Jnp  4 

&h: 

prt  "Char.  Error  »" 

t  J 

Sf*: 

fib  'Error" 

90: 

ft.o  "Scan" 

9 1 : 

"i  CC23'*n*tS3 

92.: 

3-1-S 

91: 

n*ft  I 

9*: 

ft.o  next  J" 

ci*: 

;f  flf€»l<fto  "Other" 

q*;: 

'Dif "Scf f  9i if  Tl*" 

ABORT "i Idf 

97: 

if  TIC  1 » 1  !•“♦"  and 

TIC  1 . 1 31"- 

?Z: 

if  Tit  1. 1 3I---J Jnp 

3 

QC; 

•»l  <  23  > -MIC  S3 

100: 

S-l-S 

1 0 1 

TfC  2  3-TI 

102: 

sou  ft i "  " >-r 3 

1*3: 

noi' Tf  t "E" )-r4 

104: 

l*  r3«0> l*n(T#)-r3 

107: 

it  r4#0I r*-l*r3 

1 06: 

for  I»1  to  r3 

10": 

If  TIC  1 1  I  "  1  jnp 

6 

lOSi 

it  flf9Sjnp  11 

10* 

itf  9 

no 

tf.  ■  26 < -MIC S 3 

1 1 1 

S-1-? 

112 

jnp  4 

113 

it  TIC  I.  I  3\ "  0 "  or 

TtCli  13/  ** 

114 

» ’  i  '  nun'  TIC  It  I  3  - 

32  > -WSJ 

1  ;7 

S-l-S 

1 1  * 

ne '<f  ; 

1  *  7 

:•  r3»l*nfTI<  and 

r4sOJ«tO  ‘ 

113 

i*r3  or  Ifr4!  jnB  4 

l:  ? 

b—  "Dif.  Error  #' 

«  J 

120 

fib  "Error" 

FILE  8  (cont) 

12'.:  vo  "Die" 

r*cr:-r* 

122:  :r  t*i  ’.i -4 

124:  tt;  .2?>-N*tSl 

1  t  *  .  C  *  t  ..  4 

lli':  T«23’TJ 

12':  u  T«  J.<e  3 

125:  -t!  i29>*M*CS: 

: _ f-:  i-l-S 

1201  if  T*tl.u«"-'  or  TtC  1.  ’.  J* •■--! T»: 2  ]*Tf 
::i:  t*c  i  ]-r* 

232*  ieniTfwL 

2  3  3!  it  L 3  2 *  J 0  -14 
1  34  i  tor  !»1  10  L 

13*:  :*  T«MK"r  o-  Ttll  •:  3 -••?•!  jiio  -i« 
Hi:  ■  nuniT*:  It  I :  -32'-*MtS  : 

133:  S-l-S 
2S«"  next  I 
13?:  eto  "next  J' 

14i3:  "Other": 

1 4 ; :  ;  en  ■'  US  i  -L 

142:  if  f 1*31 eto  'fill  Other" 

14":  it  tle4ieto  "Special” 

144:  it  4 

14f:  ort  "Nun.  Too  Lone  #">J 

;4<;:  e»b  "Error" 

1J3:  era  “Scar" 

143:  if  L<25.'nP  6 

14t;  nun' U*C 1» 1 3  s -43  1 ytnufl'  Of C  2' 4 I>-4Q-V 
If«:  i‘  V.'*0  and  V<»?9!eto  "Special" 

If’.:  c".  "Op»r.  Error  #"<J 
lf2:  ?sb  "Error" 

If 3:  »to  "Scan' 

If 4:  0-V 

Iff:  *or  1*1  to  26 
Iff.  v  0*»U*C  13!  I*l31tV 
1 f 3:  next  1 

if:;  ;•  "40i*to  "Special” 

If  f:  nurn  (Jf ) -4$4V 

lot?:  i*  7>»0  and  V<*?«etc  "Special" 

16'.:  p-t  "SU  Die  Op  Err  #"tJ 
162:  ?ob  "Error" 

16?:  *to  "Scan" 

1-34:  Special's 
16V.  it  tle5U28*V*V 

Hi:  if  TfP'CTO"  an  a  TM"XEQ '">  eto  “Sto" 

163:  i t  fle3*lijmp  6 
166:  ft!  (174>-»fl*CSl 
.63:  ;•  Tf ■"CTO" < V-128*V 
l*j:  ttl  iVieMftS+i: 

1 *1 :  9-2-S 

1“2:  eto  "next  J" 

it  tie4»lieto  "Sr  Sto" 

134:  i-x 
1  fa.'- 

:~i:  til  1 2? ' -M*C  S  3 

;  — :  rt*  "XEQ"  > ft  i  r30t-n*tsi 

1 3 : :  I en ' Uf  >  rL 

1'?:  *  *  i  •.24e<-L>''(1fCi*l  3 

1:9:  Char": 

1 3 1 :  *or  1*1  to  L 
122:  >8S Y*C 

1:3:  in*  ■  •  B*C3 ''21-M 

;.?4:  :*  L>«C  If  I  :»Cft  «:;onp  5 
1  ?f :  1-  CfC  I  •  I  ]>C*t  mi:m*i-b 
lii:  it  :j*C  !•  1  KCftnlin-leC 
li~:  if  e<*Cijn»  -4 
198:  9-H 
1:-?:  n-r 

1?0:  it  240! anp  4 
1  ?•  1 :  ort  "Bad  Alpha  Lbl  »"tJ 
I ?2:  ffb  "Error" 

I?3:  eto  "Scan" 

1?4:  ■  CC23.'-M$CS+l*n 

l?f:  next  I 
:?f  3-L-2-S 
1?':  no  "next  J" 

;??:  "TP  Cto'i 

1??:  it  ""14  or  Tf»"':E0"*jnP  5 

2 .?o:  fti  .ir~*t.’i-ntci3 

201 :  tf.  '01-MtSM  J 

202:  T-S-S 

2C-3:  eto  "next  J" 

204 :  fi  .20S<-N«S: 

2'?e:  it  Tf*":;EQ"«f t i  .224/«rt«SJ 
2i.it  tt,  0 <-M*C S*1  3 
203’.  fi  i  V  '-Mft  S*2  3 
2981  5*3*$ 

2'.-?:  ••  o  "ne  '  J" 

3 1 0 :  Sto': 

2:-.:  :•  Tf#“ST0" * jne  ? 

212!  it  f  ISIjnP  * 

2 13:  t*i  .  y*42  >*ri$t  S  3 

2.4:  3-1-i 

2lf:  <to  "ne-t  J" 
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FILE  8  (cont) 

216:  "i  MVW*:s: 

>m  v •-«« s*i 3 

21|:  I*;-;. 

21 8:  e*  o  "n*vt  j“ 

Ei1?!  It  T*# “RCL“  ■  )ne  8 
III!  :  *  "  l?!;. in  4 

222:  "i  t 3i*V  <*MtC  S  j 

223:  9»l-»S 

324:  no  “nen  J" 

223:  ft!  t  1 44  ■ -M*t  S  J 
226:  tti  tV-MISM  ] 

44':  2*S*S 

228 :  4to  "next  J" 

-2  9:  ;*  T*#"l6l  '  •  »t  o  'All  Other 
230:  it  f 1 o4# 1 i jno  10 
231:  1*X 
232:  59»v 

233:  ft,  a92f*(l«S3 
234:  ,  0>»M*CS*13 

2??:  l»rvU*'-L 
236:  ■- L»24 1 1 »M*t  i»2  3 

233:  tt,  .  0)4H*cSt33 
228:  i*2-S 
239:  ?t  o  "Char" 

240:  if  V  >1 41  jmo  4 
24 1 :  ft,  ( l  *v '  *M*t  S  3 
242:  1*S*S 
243:  ?to  "next  J” 

244:  tt,  C207  i  »M*t  S  3 
245:  ft,  ■  v  1  »M*C  S»  1  3 

245:  s*2*S 

247:  eto  "next  J" 

249:  -All  Other": 

249:  i-;i 

250:  164-V 

251:  x-e:Y-»C 

232:  int  <  <  3*C>'2."»rt 

233:  it  T*- 1  *t  n  3!  inp  3 

-54 :  Tt;:«n:iM4l-£t 

233:  :•  T*' I *C M 3 • M- i -C 
236:  ,t  E-  *Ci  jno  -4 
23':  0-n 
233:  (1-2 

239:  ,,  2»0;  , <ip  4 

263:  ort  Bad  Inst.  #",J 

261:  «t  "Error* 

262:  4*0  "Sean" 

263:  it  f  ret  IC  2  3>*0I  .'no  4 
264:  tt  ,  in*. ,  IC  2  3 :  )»M*[  S  3 
263:  **i  •  10O0f  rc(.IC233)-M*[S*l3 
266:  3*2'SiJno  3 
267:  ft ,  It  2  3>»M*C  S  3 
268:  S^ltS 

26?:  ;•  1123  64  or  11  2  3/142  or  *l»3i.>iie  4 
27-3:  ort  "Extraneou*  Oo  #",J 
27;:  ?j(j  "Error" 

272:  4*o  "Scan" 

273:  it  IC23  144  or  IC23>159  and  IC23<16£t4to  "n«xt 
274:  ;  t  1123,144  or  K23/191  or  V«0ijno  4 
273:  ort  "Atjjirt*  Operand  #"tj 
276:  4S0  "Error" 

277:  »*o  "Scan" 

273:  •*,  '.V3»MfCS3 
2'9:  SM-S 
260:  "next  J": 

281:  next  J 
232:  "End": 

2S3:  ft,  .l923-n*CS3 
2341  ft,  (0)rh*CS*l 3 
283:  ft,  -  47  '-.HfCS*23 
236:  S*2x« 

237:  ort  “Cone, ltd" 

238:  jf 4  8 
28s:  idf  4.0 
290:  “Error": 

291 :  r0-s 

292:  or*  In*t.  4iw*n  wai'tS* 

293:  ort  "Enter  correct" 

2?4:  on  "Inif.vHo  Line*'" 

293:  e-t  "To  abort  tt  ee" 

2-6:  eno  "880PT1- . jf 

2?':  It  rj« "AgORT" : re" 

2?::  44C-7JC,"  i  "*fl* 

299:  .  ' T  •  »R*C  J  3 

■00:  r*  1  en<  T* '"»!•» T 
JO  1 :  -e* 

302:  end 
*  1 3 1 66 
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INSTRUCTION  mnemonic 

:|  FllE 

1 

AnD  NUMERIC  Values 

1 

j/ 

75.000 

61 

CLFC 

136.000 

2 

•<c« 

77.000 

62 

CLST 

115.000 

i. 

L* 

•1.000 

63 

CLX 

119.000 

4 

J- 

72.000 

64 

COS 

90.000 

• 

:pec 

112.000 

c* 

D-R 

106.000 

i, 

* 

66.000 

66 

DEC 

95.000 

7* 

♦ 

64.000 

67 

DEC 

128.000 

3 

- 

65.000 

63 

2SE 

151.000 

? 

67.000 

69 

EMC 

150.000 

10 

:  x 

96.800 

76 

EtlTEPt 

121.000 

.1 

i 

37.000 

71 

ETX 

35.000 

1  * 

'CLF'EG 

167.203 

~2 

Etx-i 

36.000 

1  c 

tdspo 

167.204 

* 

FACT 

98.000 

14 

75  jpi 

167.205 

74 

PC’ 

173.000 

ir 

-dsf: 

167.206 

->  r 

FC’C 

171.000 

It 

7DSP3 

167.207 

76 

FIX 

156.000 

ir 

7D3P4 

167.200 

77 

PRC 

185.000 

1: 

7B5P? 

167.209 

7? 

PS’ 

172.000 

1 9 

753P6 

167.210 

?t 

PS’C 

170.000 

20 

7BSP7 

167.21 1 

30 

GRPL 

130.000 

^  1 

7BSPS 

167.212 

01 

HMS 

108.600 

75«o« 

167.213 

82 

HMS* 

73. 000 

22 

7BSPI 

167.214 

33 

HMS- 

74.000 

24 

7D62 

167.215 

34 

HR 

109.006 

■?*! 

-ss:: 

1 67. 216 

35 

I  NT 

104.000 

if 

7ENC 

167.217 

86 

ISC 

150.000 

2^ 

7P I 

167.210 

8’ 

LRSTX 

116.000 

22 

7CSEI 

167.219 

86 

LN 

68.666 

2<i 

"CTOI 

167.220 

39 

LM1»X 

101.000 

5e 

"?  i“ 

167.221 

90 

LOG 

86.006 

2i 

7is:I 

167.222 

91 

MEPN 

124.000 

32 

76.'  .c 

167.223 

92 

MOD 

75. 000 

22 

7PPPEC 

167.224 

93 

MPC 

167.193 

24 

7PPST1 

167.225 

44 

OCT 

111.000 

*?* 

7PP:: 

167.226 

9? 

OFF 

141.000 

it 

7PCL3 

167.227 

96 

P-P 

7e.@00 

v- 

7sc: 

167.220 

97 

PI 

114.000 

pbs 

97.000 

96 

PP5, 

167.082 

22 

PCS 

167.665 

99 

PRP 

167.072 

40 

PCCHP 

167. 066 

100 

PPPXIS 

167.073 

41 

ACCOL 

167. 067 

101 

PPBUF 

167.074 

42 

SCO? 

93.000 

102 

PPPLPGS 

167.075 

41 

ACSPEC 

167. 666 

ip: 

PPKEYS 

157.076 

44 

PCX 

167.069 

104 

PPOMPT 

142.000 

4T 

ADv 

143.000 

105 

PRPLOT 

167.878 

4F 

AOFF 

139.000 

106 

PPPLOTP 

167.079 

4? 

PON 

140.000 

10~ 

PPPEC 

167.080 

42 

APCL 

155.600 

1 66 

PPRECX 

167,061 

49 

ASHr 

136.060 

109 

PPSTK 

167.883 

r 

*  V 

RSIN 

92.086 

110 

pp>: 

167.084 

3: 

p?tg 

154.000 

111 

PSE 

137.000 

PTPN 

94.000 

112 

R-D 

107.000 

./•» 

pvre« 

126.000 

113 

P-P 

79.000 

?4 

BEEP 

134.000 

114 

PAD 

129.000 

33 

BLDSPEC 

167.070 

115 

PCL 

144.000 

32 

CP 

169.000 

lie 

PDN 

117.000 

37 

CHS 

84.000 

117 

RDTA 

167. 194 

52 

CU 

112.600 

ne 

PDTAX 

167. 195 

59 

CLP 

135.000 

119 

PECPLOT 

167.005 

20 

CLD 

127.060 

120 

RND 

110.000 

121 

RSUB 

167. 196 

122 

PTN 

133.000 

123 

PT 

116.000 

124 

SCI 

157.000 

125 

SEE" 

125.080 

126 

SF 

166.000 

127 

S1CN 

122.000 

123 

SIN 

69.006 

12’ 

Sl.'PCHR 

167.086 

130 

skpcol 

167.08- 

121 

SORT 

32. OOP 

no¬ 

ST* 

148.006 

ise 

ST* 

146. 000 

134 

ST- 

147.000 

•  it 
.  -V 

ST. 

149.800 

136 

stkplot 

167.086 

137 

STO 

145.000 

138 

STOP 

122.006 

139 

TAN 

91.06P 

140 

TONE 

159.000 

141 

VIEW 

152.000 

142 

WDTA 

167.199 

14? 

UDTAl; 

167.200 

144 

UNDDTA 

166. 193 

14J 

UHESTX 

166. 194 

146 

undlnk 

166. 195 

147 

WNDSCN 

166. 1 57 

1 4.5 

UNDSUB 

166. 196 

14’ 

UNDTST 

166. 198 

1  5P 

WSTS 

167.202 

15: 

X#0’ 

99.00C 

152 

X«v"> 

121.000 

is: 

X'8’ 

102.060 

154 

X'*0’ 

123.000 

153 

X 

70.006 

156 

x<  > 

266.000 

157 

X<  >Y 

113.000 

158 

X<-Y’ 

56.000 

15? 

X»0’ 

103.000 

166 

X*Y’ 

126. eee 

161 

X>0’ 

100. 000 

162 

y  >y*> 

6’. 000 

16: 

.'•7  2 

81.000 

154 

YTX 

33.000 

LOCAL 

LP8EL  AND 

S'SCI. 

FECI  STEF. 

CMPPACTEPS 

1 

A  21 

£ 

e  .22 

c 

C  23 

0 

£ 

D  24 

c 

5 

E  25 

d 

6 

F  26 

« 

7 

G  27 

8 

H 

9 

I 

10 

J 

1  1 

T 

12 

2 

13 

Y 

14 

X 

;«! 

L 

*•$ 

r 

18 

19 

20 


vain*  mp-4ic  ch*«*cte*s 
faND  CH«a*eT£.s  CODE 


!«• 


Vi 

4 ; 

4: 

44 

4  * 

4* 

4 

4: 

4  - 


5t 


5t 

?«• 


4t 

-  it 


'-t 

I.C 
1  00 


,  ooc 
.  000 
,  000 
.  000 
.  00c 
.  000 
.  000 
.  000 
.  000 
.  000 
.  000 
.  00  c 
.000 
.  000 
.  000 
.000 
.  000 
.  000 
.  0  0 
.  000 
.000 
.  000 
•.000 
.000 
.  oCC 
..000 
:  .  000 
.  C  0  c 
' .  00  c 
■ .  000 
■ .  000 
i  •  000 
**  000 
' .  000 
: .  ooo 

L.  000 

; .  00 v 
».  000 
[  ,  I*’  Mi 
:  .  00*. 
000 
: .  OOO 
*  •  0  0  0 
* .  ooo 
: .  ooo 
L.  OOO 

• ,  ooo 
s  •  00:0 
r .  ooo 

:  •  OOO 
‘.OOO 
i .  COO 

' .  ooo 

i.  OOO 
l.  ooo 
■.  000 
ooo 
ooo 
ooo 
ooo 
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file  9 

0:  "renunber  ' : 

t:  ti't  Old  3tart#">0 
2:  er,  "hewStart "  > P 
3:  ■  Inc  recent  “ .  0 

»:  t  or-  I*P*l  to  0-1 
5:  it  :•  *  -  P*C  1 3>*-l i  im>  2 
oi  brt  Error-Overwrite  Inst. "lid* 
7:  n?  '»  I 
S:  tor  1  =  1  to  750 
p*m-»K*u: 

10:  :f  I  -*0.  f  1 1  <-l>tP#m 
1 1 :  ne xt  t 
12:  0-1 *K 

13:  *  or  I=p  to  750  tor  Q 
14:  f+l-K 

15:  If  K  750i ldf  4.0 
16:  if  itf  OCSClUXeijnb  -2 
17:  k*[K3-P$CI3 
18:  ne^t.  1 

1-*:  srt  "Error-#  too  laree" 

20:  for  1*1  to  750 
21 :  K*E  1 3*P*C  I  ] 

22:  next  I 
23:  Id*  4,0 
-256S2 

FILE  12 

0:  ":a.,ebroe": 

::  •Sir.  "lave  on  whicn  f  1 le#<  >15.,  ,  i 

2:  -o*  rl0.Rs.PS. MS 
3:  Id*  4.0 

*14414 

FILE  13 

0:  " Ge *  pro ?" : 

•*  *rr  Get  *  ron  which  f  i  1  *#°  *■  >15; "  i 
2:  st  *  8 

2:  Id*  -:0.ft|.P*,rt* 

4:  Id*  4.0 
*20007 


File  15 


File  14 


C:  Born  lot  :c»  4  0 

l:  ecb  "Bar  Title’.TS 

2:  inn  "Nuiiber  o*  Bytes'5"  >6 

2:  frit  "BecinalClB)  or  Binar»(2)“,D 

4:  Pos": 

5:  enb  'Vert  .Pot.  #<  1-20)  *.  R 
6:  inn  "Hor. Pos. #< 1-15) " • H 
7:  if  t 1 ?0) 4t o  "Plot " 

8:  it  D*18!eto  "Decical" 

?:  "Binary": 

18:  for  1*1  to  8B 
11 :  ant  BCI) 

12:  next  I 

13:  tor  I«36  to  1  By  -1 
14i  6CI3*BCI*23 
15:  next  I 

16:  0+BCi:*8C2J«BC8B*4J 
17:  l*Bt  8B*3  3 
IS:  ft o  “Plot" 

19:  "iecinal": 

28:  tor  I»1  to  B 

21:  tor t  "&yte“»I 

22:  into  tllfti  tM)*MC I  3 

23:  next  1 

24:  for  1*1  to  B 

25:  ltf  frt*C  I  3)*.X 

26:  for  Y*2*8!  to  3*3':-l)  br  -1 

27:  Xmod2*EC Y 3 

28:  iX-Bt  Y  3>-2*X 

29:  next  Y 

38:  next  I 

31 :  0*BC  1  3*612  3*61  3B*4  3 
32:  1 *BC  8B+3  3 
33:  “Plot": 

34:  wrt  705. "IP  000. 080. 1 1 080, 8000" 

35:  sol  0,11.0,8 

36:  csiz  1.5,2,8x11 

37:  if  Z'  16-H)  >B‘,  jnp  2 

36:  brt  "Mo  Roon'i'iste  0ieto  "Pot' 

39:  216*C!beeblwait  400ibeeo!uait  408ibeeo 
40:  ,nt  3b, "Paper  in  Plotter",  lib 
41:  wrt  0,C.C,C,C,C,C,C*C*33,C,C,C,C,C,C,C,C.C 
42:  »en#  l 

43:  bit  .5<H-1>, .5CR-1) 

44:  urt  705*  "PD” 

43:  b«n 

46:  inn  "Position  Good'’ < Yet  o*  no'".Uf 

47:  if  oao<U*>*“N0"  or  ca»<U*  •  »  0!  »•  o 

48:  islt  0. . 35, 1 S Ibl  T* 

49:  bit  . 5< M-l ) , . 5(R-1 ) 

50:  tor  1*1  to  8B*4 

51:  If  BCI3*0,lbll  0,.3,2Stblt  6.  - .  3.  - 1 S -trb  2 
52 :  i r  1 1  8 . .  3 , 2 i  i c  1 1  8,  -  .2,  —  1  *  i c  l'  .  8 1  3 •  3 s  i r.  ’  • 
53:  iblt  .04,8,1 
54:  r,e„t  I 

55:  enb  "Bar  Conblete-ftn.  core  ,  u* 

56:  if  cob'.U*  .'*"YES"  or  oabOJt  ft  o  So'tl 

57:  Id*  4,0 
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«:  iricands " : 

1:  „ir’  6-  Cocnands  available  in  this  Broeran  are:" 

2:  wrt  6, "" 

3:  wr-  6,"  RBOPT-Error  routine  exit  oonnand* 

♦  :  •*•••  6  EPRPlOT-, Generates  a  row  of  user  defined  bar  code" 

3:  wr*  6-  '  at  a  user  defined  location  on  Plotline  surface" 

6:  ■«."  CC'MB ILE-Conbi les  the  41C  Proean  currently  entered" 

7*  w r *  6<"  DRTR-Plots  barcode  rebresentine  Ln*,Lat,and  Lon" 

e:  6-"  DELETE  nn  -Deletes  the  numbered  instruction  t ron  the  broeran " 

9:  wr,  6*  "  EXIT-  Halts  this  oroeran  or  stops  the  nunber  eenerator" 

18:  wrt  6,  GETPR0G-  Retrieves  the  conbiled  code  fron  cassette  robe" 

11:  ur*  6>"  LIST-  Lists  the  entire  41C  broeran  currently  m  nenory” 

12:  wrt  6,'  NUMBER-  Generates  41C  instruction  nunbers-stobbed  by  EXIT" 
12:  wrt  6."  REMUHBEP-  Renunbers  the  41C  Broeran  instruction  nunbers" 

14:  wrt  6,”  PROG-  Generates  the  barcode  for  HP-41C  Proeran" 

15:  wrt  6>"  PROMPT-toads  bronet  subroutine  if  calculator  is  stobbed” 

16:  wrt  6>"  RESET-Clear*  broeran  to  start  over" 

17:  wrt  6.  3ECPR0G-  Generates  the  bar-code  for  a  private  broeran" 

13:  wrt  6,"  3RVEPR0G-  Stores  the  conbiled  code  on  cassette  tape” 

IPs  * ,rf  6,  5CPRTCH-  Erases  trie  entire  41C  broeran" 

20:  wr,  6»  YES'NO-  YES  lists  the  connands  available" 

21:  cr*  6,"" 

22:  wr,  6>"S»ntax  for  instruction  entry:" 

23:  vt  6,"  R  r  Inst  ruct  ion  fornat" 

24:  !'■  6.  n  <4iC  inst ruct ion; " 

25:  ->  6-"  t Blanks  are  del  in  iters) " 

26:  6>-  B  Sbecia!  S.nboli: " 

27:  ,r"  6.  I)  Use  *5.*  instead  of  sieca  sien" 

23:  wr*  6>"  2)  use  *9’  instead  o*  anele  sien' 

2?:  w*«  6>"  3)  use  '#’  instead  o*'  the  no*  etual  sien" 

28:  wrt  6>  4;  Use  sinele  vice  double  euotes" 

2;:  ■-•*•  6*'  C  'Tert  comat:  '  Te-t  Entr, '*  'Note  sinele  tuotes  " 

22:  wrt  6*  "  or  R  *  Text  Entry)*  ,Por  abbendme  te  t ,  ♦ 

33:1  J,  4.0  PIC 

*4?3  r-lo 


